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Foreword 


This supplement to the Guidelines for Maintaining and a Chemical Fingerprinting 
Program has been developed to assist NASA personnel, contractors, and sub-contractors in 
defining the technical aspects and basic concepts which can be used in chemical fingerprinting 
programs. This material is not meant to be totally inclusive to all chemical fingerprinting 
programs; but merely to present current concepts. Each program will be tailored to meet the 
needs of the individual organizations using chemical fingerprinting to improve their quality and 
reliability in the production of aerospace systems. 
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1.0 INSTRUMENTATION 


Modern instrumental analysis laboratories utilize a variety of chemical analysis instrumentation. 
This chapter describes the most commonly used analytical instruments, the principles behind 
the corresponding analytical techniques, and the applications and requirements of these 
techniques. 

1.1. Chromatography 


Chromatography is used to separate the components of a mixture. A chromatographic system 
consists of two mutually immiscible phases: a mobile phase which can be a gas or a liquid, and 
a stationary phase which is either a solid or a liquid supported on a solid within a column. The 
mixture is introduced into the mobile phase which flows through the stationary phase. 
Chemical species within the sample interact with both the mobile and stationary phases. The 
extent of interaction depends upon the chemical and physical properties of each component. 
Each component of the mixture partitions, or distributes itself between the two phases based on 
properties such as polarity, charge, or molecular size. With proper selection of mobile and 
stationary phases, the component species are gradually separated into distinct volumes or bands 
within the mobile phase. Separated components are eluted from the column in order of 
increasing interaction with the stationary phase. A detector placed at the end of the column 
responds to the eluted species, and its signal is plotted as a function of time. The resulting plot, 
called a chromatogram, consists of a series of peaks. The location of a peak along the time axis 
can be used to identify the component, and the area under each peak provides a quantitative 
measure of the component. 


Chromatographic techniques can be broadly classified according to the physical state of the 
mobile phase (Figure 1.1). Gas chromatography (GC) refers to a technique in which the mobile 
phase is a gas, whereas in liquid chromatography (LC) the mobile phase is a liquid. 


Chromatography— 1 


Gas chromatography ■ 


"Gas-liquid (GLC) 


■Gas-solid (GSC) 


Supercritical fluid chromatography (SFC) 


Liquid chromatography 


Liquid-liquid (LLC) 

Liquid-solid (LSC) 
Size-exclusion (SEC) 
Ion-exchange (IEC) 


Figure 1.1. Classification of Chromatographic Techniques 
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Chromatographic methods can be further categorized according to the stationary phase’s 
physical state or its mode of interaction with the sample. When the separation involves 
partitioning between a gaseous mobile phase and a liquid stationary phase, the method is called 
gas-liquid chromatography (GLC). In gas-solid chromatography (GSC), separation is based on 
size exclusion or adsorption of sample components onto the surface of a solid stationary phase. 
Liquid-liquid chromatography (LLC) and liquid-solid chromatography (LSC) refer to the 
analogous methods using liquid mobile phases. Two additional liquid chromatographic 
methods are commonly used. In size exclusion chromatography (SEC), species are separated 
based on molecular size due to differential permeation of a porous stationary phase. In ion- 
exchange chromatography (IEC) ionic species are separated by selective exchange of counter 
ions with an ion-exchange resin. 

Recently, supercritical fluids have been used as chromatographic mobile phases with both 
liquid and solid stationary phases. Supercritical fluids exhibit properties intermediate between 
gases and liquids, and supercritical fluid chromatography (SFC) can be described as a hybrid of 
GC and LC. Although SFC has advantages over both GC and LC in certain applications, it 
cannot replace either method. SFC is infrequently used in industrial laboratories and will not be 
discussed further in this handbook. 

Discussions of the kinetic processes and physical forces that comprise the theoretical basis of 
chromatography are beyond the scope of this handbook. However, several definitions and 
concepts are essential for a practical understanding of chromatographic techniques. 
Chromatographic separations arise from selective retention of components on the stationary 
phase. Retention results from interactions between the component and the stationery phase. 
Four modes of interaction occur in chromatography: partition, adsorption, size exclusion, and 
ion-exchange. In the partition mode, sample components distribute themselves between the 
stationary and mobile phases on the basis of the relative phase solubility of the components. 
Components having differing phase solubilities will spend different amounts of time in the 
stationary phase, and will be eluted from the column separately. Most GC and LC separations 
are based on partitioning. In the adsorption mode, components selectively adsorb onto the 
surface of the stationary phase. At one time, adsorption was the most widely used separation 
mode in both GC and LC. Although adsorption chromatography sometimes suffers form 
irreversible adsorption or peak distortion due to slow desorption, it remains the preferred 
approach in certain applications. In the size exclusion or sieving mode, separation is based on a 
component’s ability to penetrate the pores of the stationary phase. Large molecules that are 
excluded from the pore structure are rapidly eluted from the column, whereas small molecules 
permeate the pore structure and remain in the column longer. Size exclusion is used in both LC 
and GC. Ion-exchange is based on exchange equilibria between ions in the sample solution and 
ions on the surface of an ion-exchange resin. This mode of separation is applicable only to LC. 

Parameters used to describe chromatographic separations can be understood by examining the 
hypothetical chromatogram in Figure 1.2. 


2 



Figure 1.2 Hypothetical Chromatogram 

In this example, a sample composed of Species 1 and 2 entered the system at time zero on the 
chromatogram. Solvents and sample components that do not interact with the stationary phase 
move through the column at the velocity of the mobile phase, and are eluted from the column at 
time t]yj. In this example, species 1 and 2 interact with the stationary phase, with species 2 
interacting more than species 1. The retention time, tR, of each component is defined as the 
time required for the component to elute from the column and be detected. The peak width, W, 
is obtained by drawing tangents to the sides of the peak and measuring the distance between the 
tangents as they intersect the baseline. 

The sample is initially applied to the column as a narrow band on plug. However, as the species 
present in the sample move through the column and are separated, the bands broaden. 
Chromatographic peak shapes are similar to the normal, or Gaussian, curve. 

Band broadening adversely affects the efficiency or separation capability of a chromatographic 
system. Band broadening in chromatography can be attributed to mass-transfer processes. One 
such process is eddy diffusion, which is due to the multitude of pathways a molecule or ion can 
follow through a packed column. Because these pathways differ in length, molecules or ions of 
the same species reach the end of the column at different times. A second process that can lead 
to band broadening is longitudinal diffusion, the migration of a molecule or ion away from the 
center of a band, where its concentration is highest, to regions of low concentration on the 
outskirts of the band. Longitudinal diffusion is more significant in gas chromatography due to 
the relatively high diffusion rates of species through gaseous media. Another process 
responsible for band broadening is stationary phase mass-transfer. For a liquid stationary 
phase, this involves the diffusion of a solute through the liquid stationary phase to the stationary 
phase/mobile phase interface where transfer to the mobile phase occurs. When the stationary 
phase is a solid, the stationary phase mass-transfer is controlled by the rate at which the solute is 
adsorbed onto or desorbed from the surface of the stationary phase. 

Band broadening can often be minimized by judicious selection of experimental parameters. 
Improved separation efficiency can be achieved by decreasing the stationary phase particle size, 
by decreasing the thickness of the immobilized liquid in the case of liquid mobile phases, and 
by decreasing the viscosity of the mobile phase. In gas chromatography, longitudinal diffusion 
can be minimized by lowering the temperature of the column and mobile phase. 
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Three factors control the resolution or quality of a separation: efficiency, capacity, and 

selectivity. All can be calculated directly from the chromatogram. The efficiency of a 
chromatographic column is expressed in terms of N, the number or theoretical plates, or H, the 
height equivalent of a theoretical plate. These parameters are related to each other and to the 
column length, L. 


L = NH 

Column efficiency is increased as the number of theoretical plates becomes larger or as the 
column length increases. The number of theoretical plates can be readily calculated from two 
experimentally measured parameters, tR and W: 


N = 16(^-)2 

Both N and H are used by column manufacturers and by analysts as a measure of a 
chromatographic column’s performance. To compare the efficiencies of two columns, it is 
essential that N or H be determined with the same compound and under the same experimental 
conditions (mobile phase composition and flow rate, temperature, etc.). To compare columns 
of different length, H must be used. 

The selectivity factor, a, is a measure of the relative retention of two components in a mixture. 
Selectivity can be calculated as: 


a = 




where tR 2 is measured for the more strongly retained component. Selectivity depends on the 
nature of the mobile and stationary phases and temperature. 

The capacity factor, k', is a measure of the time a component spends in the stationary phase 
compared to the time it spends in the mobile phase. This parameter is important because it is an 
indication of the column’s ability to retain a solvent. The capacity factor can be calculated from 
the chromatogram: 


k - “m 

Values of k between 1.5 and 4 are desired for adequate retention and reasonable analysis times. 

Efficiency, selectivity, and capacity factors together control the resolution, R, between two 
peaks: 


_ a -i 

R = (— )(VH 
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Resolution can also be calculated directly from the chromatogram: 


R = _ i Gu ^- Vi) 

0.5 (W2 + Wl) 

A resolution of 1.5 corresponds to 0.2% overlap of peak areas and is adequate for most 
separations. Figure 1.3 demonstrates the effect of selectivity, capacity factor, and efficiency on 
the resolution of the separation. 



Many similarities exist between GC and LC. The definitions and basic concepts discussed 
above apply to both techniques. In addition, all instruments used for either GC or LC include 
the following basic components: a mobile phase reservoir and delivery device, a sample 

introduction device, the chromatographic column, a detector, and a readout device. However, 
the types of samples analyzed by each technique, and the exact nature of the individual 
instrumental components required by each technique are different. Because of these 
differences, these techniques will now be discussed separately. 

1.1.1 Gas Chromatography 

Gas chromatography (GC) is used to separate thermally stabile volatile substances. In GC, a 
vaporized sample is carried through a column by an inert carrier (mobile phase) gas. 
Components of the sample are separated due to differences in vapor pressure and affinity for the 
stationary phase. As each component is eluted from the column, its presence is sensed by a 
detector and a response is displayed on a recording device. The major components of a GC 
system are shown in Figure 1 .4. 
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Figure 1.4 Schematic of a Gas Chromatograph 

The samples are usually introduced into the system by direct injection. The sample is injected 
by a microsyringe through a septum into a heated sample port where it is vaporized and carried 
into the column. Use of automatic samplers increase precision and frees the analyst for other 
duties. The carrier gas must be chemically inert and pure. Helium, nitrogen, and hydrogen are 
the most commonly used mobile phases in GC. Because contaminants such as water or oxygen 
can cause deterioration of column or detector performance, purity is essential. 

Two types of columns are commonly used: the packed column and the open tubular or 

capillary column. Packed columns can accommodate much larger sample volumes and are 
usually easier to use. Open tubular columns give much better resolution and are preferred for 
separations of complex mixtures. Packed columns are constructed from glass or metal tubing. 
Column inner diameters range from 1 to 10 millimeters, and lengths range from 2 to 3 meters. 
For gas-liquid chromatography, the columns are packed with a solid support material (125 to 
250 *m diameter) that has been coated with an organic liquid layer immobilized by adsorption 
or chemical bonding. For gas-solid chromatography, the packing may be a porous organic 
polymer or molecular sieves. Open tubular columns are usually constructed from fused silica. 
These columns have an internal diameter of 0.1 to 0.5 millimeter, and a length of 10 to 100 
meters. Wide-bore capillary columns made of glass have an inner diameter of 0.75 mm. The 
inner surface of the column is coated with a liquid stationary phase, 1 to 5 micrometers thick. 

The selection of the stationery phase is a crucial step in method development. Hundreds of 
materials have been proposed as stationary phases for gas-liquid chromatography. The ideal 
stationaiy phase will be thermally stable, chemically inert, and non-volatile. Selection of a 
suitable stationary phase for a specific application is based on the selectivity and polarity of the 
stationary phase. Non-polar stationary phases are used to separate components having 
significantly different boiling points. A stationary phase that will selectively interact with one or 
several of the components should be selected when the mixture contains components having 
similar boiling points. 

The most frequently used liquid stationary phases are the silicone polymers. 
Poly(dimethysiloxane) is a nonpolar stationary phase used for separations of nonpolar 
compounds based on boiling point. By replacing some of the methyl groups with phenyl 
groups, slightly polar phases capable of separating olefins, aromatics, and other unsaturated 
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species are obtained. Replacement of the silicone’s methyl groups with even more polar 
functionalities (such as cyanopropyl or trifluoropropyl) gives a more polar or more selective 
stationary phase. Most chromatographic supply catalogues offer excellent advice on the 
selection of stationary phases. Some common stationary phases, and their applications, are 
listed in Table 1.5. 

Solid materials used as GC stationary phases include the molecular sieves and porous polymers. 
Molecular sieves are porous alkali metal aluminosilicates. Pore size is uniform and depends on 
the cation present. Porous polymer packings are made from styrene divinyl benzene 
copolymers. Molecules smaller than the pore dimension penetrate the particles and are 
adsorbed. Therefore, separation is based on molecular size and shape. Gases and low boiling 
point liquids can be separated on these solid stationary phases. 

Because column temperature must be controlled to within • 1 ° C for precise work, GC columns 
are coiled and housed in a thermostated oven. Column temperature is selected based on the 
boiling range of the sample. If the sample components boil over a narrow range, it may be 
possible to achieve the separation in a single isothermal run at a temperature near the average 
boiling point of the sample components. If the sample boils over a broad range, it may be 
necessary to use temperature programming, i.e., to continuously or incrementally increase the 
column temperature during the separation. The column temperature is initially set below the 
boiling point of the lowest boiling component. The final temperature is near the boiling point 
of the highest boiling component (but within the thermal limit of the stationary phase). Thermal 
programming shortens the retention time of the later eluting components, decreasing both the 
band broadening of the later peaks and the total analysis time. 

The detector senses the presence of the separated components as they elute from the column. 
Dozens of detectors have been used for GC. The five types described below are the most 
frequently used. Figure 1.5 illustrates the range over which the most popular detectors are used, 
while a guide for detector selection is given in Figure 1.31 on page 66. 


Linear Response Range 



Figure 1.5. Useful ranges for Gas Chromatographic Detectors 
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Thermal conductivity detectors (TCD) display universal response and are rugged, relatively 
inexpensive, and nondestructive of sample. The TCD is based on changes in the thermal 
conductivity of the gas stream emerging from the GC column. The sensing element of the TCD 
consists of a metal block container and a filament that is electrically heated. When gas flows 
over the filament, heat is transported from the filament to the metal block. Heat loss from the 
filament results in decreased temperature and electrical resistance. Because the thermal 
conductivities of helium and hydrogen are six to ten times greater than those of most organic 
compounds, even small levels of organics in the column effluent cause large decrease changes 
in thermal conductivity and a marked increase in filament temperature and resistance. A 
change in filament resistance signals the emergence of a sample component from the column. 

Flame ionization detectors (FID) display nearly universal response, high sensitivity, wide linear 
response range, and excellent reliability. Compared to the TCD, the FID is approximately 1000 
times more sensitive, but also more complicated, more expensive, and destructive. In the FID, 
combustible sample components are burned in a hydrogen/air flame, generating free ions and 
electrons. The flame in the FID is located between two electrodes. Ions and free electrons 
present in the pure flame give rise to a small current when a potential is applied across the 
electrodes. As a carrier gas containing combustible sample components passes through the 
flame, additional ions and free electrons are generated and the current markedly increases. 
Response is proportional to the number of reduced carbons in the sample component. Oxidized 
carbons (present in functional groups such as carbonyl, alcohol, carboxylic acid, ether and their 
sulfur analogs) produce little or no response. Because the FID does not respond at all to water 
or to the permanent gases (N 2 , O 2 , CO, CO 2 , etc.), it is ideal for trace analysis in aqueous 
solutions and air samples. 

Thermionic emission detectors (TED) respond only to compounds that contain nitrogen or 
phosphorus. The TED is similar to the FID except that the flame temperature and electrode 
polarity are optimized to enhance ionization of nitrogen and phosphorus containing compounds, 
and to suppress ionization of other compounds. The TED is 500 times more sensitive for 
nitrogen, and 50 times more sensitive for phosphorus than the FID. The TED is frequently 
referred to as the nitrogen-phosphorus detector (NPD). 

Electron capture detectors (ECD) respond selectively to molecules containing electronegative 
functional groups such as halogens, nitrate, nitrite, and peroxide. The ECD also responds to 
unsaturated compounds such as polynuclear aromatics. In the ECD, the chromatographic 
effluent passes between two polarized electrodes, one of which is coated with a radioisotope 
that emits beta particles. The beta particles bombard the carrier gas, producing a burst of 
electrons and current flow between the electrodes. The presence of electron-capturing species 
is detected as a decrease in current. 

Mass spectrometers are used as highly specific and highly sensitive GC detectors. Combined 
gas chromatography/mass spectrometry (GC/MS) can be used to quantitate and identify the 
components of a mixture. Mass spectrometers operate under high vacuum. The gas flow rates 
from capillary columns is usually low enough to permit direct connection between the column 
and the mass spectrometer. With packed columns, an interface system must be used to remove 
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most of the carrier gas. Two frequently used interfaces are the jet separator and the membrane 
separator. The jet separator takes advantage of the faster diffusion rate of the carrier gas 
compared to analyte. When column effluent is forced through a fine nozzle into a vacuum 
chamber it rapidly expands. The carrier gas expands more rapidly than the analyte. An orifice 
aligned with the nozzle collects the core of the effluent stream which is now enriched in 
analyte. The membrane separator is based on differences between the abilities of the carrier gas 
and analyte molecules to permeate a silicone membrane separating the column effluent from the 
mass spectrometer. Organic analyte molecules pass more readily through the membrane and 
into the mass spectrometer. 

Detectors on mass spectrometers can acquire and display data in several ways. All of the ion 
currents can be summed and plotted as a function of time to give a total ion current 
chromatogram. If the analyst is interested in a particular peak the mass spectrum acquired at the 
time the peak passes through the detector can be plotted. The spectrum can then be used to 
identify the separated component. 

In selective ion monitoring (SIM), the instrumental parameters are adjusted to detect only those 
ions associated with a particular compound or class of compounds. For example, if a- 
Methylstyrene were the only analyte of interest, only the 118 and 117 m/z ions would be 
monitored. SIM gives enhanced sensitivity because the background signal due to other ions is 
filtered out and because more time can be spent collecting data for the selected ions. 

1.1.2 Liquid Chromatography 

Liquid chromatography is used to separate mixtures of high molecular weight polyfunctional 
materials, polymers, thermally unstable compounds, and ionic species. Unlike GC, LC is not 
limited to thermally stable, volatile materials. LC and GC are complementary techniques. 

Liquid chromatography was originally performed in large glass columns up to 5 m long and 5 
cm diameter. To minimize mobile phase flow rates, the stationary phase particle diameters as 
large as 200 *m were used. Even under these conditions, separations required several hours. 
Attempts to shorten the analysis time by increasing the mobile phase flow rate resulted in 
decreased separation efficiency. In the 1960s it became possible to produce highly efficient 
stationary phase packings with particle diameters near 10 *m. These small diameter packings 
lead to the development of new instrumentation capable of operating at higher pressures. This 
new technology was called HPLC, high performance liquid chromatography. Except for 
preparative applications, HPLC has replaced classical glass-column LC. 

Figure 1.6 shows the major components of a typical HPLC instrument. The components 
perform the following functions: delivery of the mobile phase, sample introduction, separation 
of the mixtures components, and detection of each component. 

A liquid chromatographic separation can be performed by isocratic elution, in which a single 
mobile phase is used, or by gradient elution, in which two or more solvents are used in varying 
ratios throughout the run. For isocratic elution, the minimal mobile phase delivery system 
consists of a mobile phase reservoir and a pump. Gradient elution requires a reservoir and 
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(usually) a pump for each solvent, and a gradient mixing device. 



frsadout^ 
.device j 


Figure 1.6 Functional Schematic of a Liquid Chromatograph 

Solvent reservoirs are constructed of glass or stainless steel. The original solvent bottle is often 
adequate. The mobile phase itself must be free of particulates which may cause blockages and 
dissolved gases which may be released as bubbles in the detector or in the pump check valves. 
Particulates are removed by in-line filters placed between the reservoir and the pump. Gases 
may be removed by sparging, sonicating, vacuum pumping, or by heating and stirring. 

Pumps used for HPLC must be able to operate at pressures up to 6000 psi and to produce 
reproducible and constant flow rates of 0.1 to 10 mL/min. Two types of pumps are used: 
reciprocating piston and positive displacement. Reciprocating piston pumps are the most 
popular. In a piston pump, a small (30 or 400 »L) cylindrical chamber is alternately filled and 
emptied by the back-and-forth motion of a motor-driven piston. Flow direction is controlled by 
means of check valves. On the backward stroke, the piston pulls solvent from the external 
mobile phase reservoir. On the forward stroke, mobile phase is pumped to the column. 
Because no solvent is pumped to the column on the backward stroke, the flow is pulsed, and a 
pulse damping system should be used. Dual-head (or triple head) pumps consist of two (or 
three) pistons mechanically coupled so that pumping and filling occur simultaneously. Flow 
pulsations are greatly reduced, but not eliminated, in the dual head and triple-head pumps. 
Advantages of reciprocating pumps include unrestricted operating time due to the use of an 
external reservoir and rapid solvent change over due to the pump's small internal volume. 
Positive displacement, or syringe, pumps consist of a large (250 to 500 mL) solvent reservoir 
equipped with a plunger. A stepping motor actuates the plunger through a screw-driven 
mechanism. Displacement pumps produce a pulse-free flow. However, it suffers from limited 
solvent capacity, requiring periodic shut-down of the system for refilling. 

Two types of gradient mixing devices are used: low-pressure and high-pressure. In a low- 
pressure system, the gradient is formed ahead of the high pressure pump. A system of 
proportioning valves accurately measure and deliver up to four solvents to a low volume mixing 
chamber. In a high pressure system, mixing of the solvents occurs after the pumps. High 
pressure systems are more costly because they require a separate pump for each solvent used in 
the gradient. However, they provide more precise control over the gradient composition. In 
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both types of gradient systems, changes in the mobile phase composition are controlled by a 
gradient programmer. Both continuous and step gradients are used. Because the mixing 
process can generate heat which can result in gas evolution thorough degassing of each mobile 
phase is essential. 

The most widely used sample introduction devices are sampling valves. These devices permit 
sample introduction at high pressure with minimal interruption of mobile phase flow. Typical 
valves can deliver 10 to 500 *L with a few tenths percent precision. Automatic samplers 
available for HPLC allow unattended operation of the equipment. 

HPLC columns are constructed of smooth-bore stainless steel tubing or glass-lined metal 
tubing. Several types of analytical columns are used: standard, short, and narrow bore. 

Standard columns are 20 to 30 cm in length with an inner diameter of 4 to 5 mm and stationary 
phase particle size of 3 to 10 *m. Short columns that are 3 to 6 cm long and packed with 3 »m 
particles give good separations with increased sample throughput and minimal solvent 
consumption. These columns are useful for relatively easy separations (N < 4000) when speed 
is essential as in quality control applications. Narrow-bore (or microbore) columns have small 
internal diameters, usually 1 to 2 mm. Narrow-bore columns offer several advantages including 
decreased solvent consumption and increased detector response. Disadvantages include a 
requirement of special low volume injection valves, detector cells, and pumps. 

Guard columns are sometimes placed between the injection valve and the analytical column. 
Guard columns are short and packed with a stationary phase similar to that used in the 
analytical column, but of larger particle size. Guard columns protect the analytical column by 
trapping particulates and sample components that would be permanently retained on the 
analytical column. 

As mentioned before, there are four modes of liquid chromatography: liquid-liquid (partition), 
liquid-solid (adsorption), ion exchange, and size exclusion. Each mode is applicable to 
different types of samples and each mode makes use of different types of stationary phases. 
Liquid chromatographic techniques are also categorized on the basis of the relative polarities of 
the stationary and mobile phases. In normal-phase chromatography, a polar stationary phase is 
used with a nonpolar mobile phase. In reverse-phase chromatography, the stationary phase is 
less polar than the mobile phase. 

Column packings used for partition chromatography are of two types. In the liquid-coated type, 
the liquid stationary phase is held on the solid support particles by adsorption. In bonded-phase 
packings, the stationary phase is covalently bonded to the support particles. Bonded-phase 
packings have almost completely replaced the liquid-coating packings due to their greater 
stability. Bonded phase packings are based on rigid silica particles. The organic stationary 
phase is bonded to the silica surface through a siloxane linkage: 


— Si — O — 
I 



R 


where R is an organic group. The polarity and selectivity of the stationary phase can be 
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controlled by varying the organic group. In the most popular bonded phase packings, the alkyl 
group is an r|-ocytl (Cg), or an T|-octadecyl (C\$) group. Both are nonpolar, however, the size 
of the alkyl group affects retention. The Cg packing will accomplish the separation faster, but 
the Ci g phase can be used with larger sample sizes. Bonded phases containing phenyl groups 
are also nonpolar, but more selectivity interact with aromatic and unsaturated hydrocarbons. 
Phases that contain cyano groups are moderately polar and are used to separate ethers, esters, 
ketones, aldehydes, and nitro compounds. Amino alkyl bonded phase are highly polar, and are 
useful in separations involving alcohols, phenols, carbohydrates, and amines. 

Column packings used for liquid-solid chromatography are composed of silica, alumina, or 
carbon. Silica packings are the most frequently used. LSC is primarily used separate nonpolar, 
water-soluble molecules. LSC can also be used to separate isomeric mixtures. 

Size exclusion packings are composed of cross-linked polymers or porous glass or silica. 
Styrene-divinylbenzene copolymers are frequently used. 

In packings used for ion chromatography, a monolayer of small (0.1 to 0.3 *m) polymeric 
beads are aminated or sulfonated and electrostatically bonded to a relatively large (10 to 30 *m) 
polymeric or glass bead. The aminated resins are used to separate anions, while the sulfonated 
material are used for cations. 

Liquid chromatography has no detector as sensitive and as universally applicable as the FID 
and TCD detectors used in gas chromatography. The most commonly used HPLC detectors and 
their properties are listed in Table 1.1. 


Detector 

Type 

Sensitivity(g/mL) 

.Linearity 

UV- Visible absorption 

selective 

io-io 

10 5 

Refractive index 

universal 

io - 7 

10 4 

Fluorometric 

selective 

io- n 

10 5 

Amperometric 

selective 

IO' 12 

10 5 

Conductometric 

selective 

IO' 8 

10 3 


Table 1.1 HPLC Detectors 

UV-visible absorption detectors are based on the absorption of ultraviolet or visible radiation by 
components as they pass out of the column. Three types of UV-visible detectors are available: 
fixed-wavelength, variable-wavelength, and photodiode-array. Fixed-wavelength detectors use 
discrete light sources which emit light at several wavelengths. Each wavelength can be selected 
by use of narrow bandpass filters. A fixed-wavelength detector with a mercury lamp as the 
source has useful emission lines at 254, 280, 313, 334, and 365 nm. These detectors are simple, 
stable, sensitive, and inexpensive. However, only the solutes that absorb at the given 
wavelengths will be detected. Variable wavelength detectors use continuous light sources in 
combination with a monochromator. These detectors offer unlimited selection of UV and 
visible wavelengths. Photodiode array detectors consist of a continuous light source, a 
monochromator, and silicon diode array detector. These detectors simultaneously monitor all 


12 




wavelengths. An entire spectrum can be collected and stored as each separated species passes 
out of the column. Photodiode array detectors are especially useful during the development of a 
new LC method because they help the analyst identify the optimal detector wavelength. UV 
visible detectors are selective in that they respond only to species that absorb ultraviolet or 
visible radiation. Detectable species include the unsaturated hydrocarbons, including aromatics, 
and compounds that contain atoms such as nitrogen, oxygen, sulfur, and halogen. 

Refractive index (RI) detectors are based on difference in refractive index between the mobile 
phase (reference) and the column eluent. Because most substances differ in refractive index, 
the refractive index detector is universal in response. However, this detector is not very 
sensitive and is seldom useful for trace-level components. Its universal response usually 
precludes its use in gradient elution procedures where the changing mobile phase composition 
results in drifting baselines. Refractive index detectors are extremely sensitive to temperature 
changes and flow rate fluctuations. 

Fluorometric detectors respond to the fluorescent emission of molecules that have been 
electronically excited by a suitable source. Fluorometric detectors are more selective and more 
sensitive than UV-visible absorption detectors. The detector's selectivity is due to the fact that 
not all molecule fluoresce. Compounds that can be detected include pollutants such as the 
polynuclear aromatics and biologically significant species such as vitamins, alkaloids, and 
catecholamines. The extreme sensitivity of fluorometric detectors results from a low 
background level of fluorescence. This detector can be used with gradient elution. 

Amperometric detectors are based on the measurement of current flow as analytes in the eluent 
undergo oxidation or reduction at an electrode. Typical amperometric detector cells have a 
three-electrode arrangement consisting of a working electrode (which detects the analyte), an 
inert reference electrode, and an auxiliary electrode. A potential is applied between the working 
and reference electrodes, while current flows between the working and auxiliary electrodes. 
The working electrode composition is important because it determines the range of potentials 
that can be applied and, therefore, the species that can be detected. Dual electrode detector cells 
have two working electrodes which can be placed in series or in parallel with the flowing 
eluent. In the series configuration, the upstream electrode produces an electrochemically active 
product that can be detected at the downstream electrode. This approach enhances the 
selectivity of the detector, and also enhances sensitivity by reducing background noise due to 
mobile phase electrolysis. In the parallel configuration, the two working electrodes are held at 
different potentials. The ratio of currents measured at the two electrodes provides an indication 
of the peaks purity and identity. Amperometric detectors are more sensitive and more selective 
than UV-visible and RI detectors. Amperometric detection is widely applicable. However, 
amperometric detectors are not widely used due to several disadvantages. Amperometric 
detectors adversely respond to fluctuations in eluent flow rate. Many compounds adsorb onto 
the electrodes, requiring frequent and time-consuming cleaning and recalibration. Operation in 
the reductive mode may require exclusion of oxygen from the system and the use of mercury 
electrodes. 

Conductivity detectors measure the ability of the eluent to carry an electric current under the 
influence of a potential gradient. These detectors respond to species that form ions in solution 
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(electrolytes). Conductivity detectors are used with ion chromatography. 

The mobile phases used in ion chromatography are ionic solutions which produce a high 
background conductivity. To detect analytes, the eluent conductivity must be suppressed. 
Membrane suppressors are used for anion (and cation) exchange separations. The eluent is 
passed over one side of a cation (or anion) exchange membrane. The membrane is continuously 
regenerated by flowing an acidic (or basic) solution over the other side of the membrane. 

Mass spectrometer detectors provide both structural and quantitative information. Interfacing 
liquid chromatography with mass spectrometry has been difficult because of the mismatch 
between the large mass flow rates used in liquid chromatography and the vacuum requirements 
of mass spectrometry. Several interfaces have been developed. In the moving-belt interface, 
the column effluent is deposited on a continuous moving belt. The belt moves through a heated 
chamber, where the solvent is evaporated, and into the ion source of the spectrometer. In the 
thermospray interface, the column effluent is passed through a heated capillary tube to produce 
an aerosol of solvent vapor and analyte molecules. When polar mobile phases containing a salt 
(such as ammonium acetate) are used, the analyte can be ionized through charge exchange with 
the salt. When nonpolar or weakly ionizable mobile phases are used, an electron beam is used 
to achieve ionization. 
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1.2. Spectroscopy and Spectrometry 

This section introduces the fundamental principles of spectroscopy, and 
describes specific instrumental methods based on these principles. Absorption, 

emission, or scattering of electromagnetic radiation alters the energy state of the 
interacting atom or molecule. Because each chemical species has characteristic energy 
states, a species interacts only with a particular region or energy range of the 
electromagnetic spectrum. The energy at which interaction occurs can be used to 
identify the interacting species, and the intensity of the interaction can be used to 
quantify its concentration. 

Spectroscopic methods are based on absorption, emission, or scattering of elec- 
tromagnetic radiation by matter. Electromagnetic radiation is a form of energy that is 
transmitted through space at the speed of light (3x10^ m/s). Electromagnetic radiation 
can be described in terms of a wave model. An electromagnetic wave consists of 
electrical and magnetic field components which oscillate in planes perpendicular to each 
other and to the direction of wave propagation. The electrical field component of 
electromagnetic radiation is responsible for most of its interactions with matter. Figure 
1.7 depicts the electrical field component and illustrates several parameters used to 
characterize electromagnetic radiation. 



Figure 1.7. Electromagnetic Wave 

The wavelength, 1, is the distance between successive maxima or minima of either the 
electrical or magnetic component. The frequency, u, is the number of waves that pass a 
fixed point in a unit of time, usually a second. Frequency is determined by the source of 
the radiation and remains unchanged by propagation of the wave through matter. The 
velocity of propagation, v, is the rate at which the wave passes through a medium. The 
velocity of electromagnetic radiation in a vacuum, c, is 3 x 10^ m/s. Due to interactions 
between the electric field component and matter, electromagnetic radiation is propagated 
or transmitted through matter at velocities less than c. The ratio of the speed of light in a 
vacuum to the speed of light in a medium is the refractive index h of the medium. 
Because the frequency is invariant, the wavelength must also decrease as radiation 
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passes from a vacuum to another medium. The wavenumber, , in units of cm'l is the 
number of wave crests that occur per centimeter. Wavenumbers are often used instead 
of frequency, and can be calculated as: 


v = \j X or v = — 
c 

The radiant power, P, is the amount of energy transmitted per second and is proportional 
to the square of the wave amplitude, A. 

Refraction, diffraction, and interference are phenomena that can readily be 
explained by the wave model. However, absorption and emission of electromagnetic 
radiation by matter can be described only by treating radiation as a stream of discrete 
particles or quanta of energy known as photons. The energy of a photon depends upon 
the frequency of the radiation, and is given by: 

E = hi) 

where E is in joules, and h is Plank's constant (6.62xl0'34 J sec). 

The quantum model is needed to describe photoionization, the emission of 
electrons from the surface of a solid when a sufficiently energetic radiation impinges on 
the surface. The energy of the emitted electrons is related to the frequency of incident 
radiation: 


E = ho - w 


where the work function, w, of the solid is the work required to remove an electron. The 
number of emitted electrons is dependent on the number of impinging quanta of 
radiation having a certain minimum energy. 

Figure 1.8 on the next page depicts the electromagnetic spectrum, the broad 
range of radiations that extend from gamma rays to radio waves. The various types of 
electromagnetic radiation differ in frequency, wavelength, and nature of interaction 
with chemical species. 
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Figure 1.8. Electromagnetic Spectrum 

Low-energy radio waves cause reorientation of nuclear spin states in materials placed in 
a magnetic field. Photons in the microwave region cause changes in rotational energy 
states of molecules. Absorption of infrared radiation results in changes in both 
vibrational and rotational energy states of molecules and complex ions. Absorption of 
visible or ultraviolet radiation changes the energy states of outer shell (valence) 
electrons. X-ray absorption results in ejection of inner shell (core) electrons (the 
photoelectron effect). These interactions of electromagnetic radiation with chemical 
species will be described in more detail in the following pages. 

When electromagnetic radiation passes through a sample of matter, select 
frequencies may be transferred to the sample's atoms and molecules by the process of 
absorption. As a result, the absorbing species are promoted from a low energy state to a 
higher energy state, or excited state. Most chemical species at room temperature are in 
the lowest energy state, a ground state. Absorption, then, usually involves a transition 
from the ground state to an excited state. An atom or molecule in an excited state may 
return to the ground state by emission, the release of energy as radiation. 

According to quantum theory, atoms and molecules exist only in a limited 
number of discrete potential energy levels. The energy of the impinging photon must 
match the energy difference between the ground state and an excited state of the 
absorbing particle for absorption to occur. Similarly, energy lost by emission of 
radiation must match the energy difference between an excited state and the ground 
state, or between two excited states. Because these energy differences are unique for 
each atom or molecule, the energies (frequencies) at which a species absorbs or emits 
radiation can be used to identify the species. 

Absorption or emission of radiation is accompanied by transition of electrons 
between fixed energy levels. Ultraviolet and visible radiation are sufficiently energetic 
to cause transitions of outer shell or valence electrons. Absorption or emission of x-rays 
leads to transition of inner shell, or core, electrons. 

Molecular spectra are more complex than atomic spectra due to an increased 
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number of energy states available. The potential energy of a molecule is the sum of the 
electronic, vibrational, and rotational energies. Normally, several rotational energy 
states exist for each vibrational energy state, and several vibrational states exist for each 
electronic state. The schematic energy level diagram in Figure 1.9 shows some of the 
electronic and vibrational states of a molecule. Lines labelled E n represent electronic 
energy states. Eq refers to the electronic ground state, Ej and E 2 refer to first and 
second electronic excited states. Several vibrational energy levels (labelled vq ? to V3) 
are pictured for each electronic state. The energy differences between adjacent 
electronic states is 10 to 100 times greater than the energy differences between adjacent 
vibrational states. The figure also illustrates two analytically useful processes, 
absorption and fluorescence. 
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Figure 1.9. Molecular Energy Levels 
1.2.1 Atomic Spectroscopy 

Atomic spectroscopy is used to determine the concentration of a particular 
element in a sample, regardless of the chemical environment and oxidation state of the 
element. Atomic spectroscopic methods are listed in Table 1.2. In each method, the 
molecular components of the sample are converted to free atoms by the process of 
atomization. In atomic emission spectroscopy (AES), absorption of additional thermal 
energy from the atomization device transforms the free atoms to excited electronic 
states. As the excited atoms return to the ground state, they emit ultraviolet or visible 
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radiation at wavelengths characteristic of the atoms present in the sample. The intensity 
of the emitted radiation is measured, and is the basis of the analytical determination. In 
atomic absorption spectroscopy (AAS), the atoms are transformed to excited states by 
absorption of radiant energy from an external ultraviolet/visible light source. The 
analytical determination is based on the amount of radiant energy absorbed. In atomic 
fluorescence spectroscopy (AFS), the atoms are excited by a radiation source placed at 
90° to the optical axis of the spectrometer. Quantitation is based on the intensity of 
radiation emitted as fluorescence. The basic components of instruments used for AES, 
AAS, and AFS are compared in Figure 1.10. All instruments used for atomic 
spectroscopy include an atomization device, a monochromator to resolve the emitted or 
transmitted light into its component wavelengths, and a detector for ultraviolet/visible 
radiation. Atomic spectroscopic instruments differ primarily in atomizer type, the 
absence or presence of an external light source, and orientation of the components. 
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Figure 1.10. Instruments for Atomic Spectroscopy 
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Atomic spectroscopic methods are frequently classified according to atomizer 
type. The most commonly used atomizers are the combustion flame, the electrothermal 
analyzer, and the plasma. Flame atomizers are used for AES, AAS, and AFS. 
Electrothermal atomizers are primarily used for AAS and AFS, whereas plasma sources 
are used for AES and AFS. 


Techniaue ( Abbreviation! 

Atomization Source (Temperature^ 

Flame atomic absorption spectroscopy 
(FAAS) 

Flame (1700-3200° C) 

Electrothermal atomic absorption 
spectroscopy 

Furnace (1200-3000° C) 

Flame emission spectroscopy (FES) 

Flame (1700-3200° C) 

Inductively coupled plasma atomic 
emission spectroscopy (ICP/AES) 

Argon plasma (6000-8000° C 

Direct current argon plasma 
spectroscopy (DCP) 

Argon plasma (6000-10000° C) 

Arc-source emission spectroscopy 

Arc plasma (4000-6000° C) 

Spark-source emission spectroscopy 

Spark plasma 

Atomic fluorescence spectroscopy (AFS) 

Flame (1700-3200° C) 

Electrothermal atomic fluorescence 
spectroscopy 

Furnace (1200-3000° C) 

Inductively coupled plasma atomic 
fluorescence spectroscopy (ICP/AFS) 

Plasma (6000-8000° C) 


Table 1.2 Atomic Spectroscopy Methods 

In the following section, the atomic spectroscopic techniques most commonly 
used will be discussed. These techniques are: flame atomic absorption spectroscopy, 
electrothermal atomic absorption spectroscopy, and inductively coupled plasma/atomic 
emission spectroscopy. 

1.2.2 Flame and Electrothermal Atomic Absorption Spectroscopy 

Atomic absorption spectroscopy is a sensitive technique for the quantitative 
determination of approximately seventy metallic or metalloid elements in solution 
matrices. The basic instrumental components were shown above in Figure 1.10. 


20 






Electrothermal and flame atomizers are used in atomic absorption instruments. 
Flame atomizers consist of a nebulizer and a burner. The nebulizer transforms the liquid 
sample into an aerosol which is introduced into the burner. Figure 1.11 illustrates the 
two main types of burners: the total consumption (turbulent flow) and laminar flow 
burners. 



In the total consumption burner, the nebulizer and burner are combined into a 
single unit. The oxidant flow around the sample capillary tip draws the sample up the 
capillary and into the burner. Fuel gas mixes with the oxidant and the sample, and helps 
to break up the sample. The flame forms at the top of the burner. Total consumption 
burners offer several advantages over laminar flow burners: (1) the entire sample is 
aspirated into the flame, eliminating error due to loss of nonvolatile components, (2) no 
possibility of flashback or explosion exists (see discussion of laminar flow burner, 
below), and (3) the burner is inexpensive and easy to maintain. Disadvantages include: 
(1) vaporization and atomization efficiency are low, (2) operators must take special 
precautions to prevent clogging of the tip, (3) short flame path length results in 
decreased signal in AAS, and (4) total consumption burners are very noisy, both 
electronically and aurally due to turbulence. 

In the laminar flow burner, the sample is drawn into a mixing chamber and 
nebulized by the flow of oxidant across the capillary tip. A series of baffles in the 
mixing chamber remove larger droplets from the sample stream. The remaining aerosol 
mixes with additional oxidant and fuel, and passes into the burner head and the flame. 
Laminar flow burners offer several advantages: (1) sensitivity is greater due to 

relatively long path length (5-10 cm), (2) burner is quiet, and (3) because larger drops 
are eliminated in the premix chamber, laminar flow burners seldom clog. Disadvantages 
are: (1) in samples containing more than one solvent, the more volatile components 
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may be preferentially vaporized in the mixing chamber, while less volatile components 
may drain off and not reach the flame (2) the mixing chamber contains an explosive 
mixture which can be ignited by a flashback, and (3) most of the sample goes down the 
drain. 


Conversion of the sample to free atoms in the flame involves a sequence of 
events. The sample enters the flame as a solution aerosol. The solvent evaporates or, in 
the case of an organic solvent, burns, leaving behind a solid aerosol. The solid particles 
undergo volatilization, and then dissociation to form free atoms. These atoms are then 
excited by radiant energy from a uv-visible light source. Atomization efficiency is the 
efficiency with which the flame produces atoms by this sequence of events. Several 
processes can decrease atomization efficiency and, therefore, the analytical response. If 
the sample drops are too large, they may pass through the flame without completely 
evaporating. Droplet size can change significantly with a change in solvent due to 
viscosity differences. The rate at which the sample is introduced into the flame also 
effects atomization efficiency because solvent evaporation requires energy and lowers 
flame temperature. 

Flame atomization surpasses other atomization methods in terms of 
reproducibility. However, sampling efficiency of the flame is low because much of the 
sample is discarded (in case of the laminar flow burner) or incompletely atomized (in the 
case of the total consumption burner), and the residence time of the analyte atoms in the 
optical path of the flame is short. The sampling efficiency and, therefore, sensitivity of 
other atomization methods is better. 

Electrothermal atomizers provide high sensitivity because the entire sample is 
atomized quickly and the residence time in the optical path is on the order of seconds. 
Compared to flame atomization, electrothermal atomization enhances sensitivity by 
factor of 100 to 4000. Electrothermal atomizers can accommodate very small sample 
volumes (0.5 - 100 *L) and solid samples. Electrothermal atomizers are electronically 
less noisy than flames. However, the precision of electrothermal methods, typically 5 to 
10%, compares unfavorably with the 1 to 2% precision obtained with the flame 
methods. 

An excitation source that emits light having an energy equal to the difference in 
energies between the ground state and excited state of the element being analyzed is 
used. Hollow cathode lamps are the most common excitation sources for AAS. Hollow 
cathode lamps consist of a cylindrical cathode constructed of the element being 
analyzed, and a tungsten anode, both sealed in a glass tube filled with an inert gas. 
When a sufficiently large potential is applied across the electrodes, the inert gas ionizes, 
forming highly energetic cations. As these cations strike the cathode's surface, atoms 
on the surface are dislodged. Some of the metal atoms are in excited electronic states 
and emit lines of radiation characteristic of the cathode element as they return to the 
ground state. Usually, a different lamp must be used for each element, although some 
multielement lamps are available. The source emission is directed through the atomized 
sample in the flame or furnace. The light not absorbed by the sample passes through to 
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the monochromator and the detector. The absorbance is calculated from the intensities 
of light detected with (I) and without (I 0 ) the analyte in the flame: 

A = log (VO 

To absorb the source radiation, an atom in the flame must be of the same element used 
in the source lamp. Therefore, the absorbance is element-specific. The absorbance is 
also directly proportional to the analyte's concentration. 

1.2.3 Inductively Coupled Plasma/Atomic Emission Spectroscopy 

Inductively coupled plasma/atomic emission spectroscopy is used for the 
qualitative and quantitative analysis for over seventy elements in solution. ICP/AES is 
capable of multi-element analysis, performed in either a simultaneous or rapid sequential 
mode. This technique is used for determination of major, minor and trace level 
elements. 

In ICP/AES the sample is atomized by an argon plasma sustained by inductive 
coupling to an rF (radio frequency) field. An ICP torch consists of three concentric 
quartz tubes as shown in Figure 1.12. An induction coil, powered by an rf generator, 
circles the top of the tube assembly. Argon flowing between the two inner tubes is 
initially ionized by a spark form a Tesla coil. An annular plasma forms when the 
resulting ions and electrons interact with the oscillating magnetic field generated by the 
induction coils. The analyte is introduced as an aerosol through the central tube. A gas 
(usually argon) flows tangentially between the two outer tubes to contain the plasma and 
to cool the quartz cylinder walls. 
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Figure 1.12. Inductively Coupled Plasma Atomizer 

Samples are usually introduced into the plasma as a solution aerosol generated 
by a pneumatic nebulizer. The most commonly used pneumatic nebulizers are the 
concentric (or Meinhard) nebulizer and the crossflow nebulizer. In both types, the 
sample is drawn through a capillary into a low pressure region generated by the flow of 
a gas past the capillary tip. The aerosol produced in the nebulizer passes through a 
spray chamber, which removes or breaks up the larger droplets, and into the ICP torch. 

When the sample is introduced into the plasma it is atomized and elevated to an 
excited state as a result of collisions with the argon ions. The excited analyte atoms 
relax to their ground state by emitting photons. The emitted radiation is dispersed by 
either a polychromator, for simultaneous multi-element analysis, or by a scanning 
monochromator, for sequential multi-element analysis. 

In the polychromator, the emitted light is focused onto a concave grating where it 
is dispersed into its component wavelengths. Selected spectral lines are isolated by a 
series of exit slits. Each line is focused onto a detection device, a photomultiplier tube. 
The polychromator permits the simultaneous detection of up to 60 spectral lines 
(corresponding to 60 elements). Simultaneous multi-element analysis is especially 
useful in applications involving routine analyses of large numbers of samples having 
similar elemental composition. In the scanning monochromator, a lens focuses the 
emitted light onto a concave mirror, which collimates the light onto a planar grating 
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mounted on a computer-controlled stepper motor. The grating disperses the light, and a 
final mirror focuses the light onto an exit slit before the detector. The wavelength is 
selected by rotating the grating. Only one wavelength can be detected at a given time. 
For multi-element analysis, the grating must be driven in a sequential manner to a 
position corresponding to the appropriate analytical wavelength for each element to be 
measured. The scanning monochromator is more flexible than the polychromator 
because it offers a greater choice of analytical wavelengths. This is an advantage in the 
development of new methods, or in laboratories where samples vary widely in elemental 
composition. 

1.2.4 X-ray Fluorescence Spectroscopy 

X-ray fluorescence (XRF) spectroscopy is used for qualitative and quantitative 
multi-element analysis. XRF can qualitatively identify all elements of atomic number 
greater than eleven, and quantitatively measure all elements of atomic number greater 
than fourteen present within a sample at the parts-per-million (ppm) or greater level. 
This technique is extremely useful because it is readily applicable to most solid and 
liquid samples with minimal sample preparation. 

XRF is based on the emission of characteristic x-ray lines by atoms following 
excitation by high energy photons from an x-ray source. The most commonly used x- 
ray source is the x-ray tube which consists of a tungsten filament cathode and a target 
anode, both sealed within a highly evacuated tube. Electrons are thermally emitted from 
the cathode when the filament is heated. These electrons are then accelerated across a 
high potential gradient to the target anode. Target materials include copper, 
molybedenum, iron, chromium, nickel, silver, and tungsten. X-rays are produced by 
bombardment of the target material. X-rays emitted by the source are then directed onto 
the sample. Upon absorption of the x-radiation, atoms within the sample become 
electronically excited. The excited atoms relax to their ground states by fluorescent 
emission of x-rays of characteristic energy. Sample fluorescence is collimated, and then 
dispersed by a crystal mounted on a goniometer. The goniometer permits control of the 
angle q between the incident collimated radiation and the crystal face. The collimated x- 
rays are diffracted from lattice planes within the crystal. The value at which a given 
wavelength 1 is diffracted is given by Braggs law, nX = 2d sin 0, where d is the lattice 
spacing of the ciystal. X-rays diffracted by the analyzing crystal are detected by gas- 
filled detectors or by scintillation counters. 

Compared to ICP/AES and AA, XRF is more rapid and more readily applicable 
to various samples . However, ICP/AES and AA are more sensitive than XRF. 

1.2.5 Infrared Spectroscopy 

Infrared spectroscopy (IR) is used to identify organic and inorganic materials, to 
elucidate molecular structures, and to quantitatively determine nontrace components of 
mixtures. This technique can be applied to solid, liquid, or gaseous samples. 
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The infrared region of the electromagnetic spectrum includes radiation with 
wavenumbers between 12,800 and 10 cm'l (wavelengths between 0.78 and 1000 »m). 
Most analytical applications make use of the mid-infrared region which encompasses 
wavenumbers from 4400 to 200 cm'l. 

Absorption of infrared radiation by a molecule results in vibration of the 
molecule's component atoms relative to each other. Two types of vibrations occur: 
stretching and bending. Stretching involves a change in the distance between two 
atoms, with movement along the bond axis. Bending can occur in any molecule having 
three or more atoms and involves a change in the angle(s) between bonds. 

A molecule will absorb infrared radiation only when the molecule undergoes a 
net change in dipole moment as a result of its vibrational motion. For example, the 
hydrogen chloride molecule possesses a significant dipole moment due to non- 
symmetric charge distribution between the hydrogen and chlorine atoms. Stretching of 
the hydrogen-chlorine bond causes a change in dipole moment, and absorption of 
infrared radiation can occur provided the frequency of the radiation matches the 
vibrational frequency. It is not necessary for a molecule to possess a permanent dipole 
moment to absorb infrared radiation. For example, due to its symmetry, the carbon 
dioxide molecule has no permanent dipole moment. However, two of its three 
vibrations produce a change in dipole moment. These vibrations absorb infrared 
radiation. Vibration of homonuclear diatomic molecules, such as N 2 and O 2 , does not 
cause a net change in dipole moment. Vibrations of such molecules are not 
accompanied by absorption of infrared radiation. 

Infrared spectra of diatomic and triatomic molecules are simple. Polyatomic 
molecules containing a number of different types of atoms exhibit complex spectra. The 
number of possible vibrations within a molecule containing N atoms is 3N-6 (or 3N-5 
for a linear molecule). Each of these vibrations is called a normal mode of vibration, and 
its frequency is referred to as a fundamental frequency. The number of peaks observed 
in the infrared spectrum does not necessarily equal the number of normal modes. Some 
normal modes do not give rise to an infrared absorption peak because no net dipole 
moment occurs, or the vibrational frequency is beyond the range of the instrument. If 
two or more vibrations occur at the same or nearly the same frequency, only one peak 
may appear. Additional peaks may be observed in the spectrum due to interaction or 
coupling of two normal modes. Overtone lines appear at approximately two and three 
times the fundamental frequency. The occurrence of multiple vibrations in a molecule 
gives rise to a complex spectrum that is uniquely characteristic of the molecule. 

The frequency at which an organic functional group (such as C-H, C=0, C=C) 
vibrates is determined by the masses of the atoms involved and the force constant of the 
bond between them. This frequency, called a group frequency, is often unchanged, or 
only slightly changed, by other atoms attached to the functional group. Over the years, 
group frequencies have been determined empirically for a large number of functional 
groups and are commonly summarized in correlation charts, as shown in Figure 1.28 
and 1.29. Group frequencies are used to establish the presence or absence of a 
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functional group in a molecule. Fingerprint frequencies are due to vibrations of the 
molecule as a whole and are characteristic of the specific molecule. 

Several types of infrared instruments are commercially available. Most widely 
used are the dispersive instruments, which use a grating for wavelength selection, and 
the popular nondispersive Fourier Transform Infrared (FTIR) spectrometer, which uses 
an interferometer. Whether dispersive or nondispersive, each infrared instrument 
contains three essential elements: a radiation source, an optical system for wavelength 
selection, and a detector. The components of an FTIR instrument are shown in Figure 
1.13. 



Figure 1.13. Schematic of an Infrared Spectrometer 

Mid-infrared sources include the incandescent wire, the Nemst glower, and the 
Globar. The incandescent wire source consists of a tightly wound coil of nichrome wire 
electrically heated to about 1100° C. The Globar is a silicon carbide rod electrically 
heated to about 1300°C. The Nemst glower is composed of zirconium, yttrium, and 
erbium oxides electrically heated to about 1500° C. Globar and Nemst glower sources 
are hotter and, therefore, more intense than the incandescent wire source. However, the 
incandescent wire is more rugged and requires less maintenance. 

Two types of detectors are used in the mid-infrared region: thermal and photon . 
In a thermal detector the infrared radiation is absorbed by a detector element. The 
resultant rise in temperature produces a measurable change in a physical property of the 
detector. A thermal detector may consist of several thermocouples, each fabricated from 
two dissimilar metals. When a thermocouple absorbs infrared energy, a measurable 
potential difference develops at the junction of the two metals. Another type of thermal 
detector is the pyroelectric detector. In this device, the sensing element is a thin crystal 
of a pyroelectric material, such as deuterated triglycine sulfate (DTGS), between two 
electrodes. Heating of the crystal by absorption of infrared radiation alters the 
polarization of charge within the crystal, and a measurable change in capacitance is 
produced. In a photon detector, infrared photons striking a semiconductor surface excite 
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electrons on the surface from a nonconducting energy state, or valence band, into a 
conducting state. The photovoltaic detector is a type of photon detector in which a small 
voltage is produced in response to infrared radiation exposure. Two examples of the 
photovoltaic detector are the lead tin telluride detector and the mercury cadmium 
telluride detector. Photon detectors respond more rapidly than the thermal detectors, but 
the two types have similar sensitivities. 

Infrared spectra can be obtained for gaseous, liquid, and solid samples. Most 
widely used as cell windows are the alkali halides, especially sodium chloride and 
potassium bromide. Unfortunately, these materials are fogged by exposure to moisture. 
Silver chloride can be used for aqueous solutions or most samples. However, silver 
chloride darkens with continued exposure to light. 

Solid and liquid samples can be analyzed in solution. No single solvent is 
transparent over the entire infrared region. To obtain the entire spectrum of a sample, 
two or more solvents must be used. Carbon tetrachloride and carbon disulfide are useful 
for many organic compounds. Carbon tetrachloride is transparent at wavenumbers 
above 1333cm' 1 , whereas carbon disulfide is transparent below 1333 cm'l Acetone, 
acetonitrile, chloroform, and methylene chloride are used for polar materials that are 
insoluble in carbon tetrachloride and carbon disulfide. Because alkali halide cell 
windows are attacked by moisture, solvents must be dried before use. Water absorbs 
strongly in the infrared and is seldom used as a solvent for infrared analysis. Infrared 
cells for solution samples are commonly constructed with sodium chloride windows 
separated by Teflon spacers or gaskets. Because solvents absorb infrared radiation, path 
lengths are short (0.005 to 5 mm). Both variable path length and fixed path length cells 
are available. 

When a suitable solvent is not available, spectra of liquids can be obtained from 
capillary films. A drop of liquid is placed between two windows, which are squeezed 
together to give a thin (0.001 to 0.05 mm) film. This technique does not give a 
reproducible path length and is limited to qualitative investigations. 

Spectra of solids not soluble in a suitable solvent are often obtained from a 
dispersion of the solid sample in a solid or liquid matrix called a mull. To prevent signal 
loss due to scattering or refraction of the infrared radiation, the sample must be 
thoroughly ground to a particle size smaller than the analytical wavelength (less than 
2*m) and dispersed in a mulling agent whose refractive index is close to that of the 
sample. For best results, a very small amount (less than 10-20 mg) of sample is ground 
in an agate mortar until the powder forms a glossy cake on the sides of the mortar. The 
resulting powder is then mixed with the mulling agent. 

Mineral oil (Nujol) is frequently used as a liquid mulling agent. If hydrocarbon 
bands interfere, a chlorofluorocarbon grease (Ruorolube) or hexachlorobutadiene may 
be used. To obtain a complete spectrum free of interfering bands, the halogenated 
agents are used for the 4000 to 1340 cm' 1 region, and mineral oil is used below 1340 
cm -1 . One or two drops of mulling oil are added to the powder in the mortar, and the 
grinding action continued until the solid is uniformly dispersed in the oil. The resulting 
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paste is then analyzed as a thin film. 

Potassium bromide and other alkali halides are used as solid mulling agents. In 
the KBr mull or pellet technique, the finely ground sample (up to 1 mg) is intimately 
mixed with approximately 100 mg powdered KBr. The mixture is then pressed in a die, 
under a pressure of 10,000-15,000 psi, to form a transparent or translucent disc. Better 
results are often achieved when the die is evacuated to eliminate occluded air. Adequate 
KBr pellets can be obtained without evacuation in a stainless steel Mini-Press. This 
device consists of two highly polished bolts within a cylindrical nut. The KBr mixture is 
placed between the bolts, and the bolts are tightened against each other to produce a 
pellet. 


The KBr pellet method has several advantages: (1) KBr has no interfering 
bands, (2) quantitative analysis is readily possible with an internal standard, and (3) the 
formed pellets can often be stored. However, the method also has disadvantages: (1) 
Some materials decompose or change form under the pressure and heat encountered 
during pellet formation, (2) inorganic salts may exchange ions with the KBr, and (3) 
KBr may absorb water which absorbs near 3450 and 1640 cm~l. 

Polyurethane and polyisocyanurate foams, cork, and epoxy resins are examples 
of materials that cannot be reduced to a powder by grinding in a mortar. A mechanical 
grinder or ball mill can be used for such materials. Grinding efficiency may be 
enhanced by first freezing the sample with liquid nitrogen. 

Polymers can be analyzed as unsupported films. A polymer solution is poured 
onto a glass or metal casting plate. The solvent is evaporated in a vacuum oven or under 
an infrared lamp, and the sample film is stripped from the plate. Samples run as 
unsupported films frequently exhibit interference fringes which interfere with the 
spectrum. In these cases, the attenuated total reflectance (ATR) technique, described 
below, is extremely useful. 

Attenuated total reflectance ATR (also called internal reflectance) greatly 
enhances the sampling capability of infrared spectroscopy. With this technique, infrared 
spectra can be obtained from solid materials that cannot be readily dissolved, dispersed 
in mulls, or cast as films. ATR is widely used to sample polymers, rubbers, cured resins, 
fibers, textiles, and papers. Figure 1.14 shows an internal reflection element (IRE), 
composed of a transparent material of high refractive index, surrounded by a sample. 
When radiation passes from a more dense medium into a less dense medium, part of the 
radiation is transmitted through the interface and part is reflected. The fraction of the 
radiation that is reflected increases as the angle of incidence at the interface between the 
two media becomes larger. Beyond a certain critical angle, all of the radiation is 
reflected. During the reflection process, the radiation penetrates a very small distance 
into the sample, and a portion of the radiation is absorbed by the sample. The loss in 
intensity due to sample absorption is sensed by the detector. 
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Figure 1.14. Attenuated Total Reflectance Apparatus 
1.2.6 Raman Spectroscopy 

Raman spectroscopy is used primarily for the identification of functional groups 
and the determination of molecular structure in organic and inorganic compounds. 
Similar to infrared spectroscopy, Raman spectroscopy is based on vibrational changes 
within molecules. However, vibrations that are observed in the Raman spectrum may 
not be seen in the infrared spectrum, and vice versa. Therefore, complimentary 
information can be obtained from using both techniques. 

Raman spectroscopy is based on a light scattering phenomenon. When light 
passes through a transparent material, a very small fraction of the light is scattered due 
to collisions with molecules present within the material. Two types of collisions occur: 
elastic and inelastic. In an elastic collision, the scattered light is of the same energy as 
the incident light. However, in an inelastic collision, energy is exchanged between the 
molecule and the incident radiation so that the energy of the scattered light is different 
from that of the incident light. This is known as the Raman effect. When this type of 
collision occurs, the energy of the scattered light is h(no • n n ) where h is Plank’s 
constant, no is the frequency of the incident radiation, and n n is a frequency in the 
infrared region of the electromagnetic spectrum. The Raman effect is inherently very 
weak. The intensity of the Raman scattered light is, at most, only 0.001% the intensity 
of the source radiation. 

Inelastic collisions give rise to two types of scattered radiation: Stokes radiation 
and anti-Stokes radiation. Stokes radiation, observed at a lower energy that the incident 
radiation, occurs when a molecule initially in the ground vibrational level is raised to an 
excited vibrational level as a result of interaction with the incident light. The energy of 
the Stokes radiation is h(nQ-n n ) where hn n is the difference in energy between the 
ground vibrational level and the excited vibrational level. Anti-Stokes radiation occurs 
when molecules already in an excited vibrational level decay to the ground vibrational 
level during the interaction with the incident light . Anti-Stokes radiation is observed at a 
higher energy than the incident radiation. The energy of the anti-Stokes scattered 
radiation is h(no+n n ). Because most molecules are initially in the ground vibrational 
level, anti-Stokes radiation is much less intense than Stokes radiation. Except for very 
specialized applications, the Stokes lines are predominantly used in analytical Raman 
spectroscopy. 
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In the previous section, it was stated that for a vibration to be observed in the 
infrared spectrum, the vibration must cause a change in the molecule's dipole moment. 
The dipole moment is the product of the distance between two centers of charge in a 
molecule and the magnitude of the charge difference. Polar function groups such as 
hydroxyl (O-H) and carbonyl (C=0) give strong peaks in the infrared. For a vibration to 
be observed in the Raman spectrum, the polarizability of the molecule must change as a 
result of the vibration. Polarizability is the measure of the ease with which the electron 
distribution within the molecule is distorted in an electromagnetic field. Nonpolar 
functional groups with symmetrical charge distribution, such as C=C, N=N, S-S, and 
C=C, give strong Raman peaks. Some vibrational modes can be observed in both the 
Raman and infrared spectra. Other vibrations can be observed only by infrared 
spectroscopy, or only by Raman spectroscopy. Therefore, although each technique 
provides considerable information about the functional groups present in a molecule, 
more information is obtained by using both techniques to characterize a material. 

Raman spectroscopy is especially useful for detecting the C=C group in olefins. 
The C=C stretching mode is usually very weak in the infrared spectrum. When the C=C 
group is symmetrically substituted, as in ethylene (H 2 C=CH 2 ), no band appears in the 
infrared spectrum because there is no change in dipole moment during the vibration. On 
the other hand, this group gives an intense peak in the Raman. Similarly, Raman 
surpasses infrared in detecting disubstituted acetylenes (R-C=C-R) and disulfides 
(R-S-S-R). Even in cases where the two R groups are not identical, but only similar in 
size, the intensity of the Raman peak will be stronger than that of the infrared peak. 
Raman spectroscopy can also be used to detect and analyze homonuclear diatomic 
molecules, such as H 2 or N 2 , which are not infrared active. Another area in which 
Raman spectroscopy is useful is the study of materials in aqueous solution. Water is a 
weak Raman scatterer and therefor does not cause interference. In contrast, it is difficult 
to obtain infrared spectra from aqueous solutions because water absorbs strongly in the 
infrared. 

A typical Raman spectrometer is shown schematically in Figure 1.15. Because 
of the inherent weakness of the Raman effect, lasers are used almost exclusively as 
sources in Raman spectroscopy. The most commonly used lasers are continuous-wave 
gas lasers. The helium/neon laser produces an intense line at 632.8 nm. The argon laser 
produces useful lines at 488.0 and 514.5 nm. The laser beam is focused on the sample 
by a series of lenses and mirrors. Because both the incident and scattered radiations are 
in the visible region of the electromagnetic spectrum, glass can be used for cell 
windows, lenses, and other optical components. The scattered light is collected at 90* 
relative to the laser beam and is focused on the entrance slit of a double monochromator. 
The scattered light is dispersed into its component frequencies by gratings within the 
monochromator, and the dispersed light passes through the monochromator exit slit into 
the detector. The most commonly used detector for Raman spectroscopy is the 
photomultiplier tube (PMT). PMT's give low background signals and high sensitivity in 
the visible region of the spectrum. Because PMT response is a function of wavelength, 
the type of PMT used depends on which laser line is used for excitation. 
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Raman spectroscopy can be performed on solids, liquids, and gases. Liquids are 
usually contained in a glass capillary tube. Solids are sampled in several ways: as pure 
powders in a glass capillary tube, as pressed pellets (either pure or mixed with an inert 
solid), and as single crystals. A visible light microscope can be coupled to the optical 
system and can be used to obtain spectra from microscopic particles. A major limitation 
to Raman spectroscopy is interference by fluorescence originating either in the analyte 
itself, or in other species present in the sample. Because fluorescence is inherently more 
intense than Raman scattering, it is difficult if not impossible to detect Raman scattering 
in the presence of fluorescence. If the fluorescence is due to other species present, it 
may be possible to remove the fluorescent component by extraction, distillation, or 
another purification procedure. If the fluorescence is due to the analyte itself, selection 
of a different laser line for excitation may be helpful. Another problem is the 
decomposition of the sample in the laser beam. This occurs frequently when the sample 
is highly absorbing (colored). This problem can sometimes be overcome by rotating the 
sample or by reducing the power of the laser source. 

1.2.7 Nuclear Magnetic Resonance Spectroscopy 

Nuclear magnetic resonance (NMR) spectroscopy is used primarily to elucidate 
the structure of organic and inorganic compounds containing atoms which possess a 
magnetic moment, such as *H, 19p, 29si t anc j 31p 

When the nuclei of certain atoms spin, a magnetic moment is generated. If these 
nuclei are placed in a strong magnetic field, the magnetic moments will assume different 
orientations with respect to the external field. Each orientation corresponds to a discrete 
energy state. Absorption of radio-frequency radiation causes a nucleus to undergo a 
transition from one orientation to another. The required energy of the radio-frequency 
radiation depends, in part, on the the identity of the nucleus and chemical environment 
of the nucleus. Therefore, this technique is extremely useful in the determination of 
molecular structure. 

The components of a Fourier transform NMR spectrometer are shown in Figure 


32 




1.16. In this instrument, all of the nuclei are excited simultaneously by a pulsed, high 
intensity radio-frequency source. This generates a time-domain spectrum that is 
converted to frequency-domain by Fourier transformation. A computer controls the 
pulses and performs the Fourier transformation. 



Figure 1.16 Schematic of a Nuclear Magnetic Resonance Spectrometer 

Two effects are important in structural analysis: chemical shift and spin-spin 
splitting. Not all nuclei of the same type (e.g., *H or 13 C) absorb at the same frequency. 
The frequency at which absorption occurs is strongly affected by the chemical 
environment of the nucleus. This effect, called the chemical shift, is due to the 
circulation of electrons around the atom being probed. Electron circulation generates a 
small magnetic field which is usually in opposition to the external field. The field 
actually sensed by the nucleus depends on the electron density and, therefore, the 
chemical environment around the nucleus. NMR peaks are usually reported relative to 
the peak obtained for a reference compound that is measured at the same time. 
Neighboring nuclei which have nonzero magnetic moments generate magnetic fields. 
This effect causes splitting of the NMR peaks and is called spin-spin splitting. 
Structural analysis is based on an analysis of the chemical shifts and spin-spin splitting 
observed in the NMR spectrum. 

1.2.8 Mass Spectrometry 

Mass spectrometry provides qualitative and quantitative information about the 
atomic and molecular composition of organic and inorganic materials. This technique is 
used to elucidate the molecular structure of unknown compounds, and to confirm the 
presence of known compounds. Combined with a separation technique, such as gas or 
liquid chromatography, mass spectrometry can also be used to identify and 
quantitatively measure the components of a complex mixture. 

In mass spectrometry, molecules within a sample are converted into highly ener- 
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getic, gaseous ions. During the ionization process, the molecule may be fragmented to 
produce smaller ions. These ions are then separated on the basis of their mass-to-charge 
(m/z) ratios. A mass spectrum is a plot of the relative abundance of ions measured at 
each value of m/z. Every compound has a unique fragmentation pattern, and therefore a 
unique mass spectrum. 

The major components of a mass spectrometer are pictured in Figure 1.17. The 
sample introduction system must permit introduction of the sample into the ion source, 
which is maintained at 10‘5 to 10'6 Torr. Sample introduction systems include batch 
inlets, direct insertion probes, and gas and liquid chromatographic inlets. Gases and 
volatile liquids can be introduced through a batch inlet system in which the sample is 
first volatilized (by applying vacuum and, sometimes, heat) and then allowed to leak 
into the ionization chamber through a pinhole in a glass or metal diaphragm. Solids and 
nonvolatile liquids can be introduced with a direct insertion probe. The sample is placed 
on the probe which is then introduced directly into the ionization chamber. The 
chromatographic inlets are described in Section 1.1. 



Figure 1.17. Schematic of a Mass Spectrometer 

The ion source ionizes and fragments the molecules within the sample. 
Commonly used ionization sources include the electron impact (El), chemical ionization 
(Cl), and fast atom bombardment (FAB) sources. The electron impact source is the 
most commonly used. Electrons emitted from a heated filament are accelerated by 
applying a potential between the filament and an anode. When these high energy 
electrons collide with the sample molecules, positive ions are produced. Some of these 
ions have the same mass as the original analyte molecule, and are referred to as the 
molecular or parent ions. Because of the high energy of the electron impact source, 
many of the molecules are extensively fragmented. The fragmentation pattern is 
characteristic of the molecule and is useful for identification. However, in some cases 
fragmentation is so extensive that the molecular ion is not observed. In these cases, it is 
more difficult to determine the molecular weight and, therefore, the structure of the 
compound. In chemical ionization, the gaseous sample is ionized by collisions with ions 
generated by electron bombardment of a reagent gas, usually methane. Chemical 
ionization produces less fragmentation than electron impact ionization. In fast atom 
bombardment, the sample is ionized by collisions with high energy argon or xenon 
atoms. Fast atom bombardment produces relatively large amounts of the parent ion, and 
is especially useful for the characterization of high-molecular-weight materials. 
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After leaving the ion source, the ions are accelerated and focused onto the 
entrance of the mass analyzer. The mass analyzer separates the ions according to their 
mass-to-charge ratios. The magnetic sector and the quadrupole mass analyzers are the 
most commonly used. Magnetic sector mass analyzers utilize an electromagnet to bend 
the paths of the ions. Ions of different mass are alternately focused on the exit slit by 
varying the strength of the magnetic fields. In quadrupole mass analyzers, four 
electrically conducting rods are positioned parallel to and symmetrically around the ion 
path. Opposite pairs of rods are electrically connected, and one pair is held at a positive 
dc potential, while the other pair is at a negative dc potential. In addition, radio- 
frequency ac potentials, 180° out of phase with each other, are applied to the two pairs. 
The combined electromagnetic fields cause the ions to oscillate as they pass along the 
axis of the quadrupole. When the oscillations are sufficiently large, the ions collide with 
the rods and are neutralized. By varying the dc potential and/or the frequency of the ac 
field, ions of different mass-to-charge ratio are allowed to pass through the analyzer. 

Most mass spectrometers use electron multipliers for ion detection. Ions leaving 
the mass analyzer are accelerated and focused onto a metal plate, or dynode. The 
dynode emits several electrons as a result of each ion bombardment. These electron are 
accelerated toward a second dynode, where more electrons are emitted. As the electrons 
strike successive dynodes, the amplification process continues. Current amplification of 
10 5 to lO? is attained. 

1.3 Microanalytical Techniques 

The advent of high technology has brought about the need to interrogate smaller 
and smaller samples of materials to solve complex materials problems. From the 
concerns of particulate contamination in high performance aerospace engines to the 
performance of magnetic storage media there is a great demand for the capability to 
analyze very small quantities of materials. And with this demand for analysis comes the 
necessity to successfully collect and test the material of interest. When the sample 
comes from an almost endless source such as in the production of commodity materials, 
the problem is not too little sample, but in contamination or failure investigations the 
samples are usually limited and hard to recover. This section will describe sampling and 
detection methods for microscopic analysis. Sampling techniques will be discussed in 
detail since this topic is seldom covered in instrumental analysis text books. 

1.3.1. Microsampling Techniques 

Microsampling techniques are usually applied to materials with limited quantity 
and can either be solids or liquids. The solid samples can include particles, fibers, films, 
laminates, inclusions, suspensions or residues. Each of these materials requires a 
different approach to their collection. Another confounding aspect to the sampling is the 
accessibility of the material. If the material of interest is under a film as an inclusion or 
in a solid mixture as a discrete particle the analyst must be creative, and sure-handed to 
retrieve the materials. 
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Some of the tools that can be used for microsampling include tweezers with 
sharpened tips, Exacto knives and single-edged razor blades, tungsten probes that have 
been formed into fine points, glass slides, capillary brushes, Roller Knife and maybe the 
most important is a prep microscope. The prep microscope is used to see samples when 
they are being transferred or prepared for analysis. The prep microscope can have 
magnification from lOx to 50x and should have the capability for transmission and 
reflectance imaging. It should have a substantial fixed stage that will support some 
force in forming samples. Any variation from the Bausch and Lomb StereoZoom series 
or equivalent would be a good choice. Pamphlets 2 supplied with the microscope 
give some instruction as to the use of the microscope but experience from others or with 
time will better develop knowledge needed in specific applications. Publications such as 
"The Microscope"^ is also a good source of information for microscopic sample 
collection. 

Collecting samples in the field requires many specialized tools to effectively 
isolate and collect the sample. Along with the techniques and tools to be described later 
there are several items that are useful and should be included in a Microscopists tool kit. 
The most important of these items is a magnifier. This is almost a necessity since to 
sample it is necessary to see the object. The magnifier can be used with a flashlight or 
can be equipped with an appropriate illuminator as an integral part of the magnifier. 
Sample transfer from the point of sampling is also important since the effort expended in 
sampling can be negated if the sample is lost or contaminated with unrelated foreign 
matter. Foresight in selecting sample containers that have been previously inspected for 
microcontamination will benefit the analyst in the long run. Plastic petrie dishes and 
self-sealing bags are convenient but require prior inspection to assure that particles are 
not present as a result of static attraction or contamination during production. 

1.3.1. 1. Particulate Sampling 

Particulate sampling is probably one of the easier sampling problems but by no 
means simple. Particulates that are loosely adhered to a surface can be easily removed 
using tweezers with very sharp, straight points. This fact needs to be emphasized; don’t 
share the sampling tools that you accumulate with anyone that is not familiar with their 
use. It might seem trivial that anyone can use tweezers but to care for a finely sharpened 
pair and then have them used as a pry to remove a stuck septa or ferrule can be very 
disconcerting. Particles that are brittle, very flat or adhere strongly to the substrate may 
require the use of a sharp blade such as a single-edge razor, or an Exacto knife. 
Cleanliness is paramount when using any sampling device and should be of special 
concern with blades. A small foreign particle on the blade prior to collection of the 
sample can become a point of contention when the sample is analyzed because this 
foreign particle can be misconstrued as part or a major portion of the suspect. The use 
of lens paper or lint-free laboratory wipes to clean the sampling tools prior to use 
followed by subsequent inspection using a microscope to observe cleanliness is a 
minimum. 
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In dry atmospheres charges tend to accumulate on sampling tools and particles. 
With loose samples this can manifest itself in the "Calivaris County Frog Syndrome" 
where the precious sample that is setting so pristinely on the substrate leaps 
uncontrollably out of sight. This charge concentration phenomena can be eliminated by 
using a static discharge gun such as a Zerostat gun (Ernest F. Fullam, Latham, N.Y). 
Two conditions should be observed when using this static discharge gun. First, never 
use the gun in an explosive environment and second, the gun discharges both positive 
and negative charges with each stroke ; a positive discharge with the squeeze of the 
trigger and negative discharge with the release of the trigger. It is important to discharge 
both the sample and the sampling tools prior to sampling. This is a trial and error effort 
to select the proper discharge combination that hopefully doesn't end in errors too 
frequently. Of course since opposites attract, static charges can also be a benefit since 
they can aid sample collection by causing particles to cling to the sampling tool, a very 
fortuitous event. 

In addition to the tweezers and blade, a tungsten probe is very useful to push or 
hold a particle while using another tool to retrieve it, or if static charges are helpful, to 
pick up the sample. A capillary brush (whose construction will be detailed later) is also 
a useful tool to collect large numbers of particles. The application of the capillary brush 
is limited in particulate collection since separation of the particles from the brush fibers 
might be problematic. An alternative to the capillary brush is the vacuum tip, a device 
constructed to provide a vacuum to pick up small particles and a filter medium to 
capture the particle for analysis. This apparatus can be constructed of a capillary and 
filter holder attached to an appropriate vacuum source. The use of adhesive tapes are 
quite useful in collecting and transporting particulates, but removing the sample from 
the adhesive can be cumbersome. An alternative is to dissolve the adhesive from the 
backing using filtered toluene and evaporating it onto the KBr window prior to analysis. 
The particle is surrounded in a matrix that can hold it securely and can be removed 
digitally from the collected spectrum by spectral subtraction. 

Particles are found not only on the surface, but can also be imbedded in objects 
as inclusions. Included particles are complex sampling problems. If the material in 
which they are contained is soft and pliable it is possible to excise them from the 
material using blades, probes, dental picks and tweezers. The best approach is to 
observe the particle under a prep microscope. A large diameter tungsten probe would be 
beneficial since it could be used to penetrate the substrate and pry the sample to the 
surface where it can be picked up with tweezers or a blade. The adherence of the 
surrounding material is of some concern, but it is of benefit and can be used in the 
analysis of the particle as a reference in spectral subtraction (the application of spectral 
subtraction will be covered elsewhere). This sampling technique is of course destructive 
in nature but can probably be limited to small areas. Substrates of an intractable nature 
are more problematic since the substrate can prevent penetration because of their 
hardness and strength. If probes and blades are not successful it might be necessary to 
cut the substrate and the particle using a diamond or jeweler's saw. The cut will have to 
be made first to remove a larger piece from the substrate to include the particle and then 
through the particle to expose it to the surface. A microtome can also be useful here to 
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take a small, thin slice of the sample. The microtome is an expensive item, so unless 
this type of sampling is a routine occurrence the effort expended in sawing is more cost- 
effective. 

Particles suspended in liquids can be handled very easily by using filtration 
media of the proper construction and observe cleanliness in sample handling. A typical 
apparatus that would be used to separate suspended particles would include a gas tight 
syringe of 5 ml or less supplied with a Luer-Lok connector. A syringe filter holder and 
filter medium constructed of materials compatible with the suspension solvent should be 
used. Filter media such as Nuclepore filters and equipment (Costar Corp. Cambridge, 
MA.) is a good selection since the filter media has uniform pore size and is constructed 
of nonfiberous, solvent resistant materials. It is important to use a filter medium free of 
fibers to avoid confusion of sample with medium. To aid in locating transparent 
samples more easily, the filter media can be sputter-coated with gold/palladium from a 
sputter coating apparatus used in Scanning Electron Microscopy if available. An 
aliquot of the suspension is passed through the filter and then flushed with an aliquot of 
the same solvent as the suspension to rinse the filter of any dissolved material that might 
crystallize after the solvent has evaporated. The flushing solvent should have been 
passed through a filter medium of at least the same pore size as that used for the 
sampling to avoid contamination of the sample with foreign particles. Once rinsed the 
filter media containing the particles of interest can be examined under a prep microscope 
and the particles can be removed from the filter by any appropriate method previously 
mentioned. 

1.3. 1.2. Fibers 

The problems that arise in fiber collection are far fewer than those experienced in 
particulate sampling. Fibers are an easier material to handle since they are flexible and 
usually have one large dimension. They can exhibit charging as do particles, that can be 
neutralized with static discharge guns. Tweezers are the most appropriate tool to use in 
fiber collection, but things as diverse as adhesive tape or vacuum filters can also be 
used. An easy vacuum method utilizes the filter support and media used to filter 
solutions and a vacuum source. These techniques are appropriate for isolated fibers or 
those that are loose and non-woven. For woven fabric or bundled fibers the sampling 
methods become more complex. In fabric the sampling of individual fibers can be 
performed by using a tungsten probe to lift the fiber from the weave, tweezers to pull out 
an appropriate length and a pair of sharp, small scissors to cut the selected fiber. With 
this method any number of individual fibers can be removed from the fabric for 
analysis. Fiber bundles can be treated in the same manner as fabric except for a need to 
immobilize the bundle. Immobilizing the bundles can be done with adhesive tape as a 
support or by holding the bundle between two glass slides so that their ends extend past 
the edge of the slide. It should be noted that some glass slides are lubricated with 
formate salts and should be removed by scraping one slide edge over the surface of 
another. This will eliminate the possibility of sample contamination. The individual 
fibers can then be selected using tweezers. If the fiber is coated or multilayered and a 
cross-section is desired so that each layer can be observed, the fiber can be mounted into 
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epoxy contained in the end of a Pasteur pipet The mounted sample can then be sliced 
or microtomed before analysis. 

Fibers are usually extruded and for that reason have a uniform cylindrical shape. 
This shape acts as a lens when exposed to the spectrometer and distorts the beam 
causing aberrations in the baseline and changes in relative absorptivities of some peaks. 
To avoid this problem it is best to flatten the fiber in the area of interest This can be 
performed using either a rigid spatula or a Roller Blade (Spectra Tech, Stamford CT). 
The Roller Blade can be used with the smallest fibers to roll a portion of the sample flat 
This flattened area is then ideal for analysis. 

1.3.1 .3. Films and Laminates 

The sampling of films and laminates can be divided into two classes that are 
dependent on the choice of analytical method. For thin films, transmission 
measurements would be the most appropriate if the film is transparent and thin enough 
to allow light to pass through. In this case the film should be treated as particulate 
material. If the film is too thick for transmission or opaque or a laminate, a thin cut can 
be taken from the edge of the sample using the following procedure. The film is placed 
between two glass slides (which have been cleaned as mentioned earlier) in such a way 
as to leave a portion of the film overlapping the edge. The exposed film edge is then cut 
flush with the glass slides. The sample is then repositioned so that a thin section overlaps 
the edge of the slide and a very thin cut is made from the edge of the sample. The use of 
a prep microscope is of benefit here. The thin section that was formed in the second cut 
should be almost translucent and can be handled as a particle. If the sample is a 
laminate each of the layers can be analyzed individually. The laminates can also be 
sampled by cutting obliquely through the sample thickness to expose wide surfaces of 
each laminate. As an alternative to sectioning, thick films of an adequate size can be 
sampled using organic-free abrasive paper such as is supplied with the Si-Carb Sampling 
Kit (Spectra Tech, Stamford, CT). The abrasive paper is used to remove a portion of the 
surface of the film and then analyzed by reflectance techniques directly on the pad. This 
sampling is limited to surfaces and would be ineffective in laminated structures. 


1. 3.1.4. Micro-pyrolysis 

Solid samples are at times intractable even to crushing. This requires harsh 
treatment of the sample in order to transform it into a more easily manipulated 
configuration. This can be done using micro-pyrolysis. Micro-pyrolysis is a technique 
that is used to convert insoluble, hard samples to a soluble liquid or solid. The method 
involves placing the sample particle or film into a capillary about 3 mm from the end. 
The sample end of the capillary is heat-sealed and the particle is then heated with a 
microbumer until pyrolysate droplets are formed. The sealed end of the tube containing 
the ash, if any, is broken off and solvent is allowed to flush the sample onto an 
appropriate surface for analysis. As was previously noted cleanliness of the solvent is 
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critical to the analysis to avoid contamination. 

1.3. 1.5. Liquid Micro-Sampling 

Liquids are probably the easiest samples to collect since they are not usually 
subject to static charges and infrequently do they adhere to the substrate. However, the 
small quantities that usually occur when collected for micro-analysis are easily lost in 
the collecting tools normally used for liquid collection, such as pipets or syringes. 
Certainly microsyringes are appropriate for this task since they hold very small volumes. 
When samples are not confined to small volumes, but spread out over a large area or in 
small droplets, it is difficult to collect samples. A solution to this is the use of the 
capillary brush. The capillary brush is a glass capillary fitted with glass fibers welded to 
the end. The process of making these brushes is tedious but worthwhile. Small groups 
of silica fibers such as those used in gas chromatography are placed in the end of the 
capillary. A microtorch is then used to melt the capillary slightly, enough to hold the 
fibers but not so much as to close the capillary. The fibers are allowed to extend out of 
the capillary and are trimmed to form an even tip. This brush can be used to collect 
samples from enclosed areas or broad surfaces by immersing the fibers into the liquid 
and allowing it to "capillary" into the fibers and capillary tube. Once the sample is 
isolated in the capillary it can be transported to the analytical station and removed from 
the capillary by flushing the capillary with a solvent by "capillary action" onto an 
appropriate surface for analysis. 


1.3.2. Micro- Spectroscopy Using Fourier Transform Infrared and Raman 
Spectroscopies 

Molecular vibrational spectroscopy is an important means of identifying the 
structure of chemical compounds. Each bond in a molecule has a specific response to 
incident light energy and will affect the energy contained in the incident light in a 
specific way. All wavelengths of the incident energy are not absorbed, but rather only 
those energies that interact with bonds in the sample will be affected. This absorption 
occurs at specific frequencies and is related to the bonds in the structure with which it 
has interacted. The interactions are not random but abide by specific rules described by 
quantum mechanics and are, therefore, reproducible and predictable. When these ab- 
sorptions are detected and plotted as a function of incident energy frequency it repre- 
sents an absorption spectrum. Each peak of this absorption plot can be assigned a spe- 
cific structure with which the incident light has interacted. These interactions cause the 
atoms and bonds to distort by stretching or bending. 

1.3.2. 1. Infrared Micro-spectroscopy 

The application of microscopic sampling to infrared spectroscopy has had a long 
history dating back to the 1940s when Burch^ applied a microscopic accessory equipped 
with Cassagranian lenses to collect microspectra using a dispersive infrared instrument 
The first application^ of microspectroscopy to Fourier transform instruments was in 
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1983 when BioRad Instruments introduced the FT-IR microscope. 

Figure 1.18 is a schematic of the beam path of a typical infrared microscope. It 
is similar to a sophisticated light microscope in that it contains both a substage 
condenser and an objective lens. These lenses are, however, reflective rather than 
transmissive lenses to accommodate the infrared beam. Standard glass lenses would be 
opaque through most of the mid-infrared region and therefore are not used. Another 
difference between the light microscope and the infrared microscope is that the latter 
includes a movable beamsplitter that allows viewing the sample during positioning and 
focusing, and variable apertures, in some cases, both before and after the sample. The 
apertures are used to mask both the incident beam and the transmitted beam. This 
"redundant aperturing" is used to assure that an undistorted baseline is collected by 
removing stray light (light which has passed around the sample). The aperture after the 
sample in the transmitted beam can be used to spectrally isolate the unwanted portions 
of complex sample from contributing to the collected spectrum. The spectral isolation 
permits analysis of small inclusions in larger matrices that would contribute to the 
sample’s spectrum. The aperture before the sample limits the light scattered by the 
substage condenser from reaching the detector. The microscope accessory can be used 
to collect spectra in both the transmission (through the sample) and reflection (from the 
surface) mode. Reflectance is more troublesome in the interpretation of the resulting 
spectrum since they are usually distorted by refractive index effects. These refractive 
index effects are manifested in the reflectance spectra by "derivative-shaped" peaks. 
The position of the peak is also shifted by several wavenumbers from typical positions 
found in transmission spectra. 
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Figure 1.18. Schematic of the Beampath of an Infrared Microscope 
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The operation of the infrared microscope is straight-forward and is very similar 
to a light microscope. Unlike light microscopy, however, the infrared radiation 
wavelength imparts problems to the collection of spectra from samples whose particle 
size is less than about 50 *m. Particles in this range will be of the same size or smaller 
than the infrared wavelength. Shape and thickness of the sample are critical since 
irregularly shaped samples are capable of greater diffraction of the incident beam. The 
diffracted beam then forms nodes that divide the incident energy and reduce the total 
throughput of the sample. The use of apertures again play a significant role here but this 
time to the detriment of total signal since they effectively eliminate some of the 
diffracted nodes. Some improvement in the spectrum is gained because the apertures 
eliminate the diffracted nodes from other areas of the sample, which is the spectral 
isolation function mentioned earlier. 

Most samples to be analyzed by infrared microspectroscopy are introduced into 
the accessory by placing the sample on an infrared-transparent material such as a KBr or 
NaCl window. The first step of the collection process is focusing the microscope on the 
sample. Focusing is important since it minimizes spectral distortions caused by irregular 
samples. It is good practice to focus deep into the sample if the sample has a rough 
surface. That is to say, rather than use the top surface of the sample use a focal plane 
that appears to be inside the sample. Flat samples are easily analyzed since they have a 
uniform surface and focusing on the top surface is adequate for data collection. 

The second step in the collection process is to select the area of interest in the 
sample by using the apertures to outline this area. This aperturing performs two practi- 
cal functions. The first is to define the portion of the sample to be analyzed and the 
second is to block radiation not passing through the sample or the area of interest. The 
radiation not passing through the sample will cause a shift or distortion of the baseline 
leading to spectral irregularities such as incorrect peak ratios within the sample and 
sloping baseline. The sloping baseline manifests itself in the display of the spectrum 
and in spectral subtraction (to be covered later). 

Once the sample is in place the third step is data collection. The FT-IR 
instrument is a single-beam spectrometer and it is necessary to collect a background 
spectrum to which the sample spectrum is ratioed. The sample in the microscope has 
been masked by an aperture and it is necessary to collect the background using the same 
aperture size as the sample. Since it is impossible to know the sample aperture size 
before-hand it is common to collect the background spectrum after the sample spectrum 
but before ratioing the sample to a background. The background spectrum is collected 
by moving the image field from the sample to an area of the supporting window that has 
no sample present. 

The collection of a good spectrum using microspectroscopy is best served by 
analyzing flat samples, but, in practice few samples are flat. Particles that are too thick 
to produce a spectrum of less than 2 AUFS (absorbance units full scale) should be 
crushed to reduce the thickness. Spectral absorbance above 2 AUFS will distort the 
relative ratios of the absorbing peaks, especially the peaks of greatest intensity. 
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Crushing can be done by placing the particle between two clean glass microscope slides 
and squeezing it until the sample is translucent. The crushed sample can then be trans- 
ferred to a KBr window for analysis. Two KBr windows can also be used to crush softer 
samples. When crushing between KBr windows, additional KBr powder should placed 
between the windows to provide a location for collecting a background. If no KBr were 
present an interference fringe would be formed by the gap between the two windows 
when the background was collected. The interference in the background would 
unnecessarily distort the sample spectrum during ratioing. In the case of elastomeric 
materials the difficulty in collecting a spectrum is not in flattening a thick sample, but 
maintaining the thickness after crushing. Because of their resilience the elastomers tend 
to return to their original shape. The solution to this problem is to compress the sample 
between two KBr windows and while holding them securely place fast-curing glue 
(Crazy Glue, for example) on the edge of the windows to hold them in place. Allow the 
glue to cure and then collect the spectrum. Samples that are difficult to secure on the 
slide can be immobilized on the KBr window by embedding them in rubber cement 
diluted 1:1 with cyclohexane. The spectrum can then be collected and the embedding 
media contribution to the spectrum can be removed by spectral subtraction. 

Single fibers are a difficult sample for classical transmission spectroscopy be- 
cause of their dimensions. In conventional spectroscopy the single fiber was difficult to 
mount reproducibly into the infrared beam and to observe without spectral distortion. 
Microspectroscopy has opened the field to capabilities that were non-existent in the past 
The infrared microspectroscopic accessory enables the analyst to collect high quality 
spectrum from single fibers as small as 5 »m in diameter. The fiber sample can be 
synthetic or natural of any dimension. The only restriction on the sample is that it 
should be capable of compression to form a flat area for analysis. This flat area can be 
formed using a hard probe to roll the fiber or a diamond anvil cell to crush the area of 
interest. In either case the fiber is flattened so that the it does not act as a cylindrical 
lens in the infrared beam. It should be noted here that modifying the fiber in such a 
manner will possibly alter any orientation information that can be obtained from the 
spectrum. 

Samples that are in the form of fine powders can pose a problem in conventional 
spectroscopy since the formation of a pellet will so dilute the sample that a complete 
spectrum is impossible to collect. The preparation of the sample as a salt dispersion or 
micropellet can be of benefit only when the sample can be observed. A salt dispersion 
can be made on a glass slide by mixing the sample with salt and then crushing the par- 
ticles to reduce their size. The micropellet provides an easier means of handling the 
sample and can be made using a microgrinder (such as a Wiggle-Bug) and dies that 
produce 3 mm pellets when compressed. The microscope provides the method of con- 
densing and focusing the beam to collect the spectrum. 

The infrared microscope is not only useful for solid materials, but is also useful 
for limited amounts of liquids. The deposition of the liquid sample onto a window by an 
appropriate technique such as a microbrush or microsyringe can form a spot that is 
dependent on sample size but can be focused on be the accessory. With careful concen- 
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tration on the window a small spot can be formed whose spectrum can then be collected 
as with any other transparent sample. Microspectroscopy is applicable as a tandem 
technique for gas chromatography since the spectrum of the liquid can be collected and 
then have the sample transferred to the GC for further investigation. It can also be ap- 
plied to HPLC eluents where evaporation of the mobile phase and resultant analysis can 
yield spectral information about the solute. 

1. 3.2.2. Raman Microspectroscopy 

The first application of microspectroscopy was reported in 1979.4 The advent of 
laser applied to Raman spectroscopy has enabled standard metallurgical microscopes to 
be coupled with the Raman spectrometer for the collection of spectra from small 
samples. Optical glass does not have a Raman interaction and is therefore readily 
adapted to Raman spectroscopy with little modification to the microscope optics. The 
Raman microscope has the same limitations as does Raman spectroscopy in that only 
samples with little color and low fluorescence can be analyzed. The spatial resolution of 
the microscope is excellent with sample sizes only limited by the focal size of the 
incident laser onto the sample. In practical terms this resolution is 5-10 »m. 

The operation of the Raman microscope is relatively straightforward for sample 
analysis. The sample is placed onto the sample stage with the microscope in the re- 
flectance mode, the area of interest is located, and focus is obtained for the sample and 
the laser spot on the sample. For samples that are thermally stable the spectrum can be 
collected using high laser energy input. The greater the incident energy the greater the 
emitted signal. Thermally labile samples are significantly more complex in the analyti- 
cal approach used. Because the high incident created in the microscope field by the 
laser generates significant amounts of heat, thermally sensitive samples must be handled 
with more conservative methods. 

The thermal deterioration is significant since repeated occurrences lead to the 
deterioration of the microscope objective through deposition of the evolved materials 
onto the objective. The first and most obvious means of reducing thermal destruction is 
to reduce the incident energy. This approach will reduce the thermal load on the sample 
but with proportional reductions in Raman effect emitted light. The compromise here is 
to increase the data collection time with either increased sampling rate or increased 
number of coadded scans. Both of these options increase data collection times but 
provide less aggressive environments. An alternative to incident energy reduction is to 
provide a means for heat to dissipate from the sample during analysis. The easiest and 
most straight forward methods are to embed the sample in water glass, a concentrated 
solution of sodium silicate that is transparent to the Raman effect and forms a hard, 
transparent mount for the sample. The water glass also protects the sample and objec- 
tives from thermal decomposition by providing a mode for the excess heat energy to be 
dissipated. Samples can also be diluted in a pellet made of materials that have no 
Raman effect such as alkali halides to minimize heating. Defocusing the incident beam 
is also an effective technique for mediation of localized heating but signal reduction may 
result in poor quality spectra. 
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Liquid samples of limited volumes can also be analyzed by Raman 
microspectroscopy. The application does not require complex sampling method but 
merely placing die sample in a glass capillary. Care should be taken in focusing the 
laser beam on the sample. It is necessary to first focus the visual image of the capillary 
surface and then while observing the laser image raise the stage to move the beam into 
the capillary. A distinct change in the laser image will occur when the beam has passed 
into the capillary contents. Data collection is identical to any other sample. 

Fiber samples can also be evaluated as with infrared microspectroscopy except 
that the fiber need not be crushed to reduce thickness but can be collected without 
modification. Samples can be immobilized on a glass slide using adhesive tape. In 
addition, surface particles on the fibers can be evaluated using appropriate pinhole 
masks to spatially separate the area from the matrix. Characteristics of the area of 
interest can also be accentuated in the final spectrum by spectral subtraction. 

1.3.3. Scanning Electron Microscopy 

Scanning electron microscopy (SEM) is used to produce high resolution and 
depth of field images of sample and to provide chemical analysis of micron sized areas 
on or near the surface. Products obtained from interactions between bombarding 
electrons and the atomic species in the specimen are ejected from an excitation volume 
in the sample and monitored during a typical SEM scan. Any partially conducting 
material or non-conducting materials coated with a thin layer of gold or gold alloy, 
which has a low vapor pressure can be analyzed with the SEM. 

In SEM analysis, the specimen is bombarded with high energy electrons (5 to 30 
kilovolts) which produces several different types of signals to provide information about 
the specimen. Both elastic and inelastic collisions result from the interactions of the 
electrons with the atoms in the sample. The elastic collisions between the bombarding 
electrons and the atomic nuclei of the specimen produce backscattered electrons (BSE), 
as shown in Figure 1.19, which provides both topographical and compositional 
information about the sample. Imaging with the BSE signals provides a means of 
distinguishing zones of different atomic number within a specimen due to the greater 
probability of elastic collisions with higher atomic number. The penetration depth of the 
bombarding electrons depends upon the accelerating voltage and the atomic number of 
the sample. Typically the BSE images provide useful information down to about 0.5 
microns. 

Inelastic collisions between the bombarding electrons and the atomic electrons of 
the specimen produce secondary electron imaging signals, x-rays and Auger electrons 
from interactions with inner shell electrons and light in the form of cathodo- 
luminescence from lower energy processes. The secondary electrons are collected to 
produce the micrograph images which are normally associated with SEM. The depth of 
penetration for these signals are usually 2 to 10 nm for metals and 5 to 50 nm for poor 
conductors. The x-ray information typically arises from a volume around 1 micron in 
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depth, while cathodo-luminescence is observed almost totally from the surface. Auger 
electrons arise from interactions in the first 0.5 to 2 nanometer layer of the specimen. 

The great depth of field of the SEM ( up to 500 times that of an optical microscope) 
allows it to produce completely in-focus images of rough surfaces at high magnification. 
However, SEM is generally inferior to the optical microscope for routine examination of 
samples prepared using standard metallographic techniques at low magnification (300 X 
- 400X). Some of the current SEM instruments allow for samples as large as 15 to 20 
cm to be analyzed. An excellent reference for using SEM in a variety of analytical 
applications is given in the monograph by Gabriel. 



Figure 1.19 Schematic showing the principal mechanisms occurring in SEM. 

1.3.4. Optical Microscopy 

Conventional optical microscopy is a relative inexpensive analytical tool for the 
examination of particulates and microstructures. Typically resolutions of 500 
nanometers are possible, which allows for observation of surface faceting and 
dislocations. Catalogs of particulates, fibers, and morphological species of many 
materials are available for use to identify a particular material. Many applications using 
polarized light and the new confocal microscopes enables optical microscopy to 
continue to provide analytical information which is not available from the other 
methods. 

A major disadvantage of current uses of optical microscopy is that experienced 
microscopists can identify many species; however, computerized databases with 
efficient pattern recognition capabilities are not readily available for the industrial 
laboratory. Such data can be included in chemical fingerprinting databases and will 
make the technique even more useful. 
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1.3.5. Transmission Electron Microscopy 

Transmission electron microscopy is extremely useful for characterizing 
materials with high resolution imaging from small (submicron) regions of the sample. 
The information obtained form TEM analysis includes microstructural analysis of 
metallic, ceramic and polymeric materials as well as chemical and crystallographic 
information. Certain crystallographic details such as crystal orientation and matrix- 
precipitate orientation relationships are also possible. The advent of scanning 
transmission electron microscopes (STEM) has expanded the role of this analytical tool 
from primarily new chemical structures and failure analysis to activities such as 
fingerprinting. 

TEM and STEM normally use an electron beam (varying from 60 to 300 
kilovolts for TEM and 2- 50 kilovolts for STEM) to pass through a thin specimen and 
the resulting intensity distibution is imaged by a photographic imaging system or 
electron energy loss spectrometer. The interactions between the electron beam and the 
sample provides several types of information. Chemical and density effects result in 
intensity variations much like an optical microscope, except that the resolutotion of the 
electron microscope is much is near the wavelength of the electron beam. In addition 
diffraction effects from the crystal lattices present in the specimen can also be used to 
identify or characterize the material. 

The capability to attach an energy dispersive x-ray spectrometer on the side of 
the instrument column allows collection of x-rays generated by interactions between the 
focused electron beam and the specimen; enabling characterization or identification of 
the chemical or elemental composition of submicron volumes possible. The limitations 
of current energy dispersive spectrometer systems is for the atomic number must be 6 or 
higher. 


1.3.6. X-ray Diffraction 

X-ray diffraction is extremely useful for identying and characterizing a large 
variety of materials. In general, x-ray analysis is restricted to crytalline materials, 
although some information may be obtained from amorphous and powdered samples. 
Most x-ray diffraction techniques are rapid and nondestructive and are capable of 
identifying the phases, grain size, texture, and crystal imperfection. Samples are 
acceptable in many different forms, depending upon the availability of the material and 
the type of analysis to be performed. 

The X-ray diffraction technique is based upon diffraction effects of x-rays 
traveling through a sample. This phenomen is characterized by the Bragg equation: 

n?L 

sin 0 = -j* 
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where X is the wavelength of the x-rays and d is the characteristic interatomic distance 
of the material. A number of spots or fringes are observed at various angles represented 
by the sin 9. These spots or fringes result from constructive interference of the diffracted 
beams and are specific for the crystal lattices of each material. Hence the ability to 
identfy or characterize a material is quite good. 

Current x-ray diffraction instrumentation can be portable. Also, there is a 
substantial amount of software available today supporting the interpretation of x-ray 
diffraction data as well as databases associated with compilations from different 
laboratories. 

1.3.7. Low Energy Electron Diffraction 

Low energy electron diffraction is a technique for characterizing surfaces and 
overlayers absorbed on surfaces. It can be considered the surface analog of x-ray 
diffraction. Low energy electrons (30 - 200 electron volts) have limited penetration into 
materials and consequently provides high sensitivity for periodic atomic structure at the 
surface. The periodic atomic structures at the surface diffracts a monoenergetic electron 
beam like a diffraction grating. The diffraction pattern and its intensity distribution can 
provide information about the atomic positions and the crystallographic features of the 
surface. 

The diffraction pattern obtained during a LEED analysis consists of a small 
number of spots displayed as a CRT image as whose symmetry of arrangement is that of 
the surface grid of atoms. Any surface contamination on the surface of the material will 
significantly affect the diffraction pattern observed and provides a useful capability for 
discriminating against undesired surface features. Catalogued LEED patterns are 
available for many single-crystalline materials. 


48 



The diffraction pattern obtained during a LEED analysis consists of a small number of 
spots displayed as a CRT image as whose symmetry of arrangement is that of the surface grid 
of atoms. Any surface contamination on the surface of the material will significantly affect the 
diffraction pattern observed and provides a useful capability for discriminating against 
undesired surface features. Catalogued LEED patterns are available for many single-crystalline 
materials. 

1.4. Surface Science 

Surface Science has evolved over the past 25 years to become widely accepted for the 
evaluation and characterization of solid vacuum compatible materials used in the aerospace 
industry. Although there are many techniques available, the most practical are: 

ESCA or XPS - Electron Spectroscopy for Chemical Analysis 

or X-ray Photoelectron Spectroscopy 
AES - Auger Electron Spectroscopy 
SAM - Scanning Auger Microscopy 
SIMS - Secondary Ion Mass Spectroscopy 
ISS - Ion Scattering Spectroscopy 

Each of these techniques exhibits certain advantages when used alone, but when two or more 
techniques are combined they exhibit a synergism that constitutes a better understanding of the 
process or material being studied. 

Except for OM or CM, the primary excitation sources used by these techniques include 
photons, electrons or ions. The primary particles interact with the surface and emit secondary 
particles (i.e., photoelectrons. Auger electrons, secondary ions or scattered ions). The emitted 
secondary particles are analyzed using either an energy or mass analyzer. 

The information gained from surface analysis includes the elemental and chemical 
composition of typically the top 2.0 - 5.0 nanometers (top few atom layers) of the sample 
surface. The surface sensitivity of surface analysis is due to the inelastic mean free path of low 
energy electrons in solids. The energy of the emitted electrons is typically between 30 and 2000 
eV. Such low energy particles allow only those originating from the very near surface of the 
material to escape and be detected. Particles originating deeper within the solid will interact 
with neighboring atoms loose their energy and be absorbed by the solid or contribute to the 
background signal. Generally, the elemental sensitivity of surface analysis is of the order of 
0.5 % of an atom layer which further makes surface science very sensitive to surface 
contamination. 

Surface science can be used on a wide variety of failure and special materials 
investigations or characterizations. Some of the materials investigations where surface science 
has proven most valuable include: 

• Chemical Staining of Processed 2219 T87 Aluminum 

• Epoxy Primer Pinhole Investigations 
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• Waste Minimization 

• Aluminum/Lithium Weld Development 

• Weld Defect Investigation 

• Epoxy Primer to Metal Adhesive Debonds 

The following is a brief discussion of each of the surface science techniques mentioned 
above, their advantages, disadvantages, types of samples analyzed, method for technique 
selection and examples. 

1.4.1 Electron Spectroscopy for Chemical Analysis and X-ray Photoelectron Spectroscopy 

ESCA is a technique which uses an X-ray source of known energy for the primary 
source of sample excitation (Figure 1.20). 



Figure 1.20. Diagram of the ESCA Process. 

Typically, a dual magnesium (1253.6 eV) and aluminum (1486.6 eV) source is used. 
However, it has been found that the use of a high energy zirconium (2042.4 eV) anode coupled 
with the magnesium offers the greatest utility when analyzing aluminum materials used on the 
ET program. The sample, once mounted and introduced into the ultra-high vacuum chamber, is 
placed at the focal point of the electron energy analyzer. The sample is than irradiated with 250 
or 300 watts each of either Mg, Zr or Mg/Zr simultaneously. These photons interact with the 
elements on the sample surface and eject core level photoelectrons which escape the sample 
surface and are passed through the electron energy analyzer and detected. Each element has its 
own set of characteristic photoelectrons of known kinetic energy. The kinetic energy of the 
ejected photoelectron is governed by the equation: 


e K.E. = hv - e B.E. - 0 
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where Ek.E. = Kinetic Energy of the Ejected Photoelectron 
hi) = Incident X-ray Photon Energy 
Eg.E. = Binding Energy of the Ejected Photoelectron 
0 = Work Function of the Electron Energy Analyzer 

By knowing the analyzer work function, incident photon energy and measuring the 
kinetic energy of the ejected photoelectron one is able to determine the binding energy of the 
photoelectron to within typically 0.1 eV. ESCA not only provides the elemental identification 
of the elements present but also offers the ability to determine the chemical state of the atoms 
on a solid surface via small shifts in binding energy. 

The ESCA technique is such that it can be performed easily on any solid vacuum 
compatible material. Consequently, another key advantage of ESCA is its ability to analyze 
both electrically conductive (i.e. metals, metal oxides, etc.) as well as insulating materials (i.e. 
powders, polymers, catalysts, etc.). ESCA is not limited by the electrical properties of the 
sample being analyzed. 

As an example, an ESCA survey scan of a chemically processed 2219 T87 aluminum 
surface was analyzed (Figure 1.21). ESCA identified the presence of carbon, oxygen, fluorine 
and aluminum on its surface. The survey scan was acquired with the x-ray source powered at 
250 watts each Mg/Zr and the spectrum normalized to zirconium at 2042.4 eV. Consequently, 
the kinetic energy at the 0 eV binding energy is the zirconium fermi edge. 



Figure 1.21. Typical ESCA survey scan of a clean chemically processed 2219 T87 aluminum 
surface normalized to Zr at 2042.4 eV. 
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Another key advantage of ESCA is the ease of quantitation with the use of developed 
peak area sensitivity factors. If the sample is assumed homogeneous within the analysis 
volume, the intensity of the photoelectron peak is given by the equation: 

I = nfo0y?iAT 

where I = The intensity of the photoelectron peak 

n = Number of atoms of a element/cnw 
f = Incident x-ray flux in photons/cm^-sec 
a = Photoionization cross-section 
G = Instrumental angular efficiency factor 
y= Efficiency of the photoelectron process 
h= Mean free path of a photoelectron in a solid 
A = Area of sample analyzed 
T = Detection efficiency 

Thus, solving for the number density of atoms yields: 

n = I/(foGyXAT) 

Defining the denominator as the Atomic Sensitivity factor (S) reduces the equation to: 


n = I/S . 


If we consider the ratio of two elements in a homogeneous matrix we have the expression: 

ni/n 2 = (Ii/Si)/(I 2 /S 2 ) . 

The generalized expression for determining the atomic fraction of any element in a sample is 
than represented by the equation: 

C x = (n x /Sni) = (I x /S x )/S(Ii/Si). 

This expression, although far from absolute, offers the ability to perform semiquantitative 
ESCA analysis to within 10 - 20 % on homogeneous samples. When more accurate quantitative 
analysis is desired, standards of known concentrations which accurately represent the samples 
of interest must be analyzed. In this way, escape depths, chemical effects, matrix effects, etc. 
are similar and thus cancel each other. 

Quite often the main disadvantage of ESCA is its large analysis area. Conventional 
ESCA instruments have a very large analysis area. For example, some instruments have a fixed 
analysis area of about 12 ram^. Although most samples are sufficiently large and unaffected by 
this constraint, some samples are far too small for accurate ESCA measurements. One such 
example is pinholes in epoxy primer. The majority of these anomalies are typically about 1 mm 
in diameter (< 1 mm^) which is 1/12 the analysis area of the instrument. Newer instruments 
have been developed which render such analyses routine. These small area ESCA instruments, 
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have variable analysis areas and are able to analyze surface features as small as 75 »m. The 
minimum effective analysis area of a small area ESCA unit is 0.004 mm 2 which is far below 
that required for the analysis of pinholes. Small area ESCA was performed on the pinhole 
anomaly and found to be extremely effective. 

1.4.2 Auger Electron Spectroscopy & Scanning Auger Microscopy 

Auger Electron Spectroscopy differs from ESCA in that AES uses an electron beam for 
sample excitation and occurs as a result of a relaxation process (Figure 1.22). 



Figure 1.22. Diagram of the Auger Process. Note that AES is a Three Electron Process. 

Whereas ESCA uses an X-ray source to eject photoelectrons for analysis, AES uses a finely 
focused electron beam to eject core level electrons thus placing the surface atoms in an excited 
state. The excited surface atoms relax to the ground state via radiation or the Auger process. 
The energy of the primary electron is typically between 3 and 10 kV. Most AES analyses are 
performed using 5 kV. The Auger process involves three electrons. The first electron (Ej) is 
ejected by the primary electron (Ep). A second electron (E2) from a higher energy level fills 
the core hole created by the ejection of Ej. The excess energy between Ej and E2 is translated 
to a third electron (E3). The third electron is the Auger electron. The Auger process places the 
atom in a doubly ionized state. The Auger process is independent of excitation source and each 
element beginning with lithium has its own set of characteristic Auger transitions. AES is also 
sensitive to 0.5 % of an atom layer or 0.1 atomic %. The Auger electrons are also passed 
through an electron energy analyzer and detected. The result is an elemental identification of 
the elements present on the sample surface. 

As an example, a typical AES survey scan of the same aluminum surface discussed 
above identified the presence of carbon, oxygen, sulfur, chlorine, nitrogen, fluorine, copper and 
aluminum (Figure 1.23). The AES spectrum was acquired using a 5 kV and 1 mA electron 
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beam. The spectra is displayed as the derivative of the signal (d(N/E)/dE) as a function of 
energy. 



Figure 1.23. Typical AES survey scan of a clean chemically processed 2219 T87 aluminum 
surface acquired using a 5 kV, 1mA electron beam. 

Quantitation of AES spectra is another advantage of performing AES analysis. The 
method employed is somewhat similar to that of ESCA, however, the derivative spectra or the 
peak-to-peak is used as a measure of intensity. This method for quantitation is less accurate than 
the peak area measurement used by ESCA. 

One of the key advantages of AES over ESCA is its small analysis area. AES has a high 
spatial resolution due to the use of a finely focused electron beam for excitation. Older 
instrument models have a minimum AES analysis area of 8 X 10*7 mm^ (~ 1 mm diameter), 
the newer commercially available instruments have analysis areas of the order of 5 X 10*10 
mm^. Regardless of the instrument, AES is orders of magnitude better suited than ESCA to the 
elemental identification of small surface features. 

Still another advantage of AES is its ability to scan the electron beam over an area of the 
sample and obtain maps showing the elemental distribution of differing elements. This is the 
SAM feature of the surface analysis instrument. Instruments can operate effectively at 
magnifications of 1000X. As an example, SAM was found very useful in the investigation of a 
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weld defect Here samples were machined from scrapped flight hardware where defects were 
known from x-ray inspection to exist. The samples were then mounted, introduced into the 
ultrahigh vacuum analysis chamber (5 X lO' 1 ^ torr) and fractured in-situ. In-situ fracture 
exposed the weld defects under UHV which were than analyzed by both AES point analysis and 
SAM. This approach allowed the identification of the contamination causing this anomalous 
condition. 

Another example of the SAM capability was the analysis of a contaminated aluminum 
surface. AES point analysis found the composition of the surface to vary from point to point 
Additionally, a 250X SEM image contained bright spots randomly dispersed on its surface. The 
elements present were aluminum, oxygen and carbon. Thus, 250X SAM elemental maps were 
acquired for these elements (Figure 1.24). 



Figure 1.24. 250X SEM (bottom-left) and SAM elemental maps for 
aluminum (top-left), oxygen (top-right) and carbon (bottom-right). 
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The SAM maps revealed aluminum and oxygen surface voids which were found to be 
rich in carbon. The carbon was found to be dispersed islands on the surface of the sample. 
These analyses were performed at a magnification of 250X. 

Although AES is a very powerful analytical tool, it does have disadvantages. A major 
disadvantage of AES is the types of samples which can be analyzed. Unlike ESCA, AES is not 
capable of analyzing electrically insulating materials. Due to the use of an electron gun for 
excitation, samples must be electrically conductive and have a conductive path to ground. 
Thus, AES and SAM analysis is limited to semiconductive and conductive materials. 
Furthermore, AES does not offer the same level of chemical state information gained by ESCA. 
However, AES is typically a quicker technique to perform. AES was not employed in the 
epoxy primer pinhole investigation due to the fact both the contamination and primer materials 
are non-conductive. 

1.4.3 Secondary Ion Mass Spectroscopy 

SIMS is the most sensitive of the surface science techniques to surface contamination. 
SIMS is typically several orders of magnitude more sensitive than either AES or ESCA. The 
lower detection limit for many elements by SIMS is of the order of 1 ppm. Furthermore, SIMS 
is the only surface science technique capable of detecting hydrogen and is sensitive to the entire 
periodic table of elements. 

Secondary Ion Mass Spectroscopy is performed using an Argon inert gas ion gun for 
excitation and a quadrupole mass analyzer (QMA) for detection. The SIMS process involves 
the Argon being ionized and accelerated onto the sample surface with typically 0.5 to 4 kV. 
The incident Argon ions interact with the surface species of the sample through exchange of 
energy and momentum. The surface species are ejected from typically the top two atom layers 
of the sample (~ 5 A) as positive ions, negative ions and molecular fragments (Figure 1.25). 
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Figure 1.25 Illustration of the SIMS Process. 

The emitted secondary ions and ionic fragments are then focused into an energy filter, 
passed through a mass analyzer and detected. The SIMS spectra are displayed as the intensity 
(log or linear scale) vs atomic mass unit (AMU). Similar to ESCA, SIMS analysis can be 
performed on both electrically conductive and non-conductive materials. One major 


56 



disadvantage of SIMS over ESCA and AES is that it is a destructive technique. Once SIMS is 
performed on a sample surface whether static or dynamic the original as received surface has 
been destroyed. Consequently, SIMS is usually the only or last technique to be performed on a 
sample surface. 


As an example of the utility of SIMS, analysis was performed on a chemically processed 
2219 T87 aluminum surface. Both positive and negative SIMS spectra were acquired (Figure 
1.26). 



Figure 1.26. Positive (top) and negative (bottom) SIMS spectra of a 
chemically processed 2219 T87 aluminum surface. 
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The spatial resolution of SIMS is governed by the spot size or raster area of the incident 
ion gun. A minimum spot size of 800 *m translates into an analysis area of 0.5 mm^. 
Typically, however, a 2 mm X 2 mm raster area is analyzed (4 mm^ ). 

Another major disadvantage of SIMS is the difficulty in the quantitation of SIMS data. 
This is due to the fact that the elemental sensitivities for SIMS cover a very wide range (i.e. 
10 4 ). SIMS is an excellent complimentary technique to both ESCA and AES. Furthermore, 
good sample to sample or good to bad comparisons can be made of similar sample materials. 

Compared to the AES and ESCA spectra discussed earlier, SIMS is far more sensitive to 
certain surface species. SIMS identified the presence of sodium, magnesium, potassium, 
calcium, silicon, titanium, vanadium, chromium, copper, iron, zirchonium, and hydroxide on a 
similar surface analyzed by both ESCA and AES. 

1.4.4 Ion Scattering Spectroscopy 

Ion Scattering Spectroscopy is similar to SIMS in that it uses an inert gas incident ion 
beam for sample excitation. ISS uses a monoenergetic low energy (0.1 to 3 kV) helium, neon or 
argon ion beam as its excitation source. A fraction of the incident ions interact with the sample 
surface and are scattered from the surface atoms by binary elastic collisions (similar to the 
scattering of billiard balls of different masses). The ISS process (see Figure 1 .27 on the next 
page) is governed by the equation: 


El = (l+p? COs6 + ' si ° 2 e ) L2 J 2 


where Ej 

E0 

Mi 

M 2 

4 

e 


= Energy of the Scattered Ions 
= Energy of the Incident Ions 
= Mass of the Incident Ions 
= Mass of the Surface Atoms 
= M 2 /Mj 

= Scattering Angle 


An ISS spectrum is a measure of the energy distribution of scattered incident ions 
leaving the sample. ISS peaks occur at specific values of Ej. Furthermore, the scattered ion 
signal is proportional to the number of surface atoms for a given element, probability that the 
incident ions will remain ionized and differential scattering cross-section. The differential 
scattering cross-section increases with atomic number such that ISS is most sensitive to higher 
mass elements. The probability that an incident ion will remain ionized following scattering is 
10^ and decreases dramatically with subsurface atoms which gives ISS its unique ability to 
analyze the very top atom layer of a material. Thus, the primary advantage of ISS is its unique 
surface sensitivity. 
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As an example, ISS was performed in the characterization of staining of chemically 
processed 2219 T87 aluminum. Both AES and ESCA identified the presence of copper in the 
area of the stain. These techniques, however, were unable to identify the elements at the very 
top of the sample surface. ISS clearly identified the presence of an elevated amount of copper in 
the stained area compared to an unstained area (Figure 1.28). 
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Figure 1.28. He3 ISS survey spectra of a stained (top) 
and unstained (bottom) chemically processed aluminum surface. 

Some of the disadvantages of ISS are similar to SIMS. Quantitation is difficult and the 
surface of the sample may be changed due to minor sputtering of the surface species during 
analysis. Still another common problem encountered with ISS is its limitation in mass range 
and resolution with any given incident ion species. Helium is able to detect the lighter elements 
beginning with lithium and offers the greatest range of sensitivity. Neon is sensitive to 
elements greater than magnesium while argon is sensitive to elements greater than titanium. 

1.4.5. Technique Selection Considerations for these Techniques 

The above discussions explained in some detail the four surface science techniques. It is 
important to keep in mind that no one technique is considered "THE" most appropriate 
technique for a given sample. AES and ESCA are by far the two most popular. However, the 
best approach to be used in a given analysis is to perform multiple techniques on a given 
sample. Thus, when possible, both ESCA and AES should be performed followed by either ISS 
and/or SIMS. To assist with the appropriate selection of techniques to be used, refer to the 
attached Table 1.3 on the next page. 
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Table 1.3 Selecting the Appropriate Surface Analysis Technique 



AES/SAM 

ESCAorXPS 

ISS 

SIMS 

CHARACTERISTIC 

Auger Electron Spectroscopy 
Scanning Auger Spectroscopy 

Electron Spectroscopy for 
Chemical analysis or X-Ray 
Photoelectron Spectroscopy 

Ion Scattering Spectroscopy 

Secondary Ion Mass 
Spectroscopy 

Excitation Source 

Focussed Electron Beam 

Mg and/or Zr X-ray source 

He, Ne, or Ar ion beams 

Primary - Ar ion beam 
Secondary - He & Ne ion beams 

Spatial Resolution 

1-5 Micronmeters 

4 millimeters 

2x2 millimeters 

800 micrometers \ 

Depth Resolution 

1 - 5 nanometers 

1 - 5 nanometers 

First atomic layer 

Beginning with first atomic 
layer 

Sensitivity 

Good sensitivity throughout the 
atomic number range. Spectra 
become more complicated with 
atomic number. 

Similar to AES/SAM 

Higher sensitivity for larger 
atomic number elements 

Excellent sensitivity for low 
atomic number elements. Also 
sensitive to molecular fragments 
and isotopes. 

Detection Limit 

0.1 Atomic Percent 

0.1 Atomic Percent 

1 atomic Percent 

~ 1 ppm for many elements. 
Most sensitive technique. 

Materials 

Excellent for metals and thin 
metal oxide films. Materials 
must be somewhat conductive. 

Excellent for all solid vacuum 
compatible materials. Least 
destructive of the four 
techniques 

Best suited for metals and metal 
oxides. Somewhat surface 
destructive. 

Best suited for metals and metal 
oxides. Applicable to both 
conductive and non-conductive 
materials. Technique is surface 
destructive by its nature. 

Elemental Identification 

Lithium through Uranium 

Lithium through Uranium 

Lithium through Uranium 

Hydrogen through Uranium 

Chemical Bonding 

Very limited, e.g. A1 vs A1?0^ 

Excellent for differentiating 
between chemical states of 
materials 

Not applicable 

Not applicable 

Ion Etching 

Excellent due to small anlysis 
area @ .5 - 50 nm/Min. 

Only fair due large analysis area 
@ .5 - 5 nm/min. 

Poor due to low energy ions 
used. 

Good due to small analysis area 
required by the very nature of 
the technique. 0.5 - 25 nm/min. 

Depth Profiling 

Excellent for elemental 
composition analysis as a 
function of depth 

Poor for elemental composition 
as function of depth. Excellent 
for chemical composition as a 
function of depth. 

Poor 

Good for elemental composition 
as function of depth. Excellent 
for molecular fragment 
composition as a function of 
depth. 
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1.5 . Thermal-Mechanical Techniques 

A number of thermal and mechanical methods of analysis of materials can be very 
useful for chemical fingerprinting purposes. Several instruments for frequently used for 
determining materials' thermal and mechanical properties using temperature programmed 
analysis include hight temperature differential thermal analysis, (DTA), thermogravimetric 
analysis (TGA), thermomechanical analysis (TMA) and differential scanning calorimetry 
(DSC). Included in this category of instrumentation is also the very useful technique described 
as rheometric dynamic spectrometry (RDS). These methods have primarily been used to 
characterize polymers and their precursors; however, numerous laboratories use them to 
characterize metals and ceramics as well. A general discussion of these methods follows. 

High temperature differential thermal analysis consists of instrumentation that measures 
the temperature of a milligram sized specimen as heat is applied at a controlled rate. The upper 
limit of these systems is generally around 1500 °C. The heat capacity of the material is easily 
obtained, as well as latent heat effects such as phase changes or transformations and 
autoignition.. The normal display of the temperature changes during the programmed heating 
provides indication of contaminants or out-of-specification material. The critical glass 
temperature, which is a measure of the purity of the material, can also be a useful characteristic 
for the material. 

Thermogravimetric analysis is performed by placing a milligram-sized sample onto a 
temperature programmed microbalance. The upper limit of most systems is around 1000 °C. 
Fingerprinting can ber performed by determining whether a material's thermal weight loss is 
characteristic of the normal material. This technique can also be used to evaluate wiring covers, 
crystallinity, stabilizer addition and the thermal stability of many types of polymeric and 
ceramic materials. 

Thermal mechanical analysis is performed by placing a sample between a vitreous silica 
platform and a movable silica rod and subjecting the sample to a variable load. In dilatory, the 
sample is wetted with a liquid and the swelling behavior followed by the measurement of the of 
the movement of the silica rod. TMA is performed at various controlled temperatures, where 
characteristic properties such as thermal relaxation phenomena and tensile compliance in 
polymers, fibers, and powders can be obtained. 

Differential scanning calorimetry is similar to differential thermal analysis, except that 
the temperature of the sample is monitored for a controlled heat input. The upper limit of most 
systems is around 500 ° C when using aluminum pans to contain the samples. Fusion between 
the aluminum pans normally occurs when measurements are made above that temperature. 
Materials properties which ar obtained in this technique include glass transition temperature, 
crystalline transition temperatures, latent heat tranformations, and specific heat measurements. 

Rheometrics dynamics spectrometry determines the viscous and elastic properties of a 
material by imposing a cyclical load on the material at programmed temperatues and 
measureing the resultant torques and body forces. Dynamic viscosity, dynamic modulus and 
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loss angle can be calculated form this data. Sample temperatures can be programmed from -150 
to 400 0 C. 

1 .7. Strategy for Instrument Hierarchy 

Analytical instrumentation and instrumental techniques frequently used for chemical 
analysis and materials characterization were described in the previous section. The theoretical 
basis, application areas, sample requirements, and limitations of each technique were discussed. 
The tables and flow charts in this section were developed to serve as guides in the selection and 
application of instrumental techniques to solve a given analytical problem. 

Before selecting specific techniques to be used, the analytical problem should be defined 
by asking the following types of questions: 

What is the nature of the sample? 

- Is it a solid, liquid, or gas? 

- Is it organic or inorganic? 

- Is it a mixture, or a pure substance? 

- What is already known about its composition? 

- Is it permissible to destroy the sample during testing? 

- How much sample is available for analysis 

- What is the material’s history, and its future? 

What kind of information about the sample is desired? 

- Is an elemental or molecular analysis desired? 

- Is complete identification of all species present required? 

- Or, is it sufficient to identify the major and minor components? 

- Which specific species or components are to be analyzed? 

- What is the required precision and accuracy of the data? 

A few of the questions will be immediately answered by visual observation of the 
sample. Others can be answered readily through discussions or correspondence with the person 
requesting the analysis. More information may be provided by persons experienced in the 
testing or processing of the material, or by the material's supplier. The remaining questions will 
require a physical or chemical measurement to be performed. 
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An analytical problem can be attacked in three ways: 

• Consult with leading experts in the field 

• Search the literature for a solution 

• Experiment based on theoretical predictions 

The first two points should not be overlooked. Do not risk reinventing the wheel by 
failing to consult others who may have had experience with the same type of problem. 
These people may be within your company, may be even in your own lab. The 
companies from whom you have purchased analytical instrumentation may have 
applications chemists available to help you. Some instrument manufacturers publish 
newsletters and product bulletins that may provide valuable information. A literature 
search may turn up the perfect solution to your problem. On the other hand, if the 
problem involves a proprietary material, literature searches may be of little help. When 
these sources of information have been exhausted, it is time to experiment. 

One source of information about most materials that are used in industrial and 
laboratory processes that may contain chemical information is the Materials Safety Data 
Sheet (or MSDS). These information packages are readily available to the user to ensure 
that personnel safety is always maintained and necessarily do have to identify any toxic 
precursors which exist in the product. 

Unfortunately, no single procedure exists for the selection of the most effective 
instrumental approach to a given analytical problem. Factors to consider include the 
capabilities and limitations of the available instrumentation, and the quality of the 
analytical results versus the cost in time and materials required to obtain them. Table 
1.4 surveys the more common instrumental techniques and summarizes each technique’s 
application areas, limitations and sample requirements. This table can be used to 
identify the techniques that could be applied to the problem at hand. The remaining 
tables and charts contain additional information that will help the analyst to "narrow the 
field" and arrive at specific instrumental methods that can be applied to the problem. 
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TABLE 1.4. COMMON INSTRUMENTAL ANALYSIS TECHNIQUES 


Method 

Applications 

Advantages 

Method Limitations 

Sample Limitations 


Atomic Absorption 
Spectrometry (AAS) 

Quantitative and/or trace 
analysis of a single element 
for each measurement 

Fast, reliable, high sensitivity 
for some 70 elements, 
relatively inexpensive 

Not applicable to most non- 
mctallic materials or 
simultaneous multi-element 
analysis, 

Small linear response range 

Rrequires time-consuming 
dissolution of sample or 
graphite furnace for atomizing 
solids. 

mg tog 

Inductively Coupled Plasma- 
Atomic Emission 
Spectrometry (ICP-AES) 

Quantitative multi-elemental 
analysis Determination of 
trace, minor and major 
elements 

Simultaneous detrmination of 
up to 60 elements. Good for 
refractory materials. 

Large dynamic range 

Limited sensitivity for 
nonmetals. 

Expensive 

Requires time-consuming 
dissolution of sample or 
graphite buffer for atomizing 
solids. 

mg tog 

X-ray Fluorescence 
Spectrometry (XRF) 

Quantitative analysis of all 
elements of atomic no. > 14. 
Qualitative determination of 
all elements with atomic no. 
>11. 

Minimal sample preparation. 
Inexpensive 

Detection limits not as good 
as AA or ICP/AES 

Solid or nonvolatile liquid 

mg 

Infrared Spectrometry (IR or 
FTIR) 

Identification and structural 
determination of materials, 
including surface adsorbants 

Applicable to most materials. 
Extensive libraries of 
reference spectra available. 

Composition limited to 
molecular species identified. 
Medium sensitivity. 

Trace and minor components 
can be masked by major 
components. 

Material must contain bonds 
which undergo dipole moment 
change during vibration 

Rg 

Raman Spectroscopy 

Identification and structural 
determination of materials, 
including surface adsorbants 

Identification of non-polar 
functional groups. Minimal 
sample preparation. 

Low sensitivity. 

Fluorescence can interfere 
with Raman signals. 

Limited libraries of reference 
spectra. Expensive. 

Material must contain bonds 
which undergo polarizability 
change during vibration 

1 to 100 
mg 

Nuclear Magnetic Resonance 
Spectroscopy (NMR) 

Structural determination and 
identification of both organic 
and inorganic materials 

Determination of molecular 
configuartion and 
conformation 

Applicable only to samples 
containing magnetic moment. 
Low sensitivity. 

Expensive 

Sample must be liquid or 
soluble solid 

1 - 100 mg 


Elemental inorganic 
identification of elements 

Low level detection 

Elements > Na 

Sample must be liquid or 
soluble solid 

- 3 cm 
diameter 
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TABLE 1.4 COMMON INSTRUMENTAL ANALYSIS TECHNIQUES (continued) 


Method 

Applications 

Advantages 

Method Limitations 

Sample Limitations 

Sample 

Size 

Mass Spectrometry 

Structural determination and 
identification of organic and 
some inorganic materials. 

Widely applicable to most 
materials. Extensive on-line 
reference libraries available. 

Slow 

Expensive 

Sample must be volatized. 

0.01 g 

Gas Chromatography (GC) 

Separation of multicomponent 
mixtures of volatile materials. 

Selectivity ranges from 
general to specific 

Not applicable to nonvolatile 
and thermally unstable 
materials 

Material must be volatile and 
thermally stable 

Pg to mg 

Liquid Chromatography 
(LC or HPLC) 

Separation of multicomponent 
mixtures of liquids and 
soluble solids 

Widely applicable to 
nonvolatile organics. 
Applicable to 

thermallyunstable materials. 
Separated materials can be 
identified by other methods. 

No sensitive universal 
detector. Subsequent analysis 
by IR or MS necessary to 
identify components. 
Moderately expensive. 

Must be soluble in one of 
many suitable solvents. 

pg to mg 

Ion Chromotography (IC) 

Separation of complex 
mixtures of ionic species for 
quantitative analysis. 
Elemental analysis of organics 
after decomposition. 

Applicable to a wide range of 
organic and inorganic anions 
and to many cations. 

Analysis of trace species in 
presence of high concentration 
species is difficult. Method 
development is time 
consuming. Moderately 
expensive. 

Must ionize in solution. 
Nonaqueous applications 
limited. Decomposition of 
organics time consuming. 

1 to 5 mg 

Size Exclusion 
Chromatography (SEC) 
Gel Permeation 
Chromatography (GPC) 

Separation of complex 
mixtures based on molecular 
size. Determination of 
polymer molecular weight 
distribution. 

Applicable to polymers. 
Determines molecular weight 
distribution. 

Calibration time consuming. 
Moderately expensive. 

Must be soluble in limited 
number of suitable solvents. 
GPC perfomed in water. 

Pg to mg. 

Combined Gas or Liquid 
Chrmatography with Mass 
Spectrometry 
(GC/MS or LC/MS) 

Separation, identification and 
quantitative anlysis of 
complex mixtures. 

Combines separation 
capability of GC or LC with 
identification and sensitivity 
of MS 

Slow 

Method develoment is time 

consuming 

Expensive 

Same as GC, LC and MS 

20-200 

ng. 
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TABLE 1.5 . OTHER USEFUL INSTRUMENTAL ANALYSIS TECHNIQUES 


Method 

Applications 

Advantages 

Method Limitations 

Sample Limitations 

HB 

Differential Thermal Analysis 
(DTA) 

Quantitative characterization 
of materials & contaminants 
Dilatometry, autoignition 

Separates materials by 
differences in thermal 
properties in inert or oxidizing 
atmosphere 

Very few; can handle a large 
variety of solids and liquids, 
such as foams, films, 
powders, or fibers. 

Small sections that fit in 
sample crucibles 
Sample containment can be 
expensive 

mg 

Thermal Gravimetric Analysis 
(TGA) 

Quantitative characterization 
of materials, wire coatings 

Measures resistance to 
thermal degradation. Inert or 
oxidizing atmosphere. 

Very few; can handle a large 
variety of solids and liquids, 
such as foams, films, 
powders, or fibers. 

Small sections; determined by 
size of sample pan and 
autobalance range. 

mg 

Differential Scanning 
Calorimetry (DSC) 

Quantitative characterization 
of polymers & inorganics by 
heat absorption, T„ 

Differences in thermal 
properties easily measured. 
Inert or oxidizing atmosphere. 

Very few; can handle a large 
variety of solids and liquids, 
such as foams, films, 
powders, or fibers. 

Small sections determined by 
sample pan size. 

mg 

Thermal Mechanical Analysis 
(TMA) 

Semi-quantitative 
characterization of polymers 
& organics 

Distinguishes materials by 
differences in thermo- 
mechanical properties. 

Very few; can handle a large 
variety of solids and liquids, 
such as foams, films, 
powders, or fibers. 

Small sections determined by 
sample pan size. 

mg 

Rheometric Dynamic 
Spectrometer (RDS) 

Quantitative characterization 
of viscous and elastic 
properties 

Measures properties with 
many geometries 

Can handle both solid and 
liquid samples. 

Small sections tha fit into test 
cells. 

mg 

Dynamic Mechanical Analysis 
(DMA) 

Quantitative characterization 
of viscous and elastic 
properties 

T„ and T m . 

Separates materials by 
differences in mechanical 
properties 

solid samples 

small sections that fit into test 
cell. 

mg 

Optical Microscopy 

Qualitative analysis of 
Particulates 

Fast 

Relatively inexpensive 

■HiliiiiHI 

Small samples 

ng 

X-ray Diffraction (XDS) 

compound identification of 
crystalline materials 

Identify crystalline 
compounds 

Slow process 

Small sections of single 
crystal or powder 


Transmission Electron 
Microscope (TEM) 

Microstructural analysis of 
metals, ceramics and 
polymers. 

Very high magnification 

Slow & expensive 

small sections of solids 
Sample preparation can be 
difficult 

Hi 

Scanning Electron Microscope 
(SEM) 

Surface Morphology and high 
spatial resolution of small 
samples. 

1 

* 

Samples must withstand High 
vacuum 

small sections of solids 
Sample preparation can be 
time consuming 
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TABLE 1.6 . COMMON SURFACE ANALYSIS TECHNIQUES 


Method 

Applications 

Advantages 

Method Limitations 

Sample Limitations 


Auger Electron Spectroscopy 
(AES) 

Compositional analysis of 
conducting surfaces 

High Spatial resolution. 
Surface sensitive in upper 10 
nm. 

Insensitive to He and H 
Semi-quantitative 

Solids, must be conductive, 
vacuum compatible 

microns 
down to 
nm. 

X-ray Photoelectron 
Spectroscopy (XPS) 
and Electron Spcetroscopy for 
Chemical Analysis (ESCA) 

Elemental analysis of surfaces 
and coatings 

Chemical state identification 

Nondestructive 
Quantitative and rapid 

Interrogates only top 10 nm. 

Solids, vacuum compatible 

cm 

Ion Scattering Spectroscopy 
(ISS) 

Identification of elements on 
surfaces 

Depth profiling of ultrathin 
films 

Depth profileing can analyze 
for elments through several 
thousand Angstroms. 

Can be time consuming to 
depth profile several thousand 
Angstroms 

Solid , vacuum compatible 

down to 
0.05 cm 

Secondary Ion Mass 
Spectroscopy (SIMS) 

Surface compositional 
analysis 

Trace element analysis of 
surfaces and thin films 

Good resolution from 1 to 5 
nm in depth 

Analysis is destructive 
and can be time consuming 

Solid, vacuum compatible 
Flat surfaces desired 

cm 

Low Energy Electron 
Diffraction (LEED) 

Surface analysis of conducting 
materials 

High resolution 

Surface preparation can be 
expensive and difficult 

Small sections of solids 
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Infrared spectroscopy is usually a good place to start. Due to the variety of 
sampling accessories available, this technique can accommodate most types of samples. 
Infrared provides information about the organic and inorganic functional groups present, 
and an estimate of the sample’s purity. 

Functional groups present in a material can be identified by infrared 
spectroscopy with the aid of correlation charts, such as the ones pictured in Figures 1.29 
and 1.30. These charts summarize the characteristic frequencies of the major functional 
groups. Frequencies associated with carbon-carbon and carbon-hydrogen bonds are 
pictured at the top of each chart. Functional groups containing oxygen and silicon 
atoms are featured in the middle of each chart. Characteristic frequencies of groups 
containing nitrogen, sulfur, phosphorus, and the halogens are shown at the bottom. 
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Figure 1.28. Correlation Chart of Infrared Group Frequencies, 1500-3700 cnr* 
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Figure 1.30. Correlation Chart of Infrared Group Frequencies, 600-2000 cm' 


When cautiously used, correlation charts can be valuable tools in the 
interpretation of infrared spectra. However, correlation charts alone can seldom 
unambiguously establish the molecular structure or identify an unknown material. Group 
frequency ranges are determined empirically from a large number of samples and are 
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often broad. Group frequency ranges may overlap each other considerably. Because the 
physical state of the sample (i.e., the sampling technique) can cause band frequency 
shifts or broadening, the sample should be run under conditions similar to those used 
during compilation of the correlation chart. Throughout the interpretation process, 
always keep in mind all of the information known about the sample, such as its physical 
state, results of other analyses, and elements known to be present or absent 

If the sample is a complex mixture, it may be possible only to identify the major 
functional groups present. If the material is a pure compound, it can be unambiguously 
identified if its spectrum exactly matches the reference library spectrum of a known 
compound. When infrared spectral libraries are not available or no match can be found, 
the compound's molecular structure can sometimes be deduced by interpretation of the 
infrared spectrum with the aid of the correlation charts. The following general 
procedure may be successful: 

I. Look at the group frequency region (4000 to 1350 cm"*) first, 
concentrating in turn, on the strong bands, then on the medium intensity bands. 

II. Look at the 2800 to 3300 cm'l region to determine the presence and 
types of carbon-hydrogen vibrations. 

A. No bands at 2800-3300 indicates no C-H bonds. Consider totally 
halogenated organics, or inorganics. 

B. Bands at 3000 to 3300 cm" 1 indicates presence of unsaturated carbon 
atoms or a halogenated compound. 

C. Bands at 2800 to 3000 cm'l indicate presence of saturated carbon 

atoms. 

III. Now look at the rest of the group frequency region. 

A. Bands at 1450-1465 cm"l indicate presence of methyl (-CH3) or 
methylene (-CH2-) groups 

1. A band at 1375 to 1380 cm"l indicates C-CH3 

2. A band at 718 to 720 cm"* indicates a string of 7 or more 

methylene groups 

B. Bands at 1470 to 1525 cm -1 and 1565 to 1620 cm* 1 indicate presence 
of aromatics. Use Figure 1.29 to determine substitution pattern. 

C. Continue to interpret the strong (and then the medium) intensity bands 
in the group frequency region. Follow-up on each interpretation by examining 
other regions of the spectrum. As examples: 
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1. A band observed at 1715 cm‘l indicates the presence of a 
carbonyl (C=0) group. Is the compound a carboxylic acid.an 
ester, a ketone, or an amide? The presence of two bands near 

3400 and 3520 cm'l along with a strong band near 1600 cm - *, would 

identify the compound as a primary amide. 

2. A strong band is observed at 1060 cm'l . Is the compound an 
ether, a siloxane, or an alcohol? A strong band at 3600 cm‘l 
would identify the compound as an alcohol. 

IV. The following generalizations should be kept in mind: 

A. Not all of the bands can be interpreted as group frequencies. Some 
bands are due to vibration of the molecule as a whole. Other bands are 
due to combinations of fundamental group frequencies. 

B. Aromatic compounds tend to give sharp bands. 

C. Spurious bands may appear (i.e., bands not attributable to the sample). 
For example, bands at 3300-3700 cm -1 and 1600-1 S00 cm -1 may be due 
to water in the sample or in the atmosphere within the sample 
compartment. Similarly, bands near 2325-2350 cm" * and 670 cm - * may 
be caused by the presence of carbon dioxide. 

D. The infrared spectrum alone may not provide enough information to 
identify the compound. Never ignore information from other analyses. 

Figure 1.32 is a guide for the selection of a chromatographic method. 
Selection is based on three sample properties: volatility, complexity, and 

polarity. Speed, resolution, and the quantity of sample should also be 
considered. In general, gas chromatography surpasses liquid chromatography in 
resolution and speed. Capillary columns (or open tubular columns) give the best 
resolution, but packed columns can accommodate larger samples. 
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Figure 1.31. Characteristic Infrared Frequencies for Benezene Derivatives 
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Physical Volatility Complexity Polarity Chromatographic Method 
State 



-gas 


simple 

complex 


GC/PC/SA 

GC/PC/SA. GC/CCVSA 


GC gas chromatography 
PC packed column 
CC capillary column 
NS nonpolar stationary phase 
P polar stationary phase 
SA solid adsorbent 


LC liquid chromatography 
RP reverse phase 
NP normal phase 
IC ion chromatography 


Figure 1.32. Selection of a Chromatographic Method 

The sample’s volatility, complexity, and polarity are often known or can be 
surmised based on information provided by the customer. If necessary, a 
microdistillation can be run to determine volatility. An infrared spectrum will identify 
functional groups and provide information about the sample's polarity and complexity. 

If it is determined that the sample is volatile enough to be run by gas 
chromatography. Table 1.7 and Figure 1.33 can be used as guides for stationary phase 
and detector selection. If the sample is nonvolatile or thermally unstable, an LC 
separation mode can be selected based on the samples solubility, polarity, and molecular 
weight according to Figure 1.34. If the sample is believed to be composed of a wide 


74 



molecular weight range of species, or if it is thought to contain species with molecular 
weights over 2000, size exclusion chromatography should be considered. 


Required Response Selectivity 

Required Detection Limit 

Appropriate Detector 

universal detection (general 
survey) 

10 ppm 

TCD, thermal conductivity 

nearly universal (most organics) 

0.1 ppm 

FID, flame ionozation 

compounds containing notrogen or 
phosphorus 

1 ppb 

TED, thermionic emission or NPD, 
nitrogen-phosphorus 

compounds containing 
electronegative atoms, halogenated 
and sulfur containing compounds, 
anhydrides, nitrates, nitriles, 
peroxides 

1 ppb 

ECD, electron capture 

highly selective for many classes 
of compounds, compound 
identification possible 

1 ppm 
1 ppb 

MS, mass spectrometer 
SIM, mass spectrometer with 
selective monitoring 


Figure 1.33. Detector Selection for Gas Chromatography. 


Types of Separations 

Required Polarity 

Stationary Phase Composition 

hydrocarbons, gases 

nonpolar 

hexamethyltetracosane 

gneral purpose 
(boiling point separation) 

nonpolar 

poly(dimethylsiloxane) 

genral purpose 
aromatics, olefins 

nonpolar 

poly(dimethyl(diphenylsiloxane) 
(80% methyl, 20% ohenyl) 

aromatics, phenols 
glycols 

semipolar 

poly(dimethyldiphenylsiloxane 
(50% methyl, 50% phenyl) 

alcohols, esters, ketones 
carboxylic acids 

polar 

polyethyleneglycol 


Table 1.7. Stationary Phase Selection for Gas Chromatography 
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anion exchange 
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polar 


molecular weight difference > 10% ^ 
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normal-phase 
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water soluble 


“►■gel filtration 


water insoluble 
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Figure 1.34. Selection of a Liquid Chromatographic Mode 

Real samples cannot always be classified as totally volatile or nonvolatile. 
Consider, for example, a polyurethane foam formulation containing extremely volatile 
components, such as CFC or HCFC blowing agents, along with nonvolatile, polymeric 
alcohols and siloxanes. In this type of situation, it is necessary to use both gas and 
liquid chromatographic procedures to isolate each ingredient. 

Once a chromatographic separation has been obtained, the separated components 
can be collected and identified by IR, MS, or NMR. Coupled, or hybrid, instrumental 
techniques such as GC/MS, GC/IR, or LC/MS permit both separation and identification 
of unknown components in a single instrumental procedure. 

If the elemental composition of the sample is of interest, a qualitative survey 
analysis can be obtained by XRF provided the sample is nonvolatile. For quantitative 
elemental analyses, or for volatile samples, an atomic spectroscopic technique can be 
selected based on the elemental detection limits shown in Figure 1.35. In selection of an 
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atomic spectroscopic method, the advantages and limitations of each method (given in 
Table 1.4) should be kept in mind. For example, the lower detection limits provided by 
furnace (electrothermal) atomic absorption are accompanied by decreased precision and 
increased analysis time. Both ICP/AES and I CP/MS are capable of multi-element 
analysis, while only one or two elements can be determined simultaneously by AA. 

As stated before, no single procedure exists for the selection of the best 
instrumental method combination to solve an analytical problem. As an example, 
consider a polyurethane foam material received as two components, one of which is a 
mixture of polyols, flame retardants, catalysts, surfactants, and blowing agent. Some of 
these ingredients are mixtures themselves, and some contain solvents. The ingredients 
differ significantly in concentration, molecular and elemental composition, and physical 
and chemical properties. Because of the material's complexity, a combination of 
instrumental methods is required. FTTR can be used to identify the functional groups 
present, and to monitor the composition of the polyols which are present at a high 
concentration level. GC can be used to monitor the more volatile constituents. Because 
of the wide range of volatilities and concentrations, more than one GC procedure is 
required. ICP/AES or other atomic spectroscopic methods can the be used to measure 
levels of specific elements, such as phosphorus, silicon, potassium, or tin to monitor the 
flame retardants, surfactants, and catalysts present in the material. 
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Figure 1.35. Atomic Spectroscopic Detection Limits (ng/mL). 
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2.0 CHEMOMETRICS 


Chemometrics is a discipline within analytical chemistry concerned with the selection 
and optimization of instrumental methods as well as the interpretation of data from these 
chemical analyses. Chemometrics makes extensive use of mathematical and statistical 
methods with the intent of producing a maximum of concise chemical information. 
Given these mathematical underpinnings, this section presents a description of a number 
of the most useful methods together with a discussion of their application to the 
resolution of problems solved by chemical fingerprinting. 

Chemometrics has undergone rapid evolution and proliferation due, in a large part, to 
the availability of computer hardware and software capable of complex calculations 
performed on large quantities of data. Although excellent statistical software packages 
(programs) running on a variety of computer platforms exist, little or no guidance is 
usually given for the appropriate application of the individual procedure(s). The 
practitioner, therefore, must employ his combined knowledge of chemistry, statistics, 
and the nature and source of the data to ensure that the correct computations are 
performed. It is the purpose of this section to offer some useful advice in these matters. 

2.1 Basic Statistics 

The purpose of any quantitative chemical analysis is to obtain a valid estimate of the 
true value of the chemical characteristic being measured. Variation is ever present and 
inevitably produces a certain amount of uncertainty in the outcome. This uncertainty is 
called error. All variation and the associated errors can be categorized into four main 
groups: 


Common causes - the collection of factors which produce relatively small and 
random changes in results and are sometimes referred to as system errors. While 
these errors can be minimized, they generally cannot be eliminated. 

Special causes- factors which sporadically introduce variation over and above 
inherent system variation. Sometimes called assignable causes because the 
source can usually be discovered and corrected. 

Structural causes - regular systematic changes due to cyclic factors such as 
day/night, moming/afternoon, and seasonal changes. 

Tampering - variation induced by unnecessary adjustments usually made in a 
vain attempt to compensate for inherent (common cause) fluctuations. 

It is imperative that the analyst be able to distinguish between these causes as the 
appropriate remedial measures for each are quite different. When viewed as a process, 
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most analytical procedures can benefit substantially from the application of modem 
statistical process control (SPC) methods. 

When an analytical chemist performs replicate analyses on the same sample using the 
same procedure(s), instrument(s), and reagent(s), the individual results will be found to 
vary from one to another due to common causes (system error). While these results do 
differ, they also exhibit a strong tendency to cluster around a certain value. This value is 
the mean or arithmetic mean, or simply the sample average . Mathematically it is sum of 
the individual results divided by the number of these results and is denoted by "x". This 
sample mean is a statistic and therefore an estimate of the parameter known as the 
population mean which is denoted by "ji". Note that there is a very important 
distinction between the statistic derived from sampled data and a population parameter. 
It is seldom if ever that we have available the complete set of observations for the entire 
population which would allow us to compute the population parameters (ji, s^, and s). 
It is, therefore, necessary to calculate estimates ( statistics ) of these parameters based on 
measurements performed on samples drawn from the population (x, s^, and s). 

[l = S(xi)/N 
x = L(xj)/n 

Where: xj= individual observations 

n = number of observations 
N = total population 

The mean is the most useful statistic. Other statistics include the mode, the median, and 
the midrange. 

The mode is defined as that observation which occurs most frequently in the data set. 

The median is defined as that observation which, when the data are arranged in order of 
magnitude, is the middle value. 

The midrange (MR) is defined as that observation which is halfway between the largest 
and the smallest observation. 


Midrange — (xmax * x min)/2 


These statistics all provide information about the location of the center of the data and 
are also known as measures of location. This information, while useful, is less than 
complete. An adequate description should also include a measure of how much the 
individual observations differ from the chosen measure of location. Collectively these 
quantities are known as measures of dispersion or measures of variation and include the 
range, the mean absolute deviation, the variance, and the standard deviation. 
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The range is the simplest measure of dispersion and is numerically equal to the 
difference between the largest observation and the smallest. 

Range = (x^* - x^/2 

The mean absolute deviation is defined as the mean (average) of the absolute values 
obtained by the subtraction of the arithmetic mean from each observation. It is therefore 
a measure of how much, on average, each observation differs from the mean. 

Mean Absolute Deviation = El(xj - x)\/n 

The variance is a measure of dispersion calculated by summing the squares of the 
differences of the individual observations from the mean, and then computing the mean 
of this sum. It must be pointed out that there are two forms of this measure: the 
population variance which is a parameter and denoted by "s 2 " and the sample variance 
which is a statistic and denoted by "s 2 ". As with the two means above, this is an 
important distinction. 


C 2 =I(xi-p) 2 /N 


s 2 = X(xi- x) 2 /(n-l) 

The standard deviation is defined as the positive square root of the variance and also 
has two forms: the population standard deviation denoted by "s" and the sample 
standard deviation denoted by "s". 


a = [E(xj- •) 2 /N] 1 '' 2 


s = [Z(x i -x) 2 /(n-l)] 1 / 2 


Another term frequently encountered is the coefficient of variation (C.O.V.) or relative 
standard deviation and is defined as the standard deviation expressed as a percentage of 
the mean. 

(population) relative standard deviation = 100 (a/ji) % 

(sample) relative standard deviation = 100 (o/x) % 

Note that measure of relative dispersion cannot be used in cases where the individual 
observations take on both positive and negative values as, for example, after a data 
transformation which has placed the mean at a value of zero. 

A word of explanation about the denominator "(n-1)" in the above equations.This 
quantity (the number of observations minus one) is the degrees of freedom (DF) for the 
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statistic and is, in general, equal to the number of observations minus the number of 
constants calculated from them. Since the mean (x) must first be calculated before the 
variance (s^) can be computed, this "uses up" one degree of freedom and therefore DF 
= (n-1). If "n" rather than "n-1" were used, it can be shown that the variance would be 
biased (underestimated) by a factor of (n-l)/n. It is intuitively obvious that as the 
number of samples increase so does our confidence in the calculated estimates. There 
is, however, a point beyond which additional samples produce a very small increase in 
confidence and this value is approximately 30 samples. 

When one looks closely at the data, there appears to be a pattern of the dispersion of 
observations about the central location. This pattern is called a distribution. A 
distribution can be also be depicted graphically by an "X - Y" plot of the frequencies at 
which given observations occur ("Y" axis) versus the observation values ("X" axis). A 
number of distributions exist and include the normal, the binomial, the Poisson, the chi- 
square, the t, Weibull, and the F. By far the distribution used most often in analytical 
chemistry is the normal or Gaussian. 

The graphical representation of a normal distribution of a large number of observations 
can be accomplished by grouping the observations into classes, each a small and 
consistent range of values, and then counting the number of observations belonging to 
each class. If one next constructs a series of rectangles whose width is the class interval 
and whose height is proportional to the number of observations in that class, then a 
histogram is created by placing the rectangles, in ascending order of class magnitude 
along the "X" axis (left to right) in a contiguous fashion (see Figure 2.1). Note that the 
individual areas of these rectangles are proportional to the density (number) of 
observations in those classes. To make additional observations from the population, the 
probability of obtaining a given value is then determined by the ratio of the area of the 
class rectangle (to which it belongs) to the total area of all the rectangles. Successive 
observations occur most frequently at values within the class intervals of those 
rectangles having the greatest area. 

More observations can be made and placed in smaller and smaller class intervals as 
shown in Figure 2.1(b) - 2.1(c). In the limit, the smooth and continuous bell shaped 
curve of Figure 2.1 (d) is formed. This is the well known shape of the normal 
distribution and is a continuous probability density function symmetrical about a central 
value which is the mean "ji". In fact, the curve (function) is described by the following 
equation: 

y = [1/ o{27t) m ]e 2 ° 2 

Where: n = 3.14159.... 

e = 2.71828.... 
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The above equation contains two constants, "p" and "a", the values of which uniquely 
determine the location of ("X" axis) and the shape (peaked or broad) of the distribution. 
A common chemical example of this distribution is the shape of the chromatographic 
peak. The mathematical term for the shape of this curve is kurtosis\ broad curves (large 
values of "a") are platykurtic while narrow curves (small values of "a") are leptokurtic. 
Note that the curve approaches the "X" axis asymptotically (never reaches zero 
probability) when moving away from the mean in either direction. 



Figure 2. 1 Histograms of Decreasing Cell Size 

A standardized form of the normal distribution, used in many statistical tables, is created 
when the value of is set to zero and the value of "s" is made equal to one. This 
standard normal distribution has been extensively studied and is well documented in 
numerous tables. Normality is a remarkably valid assumption for many distributions 
found in both analytical chemistry and science in general. 

If the standard normal distribution equation is integrated for values of "X" over the 
range of minus infinity to plus infinity, the total area under the curve is equal to one, the 
"X" axis is in units of +/- s(usually denoted by "z" or "z score"), and the maximum is 
located at zero. Any real distribution with a known "p" and "a" can be transformed into 
this standard normal format by this equation: 

z = (X-p)/c 
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This transformation is extremely useful. The area under the curve between any two 
values of "z" is the percentage of the population included within these values. For 
instance, the area under the curve between "p" and either "+la" or "-la" is 0.3413 
representing 34.13% of the population. The area between "-2a" and "+2c" is 0.6826 
representing 68.26% of the total population. Almost all (99.73%) of the population is 
bounded by the interval "-3a" to "+3a". Since these areas also represent probabilities, 
questions like: "What is the probability of obtaining an observation with a value equal 
to or less than "|i" minus "a"? can easily be answered. The answer is easily calculated 
to be 15.87%." This also means that in a sample containing 25 observations we would 
expect that four will have values less than or equal to "p" minus "a". 

The ability to transform sets of observations into their equivalent "z" values provides us 
with a method to scale sets of very different types and numerical ranges for easy 
comparison. This normalizing technique is also useful in the analysis of multivariate 
data. 

Although the mean "x" of a set of observations provides an estimate of the population 
mean "|i", it is very unlikely that "x" is exactly equal to "p". There are two reasons for 
this: (1) the random error in the measurement and (2), the number of observations used 
to calculate "x". If another equal number of observations are performed on a sample 
from the same population, the "x" from this set will probably not equal the "x" from the 
first set Repeating this process would yield a series of "x's" which would, in turn, have 
a normal distribution about a central value - "p". This is the Central Limit Theorem and 
is true even if the population from which the "x's" came is not a normal distribution! 
This distribution is known as the sampling distribution of the mean and has a mean of 
"p" and a standard deviation equal to "a/Cn) 1 / 2 " The term "a/Cn) 1 ^ 2 " is called the 
standard error of the mean and gives a measure of the uncertainty associated with 
estimating "p" from "x". Since, in practice we seldom know the value of "s" and must 
use the estimate "s"; 'V(n) 1 / 2 " is usually not identical to "a/Cn) 1 / 2 ". 

A statistic, "t" has been introduced to compensate for both "s" and the confidence limit 
we place on our estimate of "p". The value of "f" depends on both the confidence 
interval (usually 95%) and the degrees of freedom (DF) for computing "s". We can now 
state the range of values within which we are confident "p" exists as follows: 

p= x * fs/(n) 1/2 

This "f" is itself a distribution having the normal or Gaussian shape. The "t" tables are 
rows and columns for confidence interval and degrees of freedom respectively. As the 
degrees of freedom increase, the kurtosis of the "r" distribution approaches that of the 
normal curve. 
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The "t" statistic affords a means to perform tests on certain conjectures or hypotheses. 
For example: does a method have a systematic error when analyzing a sample of known 
value? In comparing two different analytical methods, do they give results (means) 
which differ significantly or are the differences due to chance alone? The test procedure 
will compare the derived 't' statistic to the table of 'tcntical' values or as usually called, 
the 't table'. If the calculated t for the sample set is higher than the critical t, then the two 
methods do differ significantly. 

Example. The values from the determination of the percentage of a known 
analyte (38.9%) are: 

38.9% 37.4% 37.1% 

Is there any evidence of systematic error? 

Calculations using the various software or calculators gives the values for the 
mean = 37.8% and for the Std. Dev. = 0.964%. Therefore the calculation for t is 
given by: 

t = (38.9 - 37.8) (3/0.964)1/2 = L9 g 

From the "t" Table for DF = 2 and 95% confidence level, the critical value for 
"r" is 4.3. Since the calculated "t" (1.98) is less than the critical "r" value (4.3) 
there is no reason to suspect a systematic error 10. 

There are also instances when, instead of comparing means, we wish to assess the 
precision of two different methods which amounts to a comparison of the variances (sj2 
vs. S2^X Another statistical test known as the "F" test permits this kind of comparison 
by considering the ratio of two variances. By convention the larger of the two variances 
is always the numerator so that the value of "F" is always equal to or greater than one. 
If the differences between variances are small, the value of "F" is close to one and the 
differences are probably due to chance alone. Larger values of "F" imply differences 
too great to be attributed to random causes. The degrees of freedom for the two 
variances need not be equal, but both have an influence on the value of "F" and the 
decisions based upon it. The tables of critical values for "F” are constructed in rows and 
columns of DF of the numerator and denominator respectively and provide for a 
selection of confidence intervals as well. Performing the "F" test involves computing 
the "F" from the variances and for a given confidence interval and DFs, finding the 
critical value from the "F" tables. Once again, if the computed "F" is greater than the 
critical value, real (statistically significant) differences in the variances exist. 
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Example . A proposed method is to be compared to an existing method for the 
determination of chemical oxygen demand in a standard sample. Eight 
determinations by each method produced the following results: 

Mean (mg/1.) Std. Dev. (mg/1.) 

Existing Method 72 3.31 (sj) 

Proposed Method 72 1.51 ($ 2 ) 

F 1J = (si)2/(s 2 ) 2 = (3.31)2/(1.51)2 = 4.8 

From the "F" Table for DFj = DF 2 = 7 and 95% confidence level, the critical 
value for "F" is 3.787. Since the calculated "F" value (4.8) is larger than the 
critical "F" value (3.787), the proposed method is more precise 5 . 

The "F" statistic is also a distribution but, unlike the "r", it is a quadratic function and is 
skewed ("tails out") to the right. 

There arise, from time to time, situations in which one or more of a set of observations 
appear(s) to be quite different from others in the set. The term for this observation is 
outlier. The outlier, upon examination, may be found to be the result of human error 
such as the transposition of numbers or the misplacement of a decimal point. Even after 
the correction of these errors, observations which appear to be outliers may still be 
present. The analyst is now faced with a difficult decision: should these observations 
be retained or rejected r ? The values of the statistics ("x", "s 2 ", and "s") computed from 
the observations will depend on whether or not the outlier(s) are included in the 
statistical calculations. One possible reason for the presence of outliers is that our 
assumption of a normal distribution is not valid. In this event, a test for distribution or 
frequency is appropriate. 

Note that in the Central Limit Theorem that multiple determinations of the sample mean 
"x" from samples of a population has a normal distribution about the population mean 
"ji" with a standard deviation of "o/(n)l/2" which can be estimated by "rs/(n)l/2." What 
about the distributions of the sample variance "s 2 " and the sample standard deviation 
"s"? Like the "F" statistic, the sample variance is a quadratic function and it follows 
what is called the "chi-square" distribution. This distribution is the basis for the " chi- 
square” test which is used to determine if an observed distribution of observations is 
drawn from a population having a certain theoretical distribution ( goodness of fit). The 
test is performed by sorting the observations into classes of observed frequencies; 
calculating the " chi-square " statistic by taking the sum of the squared difference 
between the observed and expected frequencies divided by the expected frequency for 
each class. The degrees of freedom (DF) for this calculation is equal to the number of 
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classes minus one. The calculated " chi-square " value is then compared with the "critical 
value" found in a table at the desired confidence interval and DF. If the calculated "chi- 
square" value is less than the critical value, the theoretical distribution is considered to 
be valid. 

If, however, the assumption of a normal distribution is valid, then there is a test called 
the "Q" test or "Dixon's Q" which can provide some guidance for the accept/reject 
decision. Like the other test statistics, the critical values for "Q" are given in a table. 
The test is performed by computing the "Q" ratio and then from the table finding the 
critical value of "Q" for the desired confidence interval and number of observations. If 
the computed value of "Q" is greater than the critical value, the suspect observation is 
rejected. The "Q" ratio is calculated by dividing the absolute value of the difference 
between the suspect observation and the observation nearest to it by the range of the 
observations (including the suspect value). Note that the rejection of an outlier can have 
a dramatic effect on the values of both the variance and standard deviation. The 
confidence interval implies that there is a small but finite chance of making the wrong 
decision, hence caution is advised in working with these concepts. If more than one 
observation in the set is suspect, the situation is more complex and the observations 
should be repeated or the problem referred to a statistician. 

Most modem analytical chemical instruments can usually make accurate observations 
over a considerable range of values (several orders of magnitude) requiring that the data 
be treated in a manner that is different from the methods employed in the analysis of 
repeated single measurements. The usual procedure involves performing observations 
on samples from populations whose values are known (standards), including a blank (no 
analyte) sample, and span the range expected for the "unknown" samples. From these 
data, a " calibration curve" is prepared and, by interpolation, is used to determine the 
value(s) of the unknown sample(s). Most analytical instruments are designed such that 
the signal is a linear first order function of the sample value (Figure 2.2). Ideally, this 
straight line function should pass through the origin (no analyte, no signal) and have a 
steep slope ( dyldx ) as this is a direct measure of the instrument's sensitivity. This ideal 
relationship can be expressed in this form: 

y = a + bx 

Where : a = the "Y" intercept (ideally = 0) 

b = the slope (ideally > 1) 
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A detailed consideration of this approach raises a number of important statistical issues: 

1. Since this line must be "fitted" to the experimental data, which always 
contains error, what is the "best line" through these data? 

2. What are the error estimates for the constants "a" and "b"? 

3. What are the error estimates for values predicted by the "best line "1 

4. What is the least analyte value that can reasonably be determined (limit of 
detection)? 

5. How valid is the implicit assumption that all errors are in the "Y" values and 
the "X" values are error-free? 

6. If each of the "Y" calibration values is the mean of several observations, are 
the variances of these values equal? 

7. Can a quantitative "figure of merit" be computed to describe how well the "Y" 
calibration values match the values predicted by the equation of the "best line "? 


The most universally accepted answer to Issue #1 is the method of least squares. Using 
the assumption that all error is in the "Y" values (Issue #5), the "best line " is the one 
which minimizes the differences between the observed and predicted "Y". Since these 
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differences (known as residuals) can have both positive and negative values, it is 
rcustomary to attempt to minimize the "sum" of the "squares" of these residuals. Most 
statistical software and many inexpensive pocket calculators can perform this task easily 
and will output values for the constants "a" and "b" as well as a value for the 
" correlation coefficient" (Issue #7). 

It is also possible to calculate a statistic, analogous to the sample variance, for both the 
"a" ("Y" intercept) and the "b" (slope) constants (Issue #2). These statistics are often 
provided by the statistical software also. 

The estimation of the error in obtaining an "X" value from an experimental "Y" value 
(Issue #3) is a mathematically complex undertaking and, if required, should include the 
assistance of an experienced statistician. Formulae for an approximation of this error 
do, however, exist and suggest that the error will be the least when the "Y" values are 
nearest to the "centroid" of the "best line" 

"Limits of Detection" (Issue #4) considerations arise during trace analysis or when the 
"blank” standard produces a measurable instrumental response ("Y" value). There are at 
least two reasons, one chemical and one statistical, for this "blank" response. The 
chemical reason may be due to "interferences" or "matrix effects" and may be 
ameliorated by "spiking” of the samples using the "method of standard additions .". The 
statistical reason, however, is more complex and is a subject of much controversy. A 
reasonable working definition is very method and sample specific and should, therefore, 
be documented in the reported results. 

Our confidence in the assumption of an error free "X" axis (Issue #5) is inexorably tied 
to our confidence in the accuracy of the calibration standards we use. The need for this 
assumption is fundamental to the method and consequently the great care and attention 
given to the preparation, storage, and use of these calibration standards cannot be 
overemphasized. 

If the "Y" values (signals) are the mean of several observations (Issue #6), it is not 
unusual to find that the variances associated with these means increase with increasing 
values of the means. Clearly, the values having the smallest variances should exert a 
greater influence in the determination of the "best line " than those with larger variation. 
It is possible to assign " weights " to the "Y" values which are inversely proportional to 
their variances and use these " weighted values" to compute the "best line." Note that the 
centroid, and along with it, the zone of least error will be moved down the "best line " 
closer to the origin, thereby improving the accuracy of the smaller sample observations 
(Figure 2.3). 
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The " correlation coefficient" (Issue #7) is usually provided along with the values for "a" 
and "b” in the "best line" equation by the statistical software and is denoted by "r". The 
values for this statistic range from -1 to +1 with "r = +1" indicating a perfect point-by- 
point correlation along a "best line" having a positive slope (b > 0). For "r = -1" a 
perfect negative correlation (b < 0) exists. For "r = 0" there is no correlation ("X" and 
"Y" are independent). The square of the " correlation coefficient" (r^) is known as the 
" coefficient of determination" and is a measure of the proportion of the variation of the 
"Y" values accounted for by the "X" values. 

2.2 Design of Experiments 

The process of formulating an accurate, reliable, and efficient plan for the investigation 
of the effects of certain " factors " on the performance of a "system" is the focus of 
Design of Experiments . An analytical chemist involved in "fingerprinting" is faced 
with a task broader than just the measurement of a chemical or physical property of a 
material. The proposed measurements, in many cases, must be evaluated in the context 
of their ability to provide some or all of the information necessary to successfully 
fingerprint the product or material under investigation. Design of Experiments 
methodology is especially valuable for systems with multivariate relationships. Design 
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of Experiments is the planning activity which precedes and interacts with the phases of 
method development, data acquisition, data treatment/processing, and data 
interpretation. Design of Experiments can, and has been, successfully applied to a broad 
spectrum of scientific and engineering investigations. 

The approach to a scientific investigation is composed of answering a number of well- 
phased concise questions. These questions define the analytical problem and are 
formulated on the prior development of a model of the system under investigation. A 
complex system may require more than one model, each of which may constitute an 
independent investigation. The models most commonly used are polynomials of the 
first or second degree which have been demonstrated to be adequate descriptions for the 
vast majority of systems. The model will, in tum, define the requirements for obtaining 
the necessary information including the data and its quality (accuracy, precision, 
repeatability, etc.). 

The Design of Experiments strategy is generally a sequence of plans created to identify, 
analyze, and optimize the effects of factors on a system or analytical method under 
investigation. The first phase is a screening design and will identify those factors (and 
any interactions between them) which have a statistically significant effect on the 
performance of the system. The screening experiments segregate the " vital few " from 
the "insignificant many." The next phase will separate and elucidate the "main effects" 
from the " interactions ." The third phase will provide a mathematical model (see above) 
of the " response surface " for the system and thereby guidance for system optimization 
Figure 2.4 illustrates the methodology one might perform is carrying out such a design. 
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Figure 2.4 Evolution of an Experimental Program 

Although all experimentation requires careful planning, some simple investigations need 
no special considerations (as in the case where only one factor has an effect on the 
system's performance). Design of Experiments, however, offers considerable 
advantages in both the economy of experimental resources and maximum production of 
information especially where there are a number of factors to be investigated. The 
" classical " approach in which the system's response is explored for each factor in turn 
while all other factors are maintained at some constant level is not as efficient nor as 
informative as the "factorial" method of Design of Experiments. The factorial method 
measures the system's response at all or some fraction of the possible combinations of 
the chosen (usually two or three) levels of the factors. There are two convincing reasons 
for the choice of the factorial over the classical method: 

1) The factorial detects and estimates any interaction between factors while 
the classical method cannot. 

2) The factorial requires fewer experiments than the classical method for the 

same precision. 

It can be demonstrated that for "k" factors, a classical approach involves "k” times as 
many experiments as the factorial approach for the same precision. The factorial design 
approach involves measures for "k" factors with "p" levels for a total of "p*" 
experiments. For example, for three factors at two levels each, the factorial method 
requires 8 experiments, while the classical method requires 24 experiments (see Figure 
2 . 6 ). 

In order to understand the factorial design concept, it is necessary to introduce the idea 
of "measurement space." Each of the factors is assigned an axis orthogonal to the rest 
resulting in a " space filling " model encompassing the " experimental volume" of the 
system. For example, a system involving three factors would be represented as a three 
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dimensional cube (Figure 2.5). It is not possible to visualize a cube in more than three 
dimensions (a hypercube), however, the mathematics generalized from three dimensions 
are valid. An additional benefit of factorial designs is the " hidden replication" inherent 
in the geometry of the design. For obvious reasons, it is desirable to maximize the 
"experimental volume" by choosing large ranges (tempered by good judgment) for the 
factors. This is known as "bold" experimentation. 

Factorial experiments are generally performed at either two or three levels of each 
factor. The two level designs are represented as 2^ factorials where "k" is the number of 
factors. These two level designs will not detect any "curvature" in the data and are 
limited to linear first order applications. In many cases a "fractional factorial" design is 
adequate; it 



Figure 2.5 Factorial Design 

is not necessary to perform the full factorial number of experiments. The fractional 
designs require 1/2 to 1/8 the number of experiments of the full factorial design and are 
well suited to the "screening" of potential factors to identify those which are critical to 
the selection and development of an analytical method. 

The three level designs are represented as 3^ factorials where "k" is the number of 
factors. These designs will detect "curvature" in the data and are sometimes called 
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" limited response surface" models. The principle of fractional factorials can be 
extended to the three level designs and require 1/3 to 1/27 the number of experiments of 
the full factorial design. 

A special form of two level fractional factorial design known as "Plackett-Burman" is a 
very economical screening design frequently useful when the number of factors are one 
less than multiples of four (3, 7, 11, 15, 19, 23, 27, 31, etc.). 

The design strategy uses some of the experiments from the screening phase as elements 
of the limited response surface design (by inclusion) and further reduces the number of 
required experiments. This process of "overlaying" the designs is carried on to the 
response surface designs with similar economy. 

The preparation of any of these factorial designs is not a trivial undertaking and is 
usually a task for the statistician. The individual experiments are dictated by the 
geometry of the design and must be performed exactly as specified (even though the 
analyst may think the specified conditions are unrealistic). Clear written instructions 
and close supervision should be provided. The data analysis and interpretation are a 
form of "Analysis of Variance" (ANOVA) and should be performed by, or under the 
supervision of, a statistician. Experience has shown that success with design of 
experiment is only achieved by close cooperation between the anaytical chemist and the 
statistician. Figure 2.6 presents a convenient guide for designs involving up to eighteen 
independent variables. 


18 INDEPENDENT VARIABLES 


LESS THAN 7 VARIABLES 7 OR GREATER VARIABLES 

I I 

INTERACTIONS (7) NO INTERACTIONS 


1 

1 

1 

1 

1 

1 

YES 

NQ 

2 

L£ 

2=12— 

13-18 

2 Level Designs 

12 Ron PB 

12 Ran PB 

12RunPB 

12 Ran 

28 Ran PB 

32 Fnc. 

29 Run PB 

20 Run PB 

28 Run PB 

PB 


Factorial 

16 Run 

16 Ron 

16 Run 

28 Run 



Frac. 

Frac. 

Frac. 

PB 



Frac tonal 

Fractorial 

Fractorial 




Figure 2.6 Selection Guide to Experimental Designs 

It is the intent of design of experiment to provide information for the selection of factor 
values which will maximize the response of the dependent variable. Performing this 
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task is known as "optimization." If no significant interaction between factors exists, this 
task is relatively simple. In the presence of interactions, however, the task can be quite 
complex. As stated earlier, the description of the response surface takes the form of a 
first or second degree polynomial equation in the "k" factors. Visualizing the surface 
when "k" >2 is not possible and solving the polynomial equation for maximum response 
in "k" factors is computationally intense. In those situations where the objective is to 
reach the optimum response and a "model" of the response surface is not required, a 
procedure known as the "Simplex Method" is an excellent and economical approach. 
This "Simplex Method " should not be confused with the "Simplex Technique" of linear 
programming or the "Simplex Mixture " design of formulation problems. A "simplex" is 
a regular geometric figure whose comers (vertices) are all equidistant. For the 
optimization of "k" factors, a simplex of "k +1" vertices is required. Thus for one 
factor, the simplex is a line segment; for two factors, it is an equilateral triangle; and for 
three factors, it is a tetrahedron. 

A useful concept for determining where to start in a Design of Experiments approach to 
improving quality in manufacturing operations has evolved through the Taugchi 
approach. This methodology has become popular in today's quality engineering 
strategies for achieving target values in production and identifying variables which can 
be controlled so as to reduce performance variations. A more extended analysis usually 
requires the use of the approaches shown in Figure 2. .6 to optimally determine the 
responses of factors in an analytical program. 
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2.3 Multivariate Methods 


As the title implies, multivariate methods deal with bodies of data representing 
observations of several (usually more than three) variables on each sample. These data 
may be composed of measurements of very different kinds including both continuous 
and discrete variables and may be expressed by a variety of units encompassing a wide 
range of values. These observations are usually summarized in a data matrix of "p" 
variables (columns) measured on "n" samples (rows). Thus the data can be considered 
to describe "n" points (objects) in "p" dimensions (space). Obviously, if "p" is greater 
than three, visual display of the data becomes a challenge. 

"p" variables 
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"Chemical fingerprinting" is a methodology designed to provide a detailed 
description of complex chemical formulations for which no single measurement 
provides adequate information. It is, therefore, a multivariate approach based on 
measurements developed in designed experiments and addresses four important issues: 

reduction of dimensionality 
multivariate correlations 
multivariate classification 
data summarization. 

It is important to keep in mind that the multivariate methods applied to these issues are 
not always exclusive in that, for example, a method for the reduction of dimensionality 
may also provide useful information on possible classifications. Graphical presentation 
of the results of an analysis is a vital and common feature in all of these methods and 
can often provide the analyst with important new insights. 

It is the aim of all multivariate analyses to gain some insight into the structure and 
information content of these data. The efforts to achieve this aim are generally directed 
to a few broad issues: 

1. Reduction of "dimensionality’’ - can fewer than "p" dimensions (variables) 
describe (graph) the data without significant loss of information? 
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2. Multivariate correlations - are there significant correlations between some of 
the "p" variates? 

3. Multivariate classification/discrimination - are there "natural" groupings in 
the data and can individual samples be "assigned" to one of these groups? 

4. Summarization and presentation of results - what are the "best" methods for 
conveying the information content of the data? 

Prior to actually performing multivariate analysis on experimental data sets, one should 
consider the following issues: 

1. In the beginning, it is not always known exactly which direction the 
analysis should take and a good deal of data exploration may be required. 

2. Given data of many different types and value ranges, should the data be 
"normalized" prior to analysis and, if so, by which means? 

3. If some of the "p" variables are more important than others, a scheme for 
"weighing" the variables must be selected and applied. 

4. In most cases the selection of the "p" variables is based on scientific 
intuition and imposes a limit on how much analysis can be performed. 

5. The nature of the problem and available experimental resources determine 
the ratio of "n" to "p" (ideally 10:1). This places constraints on the extent of 
the analysis. 

6. Quite often, the human inability to visualize graphically more than three 
orthogonal axes prevents the analyst from obtaining a necessary "feel" for 
the data and, hence, a notion of what to do next. 

7. Multivariate methods are computationally intense and the "number 
crunching" capacity of most available computers sets an upper limit on "p", 

"n", or both. 

8. Unlike univariate points on a line, points in "p" space do not have a unique 
linear ordering and the human desire to impose this condition may result in a 
misleading view. 

Close cooperation between the analytical chemist and the statistician can overcome 
most, if not all, of the above difficulties and produce insights and understandings of 
complex multidimensional objects (samples) which would, in the absence of these 
methods, be obscured. 

In the section 2.2 Design of Experiments , the idea of "measurement space" was 
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introduced and this concept is a fundamental component of the multivariate approach. 
The individual objects (samples) can be seen as " points " in a space of "p" dimensions. 
The location of the points (in "p" space) is specified by a set of coordinates which are 
the values of the "p" variables for that object (sample). An equivalent description is that 
the " points " are the terminations of "vectors" whose components are the same "p" 
variables. Now it is possible to understand that "points" which are close to each other in 
this "p" space are also "close" in their measured variables; the " distance ” between them 
being a measure of their " similarity ." If all of the "p" variables are of the same kind 
(concentration, absorbance, mass abundance, wavelength, etc.) there is no difficulty 
with this simple model. If, however, the variables are of a mixture of kinds with quite 
different "values," the simple model may be distorted by the larger " values " and 
consideration should be given to "scaling" the variables by some " normalizing " 
technique. By far the most common normalizing technique is the u z score" method 
whereby the variables are each transformed to have a mean of zero and a variance of 
one. The effect of this transformation is to place all of the variables on an equal basis 
and if this is not the case (all variables equally important), the variables must then be 
individually assigned " weights " by the analytical chemist (here, the statistician can only 
provide advice). These decisions concerning " normalization " and "weighing" of the 
variables are important data " pre-treatment " considerations, but, given the preservation 
of the original " raw data" are not irrevocable. 

For those conditions where "p" is greater than three, visual examination of the vector 
terminations (or points) is not possible. However, the principles of Euclidean Geometry 
are also equally valid for vectors in ”p" dimensions. The true length of the vector is 
equal to the square root of the sum of the squares of its coordinates. The Euclidean 
distance "d" between the terminations of any two vectors is simply the square root of the 
sum of the squares of the difference between their respective "p" coordinates. These 
relationships are summarized for 2, 3, and "p" dimensional space in Figures 2.7, 2.8, and 
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Figure 2.7 Two Dimensional Case (Plane Geometry) 
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Figure 2.8. Three Dimensional Case (Solid Geometry) 
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2.9. Note that the two vectors ("Ri" and "R 2 ") and the distance line "d" in each case 
form a triangle and only requires two dimensions to display the graphical features. This 
holds for the comparison of any two objects (points) in "p" dimensional space. Given a 
set of "n" objects (samples) defined by "p" variables, it is possible to calculate all of the 
n(n-l)/2 pairs of "d" distances between the objects. 
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Figure 2.9 Multi (p) - Dimensional Case (Hyperspace) 

A distance between variables can also be calculated by means of a transformation of the 
correlation matrix. The correlation matrix is a square symmetrical matrix providing for 
each pair of variables, a correlation coefficient "r". The diagonal of the matrix is unity 
for all self-correlations. The most direct transformation is made by subtracting the 
absolute value of the individual correlation coefficients from one (1). If negative 
correlation coefficients are to be considered an indication of dissimilarity, then one 
should use the signed values in place of the absolute values. Having done this we now 
have a pseudo distance table for the variables. 

The Euclidean distance measure discussed above is only one of many possible metrics 
or measuring systems. If each of the "p" variables are normally distributed (which is 
most often the case), then the points in "p" space have a characteristic form or pattern. 
This pattern is dependent upon both the individual variances and the correlations 
between the variables. If, for instance, "p" = 2, "r" (correlation coefficient) = 0 (no 
correlation) and the variances are equal, the locus of points equidistant from the centroid 
of the pattern is a circle. Given some degree of correlation between variables and/or 
unequal variances, the locus becomes an ellipse. With "p" > 2, the locus of points 
equidistant from the centroid is the surface of an ellipsoid ("p" = 3) or hyperellipsoid 
("p” > 3). Note that the major diameter of the ellipse (or ellipsoid) is in the direction of 
greatest variance. The Euclidean metric defines the locus of points equidistant from the 
centroid to be a circle ("p" = 2), the surface of a sphere ("p" = 3), or a hypersphere ("p" 
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> 3) and distances (to the centroid) based on this metric can, in a statistical sense, be 
misleading. A metric which does take into account these statistical considerations 
(unequal variances and correlations) is the Mahalanobis distance. Although this metric 
is often superior to the Euclidean, its computation is more complex and involves matrix 
algebra to find the matrix inverse of the covariance matrix: a job best left to a statistician 
and computer. 

A closer look at the multivariate ideas and techniques having the most relevance to the 
search for chemical "fingerprints" and signatures is now in order. 

Issue #1. Reduction of dimensionality 

Two techniques, Multi-Dimensional Scaling (MDS), and Principal -Component Analysis 
(PCA), have direct applications to the "fingerprinting" problem in that both provide a 
picture in two or three dimensions of the similarities (or differences) between the 
signatures of the samples. These methods not only facilitate interpretation, but simplify 
the mathematics (fewer dimensions or variables to handle). Additional benefits include 
the elimination of redundancies (correlations) in the original data and the possibility that 
the abstract "new" dimensions (variables) may have some "real" physical significance 
(reification). 

The technique of MDS is designed to produce a "map" depicting the relationships 
between the objects (samples) based on a table (see above) of distances between all 
pairs of objects. As noted above, these distances are usually expressed in the Euclidean 
metric, but not exclusively so. The MDS procedure first finds from the table the two 
objects which are maximally distant. These objects are then plotted on the map’s central 
horizontal axis at maximum separation. Next, the object at maximum distance from the 
first two is identified and, by an arbitrary convention, is plotted either above or below 
the central horizontal axis to for a triangle whose sides are the distances between the 
three objects. The remaining objects are then in turn ' fitted " into this triangular 
configuration. An initial " goodness of fit" between these configuration distances and the 
distance table is then calculated and used iteratively to change each object's map 
coordinates until no further improvement in "fit" can be made. At this point the final 
map is drawn. Most MDS software give the user some control of: 

1. Map dimensionality (usually 1, 2, or 3) 

2. Assignment of " weights " for the variables 

3. Choice of "normalizations" for the variables 

4. Choice of the distance metric 

5. The maximum number of allowed iterations 

6. Map " convergence " parameters 


103 



7. 


Map editing - size, position, and rotation 


The method of Principal Component Analysis (PCA) is a transformation of the original 
variables into a completely new set with reduced dimensionality, the principal 
components. While there exist the same number of new principal components as original 
variables, there is also the expectation that the first few principal components (2 or 3) 
will account for most or possibly all of the information content in the original data set. 
PCA is general in its scope, makes no assumptions about the nature of the original 
variables, and uses no mathematical model. The principal components (PCs) generated 
by the procedure have the following properties: 

1 Each PC is a linear combination of the original variables. That is, each is 
equal to a sum of the original variables each having unique coefficients. 

2. For each PC the sum of the squares of the coefficients is unity. 

3. The first PC is that linear combination having the greatest variance. 

4. The succeeding PCs are uncorrelated with previous PCs and contain the 
greatest amount of the remaining variance. 

A consequence of 3 and 4 above is a set of variables (PCs) which are uncorrelated with 
one another and, therefore, mutually orthogonal and are arranged in decreasing order by 
the percentage of the total variance they contain. 

The PCA procedure can be explained in geometric terms. If the original "p" variables 
are all normally distributed, the collection of the "n" points (samples) in "p" space will 
be a hyperellipsoid having a definite centroid located at the coordinates of the common 
mean. Finding the PCs corresponds exactly to finding the principal diameters of the 
hyperellipsoid and placing them in decreasing order of length. The locations of the 
projections of the "n" points onto these diameters as measured from the centroid arg the 
individual PCs. 

One should understand that the PCA method is not independent of the scale of the 
original variables and the application of " weightings ” to some of these variables can 
have a profound effect on the results. If the original variables are of quite different 
types, of different units, or of different scales, consideration should be given to 
normalizing or standardizing these data prior to PCA. If a correlation analysis reveals 
that the original variables are all uncorrelated, or nearly so, there is little to be gained 
from a PCA as the probable result would be a computationally expensive coordinate 
transformation. 

In practice, the entry point for the PCA is the matrix of correlations between the 
original variables (or a variance/covariance matrix). The matrix is a square ("p" x "p") 
symmetrical about the diagonal. Mathematical operations are performed on this matrix 
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and result in "p" eigenvalues together with "p" eigenvectors whose "p" components are 
the coefficients described in 1 and 2 above. The percentage each eigenvalue has of the 
sum of all the eigenvalues is the percentage of the total variance it contains. 

Issue #2 - Multivariate Correlations 

In the section Basic Statistics we discussed briefly the calibration of an analytical 
instrument using the method of least squares. That is we developed the equation: 

y = a + bx 

where: a = the "y" intercept, b = the slope 

This is an example of a univariate linear regression. We also calculated a "figure of 
merit" for this equation called the correlation coefficient "r" with values ranging from -1 
through 0 to +1. Recall that large absolute values of "r" indicate a close relationship 
(correlation) between "x" and "y”. 

In the section Design of Experiments an extension of this idea was presented when the 
"fitting" of a "model" was used to define the response surface. In fact, using a first 
degree polynomial equation in the "k" factors to predict the value of the dependent 
variable (response) is an example of a Multiple Regression Analysis. Multiple 
regression analysis is a very useful and general statistical procedure for the study of the 
relationship (correlation) of a single variable (response) to a linear combination of two 
or more other variables (predictors), including a constant term and a random error term. 

yi = bo + bix! + b2X2 + b3X3 + b4X4 + b p x p + e p 

Where:yi = response variable 

bq = constant term 
x l— x p = predictor variables 
bj...b p = regression coefficients 
ep = random error term (residual) 

An example of a chemometric application of this method is the study of possible 
relationship of a measured physical property of a material (tensile strength, viscosity, 
hardness, etc.), to a group of chemical measurements made on the same sample of the 
material. 

As with the univariate regression above, a "figure of merit" can be computed for the 
result of a multiple regression analysis and it is known as the coefficient of 
determination ("r^") and its values range from 0 to 1. The exploration of these 
correlations is complicated by the large number of possible linear combinations of the 
predictor variables and many statistical software programs provide a means to automate 
the search for an optimum combination. There are usually three approaches to the 
problem and all are based on a user supplied criterion (rule) for the determination of the 
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optimum combination: 

1 . A forward selection of predictors. 

2. A backward elimination of predictors. 

3. A stepwise regression (a mixture of #1 and #2). 

Having completed these searches, it is left to the user to make the important decision on 
which linear combination of predictors and coefficients will become the model. The 
software can provide some assistance in the form of "F" statistic values and sometimes 
graphical plots of the predicted values versus the measured values for each combination 
can be examined. 

Example: Determination of the component ratio for polyurethane foam insulation for an 
aerospace system. 

An important requirement in the preparation of a cured thermal protection polyurethane 
foam material is the control of the two part component ratio A and B. A study to 
determine the true relative amounts of the 'A' Component (isocyanate) and the 'B' 
Component (polyol)was conducted using instrumental analytical methods (FTIR, ICP, 
and TGA) to determine the true ratio of these components after cure has been 
completed. The application of these analytical methods together with the use of 
multivariate linear regression techniques permit the determination of the true 'A':'B' ratio 
with an error of less than ten percent (10%). 

Fingerprinting work performed on this foam material has demonstrated that FTIR, ICP, 
and TGA are the instrumental methods best suited to the ratio determination. 

A single lot of foam material for which the results from the full complement of 
acceptance tests were available was used to prepare accurately weighed "cup" samples 
(in triplicate) of 0.60:1.00, 0.80:2.00, 1.00:1.00, 1.20:1.00 and 1.40:1.00 'A’:’B' ratios. A 
complement of "production sprayed" samples were also prepared for correlation 
analysis. The ability to produce these samples (cured foam material) from such extreme 
ratios is a testimonial to the robust mature of this formulation. The individual samples 
were then each subjected to FTIR, ICP, and TGA analyses. Both a forward and 
backward stepwise multiple linear regression analysis was performed on the results from 
each method and on a composite of all of the results. The resulting predicting equation 
are of the form: 


Z = a + bjxj + b 2 X 2 . . . + bjXj 

where: Z = A':'B' ratio 

a = intercept (forced to zero) 
b = regression coefficient 
x = predictor variable 
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From the analysis of fourteen wavelengths in the inffa-red spectrum, 1510 cm _ l and 
1067 cm* * were identified as the critical predictor variables. These two intensities 
remained the critical predictors in the composite regression analysis. Partial ”f ’ tests and 
"t" tests confirmed the significance of these tow predictors at the 95% confidence level. 
The final prediction equation defines a three dimensional plane (See Figure 1.): 

Z = 0.02349xi - 0.02494x2 

where: xj = 1510 cm -1 

X2 = 1067 cm"l 



Figure 2.10 Response surface for two wavelengths used in this study. 

This equation produces an adjusted coefficient of determination r^ = 0.999 and 
an interval error estimate of *0.071 meaning that 99.9% of the 'A'r'B' ratio is explained 
over the entire ratio interval with a probable error no greater than 7%. The ICP and 
TGA regressions each had a coefficient of determination of approximately 0.95 and 
interval error estimates of 0.10 (10%). A graphical representation of the response surface 
is shown below. 

Conclusions and Recommendations 
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An on-going study to determine if foam aging has a significant effect on the analytical 
results is being completed. Consideration will also be given to the expansion of this 
effort to include the correlation of the foam's physical and mechanical properties to both 
the ’AVB’ ratio and its analytical predictors. Alternatively, the knowledge gained in this 
study could be applied to the adjustment of some or all of the major material or process 
factors to produce the desired ’AVB; ratio. It is now possible to accurately determine 
whether, in the foam application process, the desired 'A'r’B' ratio has, in fact, been 
achieved. FTIR analysis, in particular, revealed that the ratio of infra-red light 
absorption be a pulverized sample of foam at 1410 cm** and 1067 cm _ l is a direct 
measure of the relative amounts of the A' and 'B' components reacted to form the cured 
foam. 


Example. Using factor analysis to determine the reliability of a spectroscopic analysis, 
by determining the number of significant factors 

Consider the case where a spectroscopic analysis is performed in the presence of other 
potential absorbants and system noise. The analyst records measurements for two peaks 
in the absorption spectrum of the desired compound at different concentrations . This 
approach can be used to answer the following questions: 

(a) How many compounds are in the mixture? 

(b) What are the intensities at these wavelengths for these compounds ? 

(c) What is the concentration of each compound in each experiment? 

The theoretical absorption for peak 2 should be twice that of peak 1; however, the data 
is recorded as: 


Wavelength 1 Wavelength 2 


0.8404 

2.0296 

1.8982 

4.1676 

3.1831 

5.8166 

4.1665 

7.9562 

5.3620 

9.8650 

5.9200 

11.6624 

7.0987 

13.8794 

8.1588 

16.3139 

8.9828 

18.0685 

10.3643 

20.3791 


The data can be graphed according to the figure below. Note that the linear fit is not 
very good. Variations from the linear relationship may be due to system noise or the 
effects of unknown absorbants. 

Ordinarily the analyst would determine the best linear fit to the data by performing a 
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linear regression on this data. One would then obtain the equation: 


II = 0.5063I 2 +0.02125. 

This equation is represented by the line drawn through the data points in the graph . In 
most cases this approach is satisfactory; however, there is always a lingering question 
about if this one fit is sufficient. 

Another useful approach is to reduce the dimensionality of the data using principal 
components and see how many factors are needed to accurately model the data. A 
regression and principal component analysis using a commercial software package 
provides the following information: 

Whole-Model Test 
Analysis of Variance 

Source DF Sum of Squares Mean Square F Ratio 

Model 1 87.016642 87.0166 2101.656 

Error 8 0.331231 0.0414 

C Total 9 87.347873 0.0000 

Note that the sum of squares for the one degree of freedom model is much larger than 
the error variance., hence the F ratio is very high. Obviously the analyst is able to apply 
the least squares fit to the data with confidence. 

Further results from the analysis gives: 

Response: Wavelength 1 

Summary of Fit Rsquare 0.996208 

Root Mean Square Error 0.203479 
Mean of Response 5.59748 
Observations (or Sum Wgts) 10 

The response parameters show an r^ = 0.996 for the 10 data points. Hence the linear fit 
is good. The response curve for this data set also shows a good fit in the plot: 
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Figure 2. 1 1 Regression fit for spectrometric calibration. 

For the principal components portion of the analysis, one then decomposes the original 
data into principal components of the form: 

Xj = xjcos f + X 2 sin f 
X 2 = -xjsin f + X 2 cos f. 

This analysis reveals the following information: 

Eigenvalue 1.9981 

Percent: 99.9051 

which indicates that the principal component at wavelength 1 accounts for 99.9 % of the 
information, thereby confirming the earlier observations. 


Issue #3 - Classification and Discrimination 


The method of Cluster Analysis is directed to the task of identifying groups of similar 
objects based on the observational data. This analysis is usually performed to search for 
the existence of some natural grouping in the data. The method is relatively simple, 
straight forward, and produces an easily understood dendrogram or hierarchical tree. 
Cluster analysis, like MDS, is based on distances and can be applied to the study of the 
relationships between individuals (samples) as well as variables. 

There are two basic approaches to the clustering problem: 
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1 . The divisive approach which begins with all of the objects in a single 
cluster and proceeds to subdivide this cluster into separate clusters based 
on dissimilarities. 

2. The agglomerative approach which begins with each object in its own 
cluster and proceeds to form other clusters based on similarities. 

By far, the agglomerative approach is the most popular and, therefore, the one we will 
examine. The method requires measures of similarity between all objects being studied 
and this is usually a distance matrix. Any data which can be used to create such a 
matrix is, therefore, suitable and the previous comments on the distance metric apply. 
The four general steps in the clustering procedure are: 

1. Assign each object to its own unique cluster. 

2. Locate the shortest distance in the matrix and merge these two clusters. 

3. Update the distance matrix for this cluster's distance to the remaining 
clusters. 

4. Repeat steps #2 and #3 until there is only one cluster. 

During this procedure a linked list is generated which contains, in order, the mergings 
(fusions) and the distance at which they occurred. This list is then used to draw the 
dendrogram or tree. The rules by which the distance matrix is updated in step #3 above 
are known as the clustering algorithm and there are five popular ones: 

1. The single linkage algorithm 

2. The complete linkage algorithm 

3. The average linkage algorithm 

4. The centroid algorithm 

5. Ward's method 

Implementation of these algorithms on a single data set will generally produce five 
distinct dendrograms and this raises the question "Which is best?" The best will usually 
be the dendrogram with the most faithful reproduction of the original distance matrix. 
The distance between any two clusters (objects or groups) is the distance level at which 
they both first appear in the same cluster (fusion level). With this understanding, a 
distance matrix can be constructed from the dendrogram and then compare the 
corresponding entries to the ones in the original distance matrix and even calculate a 
goodness of fit in much the same (RMS Error) way as for the MDS map. In fact both 
cluster analysis and MDS constitute summaries (in different visual forms) of the 
information contained in a distance matrix. 

The analyst must make many of the same basic decisions prior to a cluster analysis as 
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was necessary for MDS, namely: 

1. Should the variables be normalized? 

2 Should "weights" be assigned to the variables? 

3. Should principal components be used and if so, how many? 

4. Which algorithm should be used? 

5. If variables instead of samples are being studied, how should the distance 

matrix be calculated? 

6. Which "goodness of fit" calculation should be used and what is the 
decision value ? 

A major criticism of cluster analysis is that it will by its very design always find clusters 
even when "homogeneous" or "clusterless" data are analyzed. Any clusters identified 
by this method should be confirmed by MDS and scatter plots. The analyst should be 
aware that the results of cluster analysis depend a great deal on the distance matrix and 
the particular algorithm used. 

The method of Discriminant Analysis deals with the problem of developing an optimal 
separation rule to distinguish between two or more known ( a priori ) groups of objects 
(samples) based on the measurements of several variables on all of the objects. The 
method could be used, for example, to investigate how well samples of acceptable and 
unacceptable materials can be separated using a number of chemical measurements 
applied to both. The intent is to develop a rule which would, on the basis of the 
chemical measurements, correctly allocate a subsequent sample to one or the other 
group (acceptable or unacceptable). While the advantages of such a rule are great, there 
is a certain probability of misallocation and the rule must be chosen to minimize this 
risk. 

There are a number of ways to view this task. First, it can be seen as an extension to 
multivariate observations of the method of analysis of variance wherein the "between 
group" variances are maximized while at the same time the "within group" variances are 
minimized. Like multiple regression analysis, discriminant analysis is concerned with 
finding an optimal linear combination of the original variables (or a subset thereof) 
which satisfies these variance objectives. That is to say, the optimal linear combination 
(an equation having a constant term and coefficients for all included variables) would, 
when applied to all objects in both groups, yield means for the groups which have the 
greatest possible numerical difference and variances for each group that are as small as 
possible. Application of this rule would then create a condition where the centers 
(means) of the groups are far apart and there is minimal overlap between their 
distribution curves (minimum variances). The calculated numerical result for any 
individual object (sample) is known as a discriminant score and the equation's individual 
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coefficients can be viewed as weights for their respective variables. 

DF = bo + bixi + b2X2 + b3X3 + b4X4 + bpXp 

WhererDF = discriminant function 
bg = constant term 
xj.-.Xp = variables 
b^...bp = discriminant coefficients 

An equivalent geometrical interpretation is to find the centroids of the hyperellipsoids of 
the groups in "p" (or less) space which are maximally distant from one another. It is 
preferred that the distance be the Mahalonobis distance based on the normalized (z- 
scored) variables. Then the degree to which the hyperellipsoids interpenetrate each 
other sets the limit on the method's ability to properly allocate an object to a given 
group. Having located the group centroids and objects in "p" (or less) space, the simple 
method of "k" nearest neighbors can be used to assign a new object to the proper group. 
The method assigns the new object to the group having the largest percentage of the "k" 
neighbors closest to it. The value of "k" is usually taken to be a small odd number in 
order to prevent "tie" votes. Given a rather large number of assigned objects in the 
groups, this method is rapid, straight forward, and makes good statistical sense. The 
other methods of making assignments require the calculation of discriminant score value 
known as the cutting score and this can be a complicated matter if the population size of 
the groups differ significantly. 

Since the mechanics of discriminant analysis and multiple regression analysis have 
mathematically much in common, the statistical software programs usually provide the 
same types of options for exploring the data: "forward" selection, "backward" 

elimination, and "stepwise" operations on the original variables. Here too, it is 
sometimes possible to examine a graphical plot of the separation results. 

Issue #4 - Presentation of Results 

While it cannot be denied that tables and listings of calculations and results efficiently 
present the essential hard facts, graphical displays or other visual representations of 
multivariate data are extremely useful tools for their examination and presentation. 
Contour maps (response surfaces), similarity maps (multidimensional scaling and 
principal components analysis), dendrograms/trees (cluster analysis), and scatter plots 
(regression analysis) are all intuitive and easily understood by everyone. A deep 
understanding of the mechanics behind the generation of these visualizations is not a 
prerequisite to reading the information they contain and this should be a major 
consideration when presenting results to the non-statistician. Frequently a picture is not 
only the best , but the only way to convey the complex relationships between 
multivariate objects. Beyond their value in data presentation, the graphical techniques 
are indispensable in the exploring and formulating stages of analysis. A great deal of 
good and practical advice on the design and preparation of a variety of charts and graphs 
is contained in the book by W. S. Cleveland and the two books by E. R. Tufte (see 
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Reference Materials). 

Final Comments and Advice 

With most multivariate data sets there are almost always a variety of alternative 
analytical approaches. The choice of which method is most appropriate depends on the 
data type(s) and the objective(s) of the analysis. Unfortunately, no one method or 
technique is necessarily the best choice and it is often wise to use several methods to 
explore different facets of the data. There are two wavs to view the data set: as the 
relationships between the individuals (samples) as defined by the variables and their 
values, and as the relationships between variables defined either by their pairwise 
correlations or their variance/covariance structure. In the first case we examine the "n" 
individuals in the "p" variables space while in the second we see the "p" variables in ”n" 
space. 

Before undertaking any multivariate analysis, the univariate summary statistics for each 
variable should be thoroughly studied. The mean , variance , standard deviation, range, 
the skew, the kurtosis, and the z-scores can all be easily calculated using available 
software. This variable by variable "quick look" should be supplemented by histograms 
to verify distribution assumptions and to screen for "outliers." Next, the correlation 
matrix should be created and examined, followed by a scrutinizing of the scatter plots of 
all pairs of variables. The correlation matrix contains only the simple linear correlations 
and if the relationships are more complex, the scatter plots together with "eye-brain" 
device may uncover them. Recall that if the variables are all relatively uncorrelated, a 
principal components analysis is probably pointless. The degree to which some of the 
variables are correlated will give some guidance in the investigation of possible linear 
combinations. At this point it is advisable to consider the question of normalization of 
the variables and whether some of the variables should be assigned weights. 

A statistical computer software program will obligingly perform any of its analyses on a 
data set whether or not the method is appropriate. Consider the problem to be solved 
and your objectives, then use common sense (the computer has none). 

Just as in the univariate case, multivariate outliers can exist and if scanning the data by 
eye fails to detect them, consider using the Mahalanobis distances from the common 
centroid as a discriminator. 

Finally, recognize that multivariate analyses do not always give " text book." answers, 
even after a great deal of exploration. It is also very difficult to know when one has 
exhausted all of the many possible approaches. One of the best ways to develop skills 
and gain experience in these methods is to use available software to analyze the 
classical (known) data sets (many of the data sets are included in the books listed in the 
Reference Materials). 
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AAS 

Atomic absorption spectroscopy 

AES 

(1) Atomic emission spectroscopy, (2) Auger electron spectroscopy 

AFS 

Atomic fluorescence spectroscopy 

AA 

Atomic absorption 

AI 

Artificial Intelligence 

AMU 

Atomic mass unit 

ANOVA 

Analysis of variance 

ASCII 

American Standard Code for Information Interchange 

ESCA 

Electron spectroscopy for chemical analysis 

ATR 

Attenuated total reflectance 

AUFS 

Absorbance units full scale 

CC 

Capillary column 

CFC 

Chlorofluorocarbons 

Cl 

Chemical ionization 

CILO 

Computer Integrated Laboratory Operations 

COV 

Coefficient of variation 

CPU 

Central processing unit 

CRT 

Cathode ray tube 

DF 

Degrees of freedom 

DOE 

Design of experiments 

DSC 

Differential scanning calorimetry 

DTGS 

Deuterated triglycine sulfate 

ECD 

Electron capture detectors 

EDS 

Energy dispersive X-ray Spectroscopy 

El 

Electron impact 

ES 

Expert Systems 

ET 

External Tank 

FAAS 

Flame atomic absorption spectroscopy 

FAB 

Fast atom bombardment 

FES 

Flame emission spectroscopy 

FID 

Flame ionization detector 

FTIR 

Fourier transform infrared spectroscopy 

GC 

Gas chromatograph 

GC/FID 

Gas chromatographic with a flame ionization detector 

GC/MS 

Gas chromatography/mass spectrometry 

GLC 

Gas-liquid chromatography 

GSC 

Gas-solid chromatography 

HCFC 

Hydrochlorofluorocarbons 
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High performance liquid chromatography 


IC 

Ion chromatography 


ICP 

Inductively coupled plasma 


ICP/AES Inductively coupled plasma/atomic emission spectroscopy 


IEC 

Ion-exchange chromatography 

— 

IR 

Infrared spectroscopy 


IRE 

Internal reflection element 
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Ion scattering spectroscopy 

— 

LAN 

Local area network 


LC 

Liquid chromatography 


LLC 

Liquid-liquid chromatography 

— 

LIMS 

Laboratory Information Management System 


LHDS 

Laboratory Host Data System 


LNTB 

Laboratory Network Test Bed 


LSC 

Liquid-solid chromatography 


MAD 

Mean absolute deviation 


MAPTIS 

Materials and Processes Technical Information System 


MCT 

Mercury cadmium telluride 


MR 

Midrange 

— 

MS 

Mass spectrometer 


NMR 

Nuclear magnetic resonance 


NPD 

Nitrogen-phosphorus detector 

— 

NP 

Normal phase (LC term) 


NS 

Nonpolar stationary phase (GC term) 


P 

Polar stationary phase (GC term) 


PC 

(1) Personal computers; (2) Packed column (GC term) 
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Photo multiplier tube 


QC 

Quality control 


QEL 

Quality Evaluation Laboratories 
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Quality Laboratory Instructions 

— 
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Quadrupole mass analyzer 


RDS 

Rheometrics dynamic spectrometry 


RI 
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— 

RP 

Reverse phase (LC term) 


SA 

Solid adsorbent (GC term) 
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Scanning Auger microscopy 

— 
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SFC 

Supercritical fluid chromatography 


SIMS 

Secondary ion mass spectroscopy 
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Selective ion monitoring 
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Standard operating procedures 
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Statistical process control 


SRM 

Standard Reference Materials 

— 
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Thermal conductivity detector 
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TGA Thermo-gravimetric analysis 

TMA Thermal Mechanical analysis 

UV Ultra-violet 

XPS X-ray photoelectron spectroscopy 

XRF X-ray fluorescence spectroscopy 



